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Chapter  I 
INTRODUCTION 

In-tube  condensation  occurs  in  a  variety  of  engineering  situations 
such  as  airconditioning  and'  refrigeration,  petrochemical,  process  and 
numerous  other  industries.   The  use  of  more  efficient  heat  exchangers 
in  the  fields  of  nuclear  and  conventional  power  plants,  space  ships,  and 
naval  ships  has  become  a  necessity  in  the  modern  world.   Economic  and 
energy  considerations  have  motivated  the  interest  of  different  investi- 
gators to  explore  the  use  of  various  techniques  to  augment  heat  transfer 
in  heat  exchangers.   The  goals  of  these  investigators  were  either  to  reduce 
the  size  of  heat  exchangers  for  a  given  heat  load,  thereby  reducing  their 
cost,  or  to  increase  the  heat  capacity  of  existing  heat  exchangers.   In 
a  report  published  in  1975  by  Bergles  [1],  it  was  estimated  that  the  in- 
vestment in  heat  exchanger  equipment  in  the  U.S.  was  approaching  one  bil- 
lion dollars  yearly.   Ten  to  twenty  percent  reduction  in  capital  cost 
could  produce  savings  in  the  100  million  dollar  range. 

Usually,  in-tube  condensation  results  in  high  heat  transfer  coef- 
ficients.  However,  fluids,  such  as  certain  organic  liquids  as  well  as 
fluro-carbon  refrigerants,  have  thermophysical  and  transport  properties 
which  usually  result  in  low  condensation  heat  transfer  coefficients. 
Augmenting  the  condensation  heat  transfer  of  these  fluids  could  signifi- 
cantly improve  the  efficiency  and  performance  of  condensers. 

During  the  past  several  decades  numerous  techniques  have  been  tried 
to  improve  the  heat  transfer  in  heat  exchangers.   These  techniques  can 
be  broadly  classified  into  active  and  passive  techniques.   Vibration  of 


2 
the  heat  transfer  surface  and  the  application  of  electrostatic  fields  are 

among  the  active  techniques .   As  active  techniques  of  augmenting  the  heat 

transfer  required  external  power  application,  they  received  little  attention 

by  investigators.   Roughened,  corrugated,  fluted,  and  internally  finned 

tubes,  and  tubes  with  vortex  generator  inserts  such  as  twisted  tapes  and 

static  in-line  mixers  are  a  few  among  the  passive  techniques  for  augmenting 

heat  transfer.  Most  investigations  in  the  past  were  directed  toward  passive 

techniques  because  they  do  not  require  external  power.   These  passive 

techniques  are  easy  to  incorporate,  relatively  simple  to  operate  and 

pose  less  problems  in  maintenance.   Commercial  production  of  most  of  these 

passive  augmented  tubes  are  currently  in  vogue.   One  of  the  factors  which 

is  limiting  the  use  of  some  of  these  techniques  is  the  lack  of  reliable 

design  correlations  for  predicting  the  heat  transfer  and  pressure  drop. 

Considerable  amount  of  testing  is  needed  before  most  of  these  techniques 

can  be  put  into  full  scale  industrial  use. 

The  main  concern  of  the  present  study  was  to  investigate  the  enhance- 
ment of  in-tube  condensation  of  R-ll  by  internally  finned  tubes.   It  was 
an  extension  of  the  investigation  of  Said  [2]  in  which  he  studied  the 
augmentation  of  in-tube  condensation  of  R-113  by  internally  finned  tubes. 
The  tubes  tested  by  Said  [2]  were  also  tested  in  the  present  study  using 
R-ll  test  fluid. 

The  objectives  of  the  present  study  were: 

1.  To  take  different  sets  of  heat  transfer  and  pressure  drop  data  for 
condensation  of  R-ll  inside  internally  finned  tubes  at  different 
flow  parameters. 

2.  To  compare  the  experimental  results  with  the  predictions  of  existing 
correlations  in  order  to  identify  the  correlation  which  best  predicts 
the  experimental  results  of  the  present  study. 


Chapter  II 
LITERATURE  SURVEY 

A  bibliography  of  world  literature  on  augmentation  was  published  by 
Bergles  et  al.  [3].   Over  1900  references  were  cited  in  this  bibliography. 
The  references  were  grouped  under  two  major  classifications: 

1.  Techniques  of  Augmentation:   Active,  requiring  external  power,  or 
passive  requiring  no  external  power. 

2.  Modes  of  heat  transfer. 

Considering  the  vast  scope  of  augmentation,  the  present  literature 
survey  had  to  be  limited  to  in-tube  condensation  (mode  of  heat  transfer) 
inside  internally  finned  tubes  (technique  of  augmentation).  A  prelimi- 
nary review  of  the  subject  matter  revealed  that  design  correlations  of 
internally  finned  tubes  for  condensation  heat  transfer  and  pressure  drop 
are  only  the  extensions  of  smooth  tube  correlations.   Also,  the  develop- 
ment of  design  correlations  for  condensation  is  closely  related  to  single- 
phase  flow  inside  tubes.   So,  a  few  single-phase  flow  studies  inside  tubes 
were  also  included  in  the  present  study. 

As  a  first  step,  the  literature  surveys  of  recent  studies  of  different 
investigators  such  as  Royal  [4],  Luu  [5]  and  Said  [2]  were  reviewed.   Royal 
[4]  studied  experimentally  the  augmentation  of  in-tube  condensation  of 
steam  with  internally  finned  tubes  and  twisted  tape  inserts.   Luu  [5] 
and  Said  [2]  experimentally  studied  the  augmentation  of  condensation  heat 
transfer  of  R-113  by  internally  finned  tubes  and  twisted  tape  inserts. 
Luu  [5]  also  conducted  experiments  with  rough  tubes. 
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Royal  [4]  surveyed  the  condensation  literature  in  general  and  aug- 
mentation of  condensation  inside  tubes  in  particular.   He  classified  the 
general  condensation  literature  into  two  categories,  dropwise  and  filmwise 
condensation.   Filmwise  condensation  is  the  most  common  in  condensation 
processes.   He  classified  filmwise  condensation  literature  into  conden- 
sation outside  tubes  and  condensation  inside  tubes.   The  literature  on 
augmentation  of  inside  film  condensation  was  again  classified  according 
to  augmentation  techniques  into  passive  and  active  techniques.   Passive 
augmentation  techniques  cited  were  surface  roughness,  extended  surfaces, 
and  vortex  generators.   Also,  he  categorized  the  references  according 
to  flow  conditions,  such  as  laminar  or  turbulent,  and  according  to  whether 
the  design  correlations  developed  were  experimental  or  analytical.   He 
cited  over  40  references  dealing  with  predicting  in-tube  condensation 
heat  transfer  coefficients  and  over  70  references  on  augmentation  of  con- 
densation.  Royal's  [4]  literature  survey  was  a  comprehensive  attempt  at 
reviewing  the  condensation  literature  available  prior  to  1975. 

Luu  [5]  studied  the  literature  systematically  according  to  subject 
matter.   The  broad  headings  under  which  he  reviewed  the  literature  were: 

1.  Flow  regime  studies  in  two-phase  flow 

2.  Heat  transfer  studies  of  single  and  two-phase  flow 

3.  Augmentation  of  heat  transfer  using  twisted  tapes,  internally  finned 
tubes,  and  rough  tubes. 

In  the  heat  transfer  studies  of  two-phase  flow,  Luu  included  references 
pertaining  to  vertical  film  condensation  and  in-tube  condensation.   In-tube 
condensation  literature  was  then  critically  reviewed  to  include  various 
parametric  effects  on  heat  transfer  and  pressure  drop,  such  as  inlet 
superheat,  heat  flux,  tube  inclination,  pressure  level,  diameter  to  length 
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ratio  of  condenser,  and  fluid  properties.   The  literature  of  augmentation 

heat  transfer  using  internally  finned  tubes  included  single  and  two-phase 

flow.   Luu's  [5]  literature  survey  cited  over  200  references  published 

prior  to  1980. 

Said  [2]  attempted  to  update  the  condensation  literature  incorporating 
work  done  up  to  1982.   Omitting  some  of  the  basic  references  which  were 
already  reviewed  thoroughly  by  Royal  [4]  and  Luu  [5],  he  included  the 
references  relevant  to  his  study. 

Considering  the  objectives  of  the  present  study  and  realizing  the 
extent  to  which  the  literature  has  been  reviewed,  an  attempt  to  review 
the  widely  used  design  correlations  for  heat  transfer  and  pressure  drop 
during  condensation  inside  smooth  and  internally  finned  tubes  was  made. 
More  interested  readers  are  urged  to  refer  to  references  [2,3,4,5]. 

SMOOTH  TUBES 

Several  correlations  are  available  in  the  literature  for  predicting 
local  and  overall  average  condensation  heat  transfer  coefficients  inside 
smooth  tubes.   In  general,  these  correlations  were  either  based  on  single- 
phase  flow  heat  transfer  correlations  or  on  analogy  between  heat  and 
momentum  transfer  in  the  condensate  film.   The  correlations  of  Akers 
et  al.  [6],  Boyko  and  Kruzhilin  [7],  Rosson  and  Myers  [8],  Traviss  et  al. 
[9],  Soliman  et  al.  [10],  Azer  et  al.  [11],  Cavallini  and  Zecchin  [12], 
and  Shah  [13]  are  a  few  of  the  several  correlations  available  for  pre- 
dicting heat  transfer  coefficients  during  in-tube  condensation.   Reference 
to  the  correlations  of  Akers  et  al.  [6],  Boyko  and  Kruzhilin  [7],  and 
Shah  [13]  will  be  made  in  a  later  chapter. 

Pressure  drop  during  two-phase  flow  inside  a  horizontal  tube  is  a 
combination  of  frictional  and  momentum  parts.   The  correlations  of 
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Lockhart-Martinelli  [14],  Dukler  et  al.  [15,16],  and  Hughmark  [17],  are 

widely  used  for  the  estimation  of  the  frictional  pressure  drop.   The  mo- 
mentum pressure  drop  requires  the  knowledge  of  the  void  fraction  variation 
in  the  direction  of  flow.   The  void  fraction  correlations  of  the  homogeneous 
model  [18],  Lockhart-Martinelli  [14],  Baroczy  [19],  Thorn  [20],  Zivi  [21], 
Turner-Wallis  [22],  and  Hughmark  [23]  are  among  the  correlations  available 
in  the  literature.   Reference  will  be  made  to  these  correlations  in  a 
later  chapter. 

INTERNALLY  FINNED  TUBES 

In  recent  years,  manufacturing  techniques  have  been  well  developed  to 
produce  a  wide  variety  of  internally  finned  tubes.   Numerous  experimental 
and  analytical  studies  on  the  effects  of  these  tubes  on  heat  transfer  and 
pressure  drop  in  single-phase  as  well  as  two-phase  flow  were  conducted  by 
different  investigators. 

Single-Phase-Studies. 

Watkinson  et  al.  [24,25,26]  experimentally  investigated  the  heat 
transfer  and  pressure  drop  with  internally  finned  tubes  of  different 
diameters  and  fin  geometries  in  single-phase  flow.   An  enhancement  of  heat 
transfer  as  high  as  170%  at  the  same  Reynolds  number  and  up  to  80%  at  con- 
stant pumping  power  was  reported  [24]  for  turbulent  flow  of  water  in 
internally  finned  tubes.   A  95%  enhancement  of  heat  transfer  for  air  [25] 
and  224%  for  laminar  flow  of  oil  [26]  was  obtained  using  internally  finned 
tubes.   In  the  above  papers,  it  was  also  reported  that  spiralling  the  fins 
improved  the  heat  transfer  with  only  moderate  pressure  drop  increase. 

Carnavos  [27,28,29]  conducted  experiments  to  cool  air,  and  to  heat 
water  using  internally  finned  tubes.   In  his  study  of  cooling  air  [27],  he 
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used  21  tubes  of  integral  spiral  and  longitudinal  fins.   The  finned  tubes 

were  reported  to  perform  better  than  smooth  tubes  by  factors  of  1.2-2.0 
at  constant  pumping  power.   He  also  presented  correlating  equations  to  pre- 
predict  the  heat  transfer  and  friction  factor  for  internally  finned  tubes. 
In  [28],  Carnavos  experimentally  investigated  the  heat  transfer  performance 
of  five  composite  tubes,  made  by  mechanically  coupling  in  parallel  indi- 
vidual copper  tubes  having  continuous  integral  internal  spiral  and  longi- 
tudinal fins,  for  cooling  air  in  turbulent  flow.   He  reported  that  at 
constant  pumping  power,  the  composite  tubes  performed  better  than  smooth 
tubes  by  factors  of  1.7  to  10.   Correlations  to  predict  heat  transfer  and 
pressure  drop  were  also  presented.   In  [29],  Carnavos  experimentally  investi- 
gated the  heat  transfer  performance  of  heating  water  and/or  a  50%-50% 
ethylene  glycol-water  solution  in  turbulent  flow  by  tubes  having  internal 
spiral  and  longitudinal  fins  and  found  that  these  tubes  performed  better 
than  smooth  tubes. 

Compared  to  the  experimental  investigations,  the  analytical  studies  of 
heat  transfer  and  pressure  drop  were  limited. 

Hu  and  Chang  [30]  studied  the  heat  transfer  of  fully  developed  laminar 
flow  in  internally  finned  tubes  analytically.   The  fins  considered  had 
fictitious  zero  thickness.   They  concluded  that  the  presence  of  internal 
fins  improves  the  heat  transfer  performance  for  laminar  flow  more  than 
turbulent  flow  and  the  Nusselt  number  of  laminar  flow  in  tubes  with  optimum 
number  of  internal  fins  of  particular  height  can  surpass  that  of  many  cases 
of  turbulent  flow  in  f inless  tubes . 

Nandakumar  and  Masliyah  [31]  used  the  finite  element  method  of  analysis 
to  solve  the  momentum  equation  describing  the  laminar  fluid  flow  in  a 
finned  tube.   They  developed  empirical  correlations,  covering  a  wide 
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range  of  fin  parameters,  for  the  friction  factor  f  and  the  Nusselt 

number.   In  a  later  paper  [32],  they  obtained  the  heat  transfer  coef- 
ficients for  forced  convection,  fully  developed  laminar  flow  in  an 
internally  finned  tube,  with  axial  uniform  heat  flux  and  with  periph- 
erally uniform  temperature,  using  the  finite  element  method.   They  con- 
cluded that  the  Nusselt  number  of  a  triangular  finned  tube  was  higher 
than  that  of  a  finless  tube,  and  was  a  strong  function  of  fin  length 
and  fin  thickness. 

Soliman  and  Feingold  [33]  conducted  an  analysis  of  fully  developed 
laminar  flow  in  internally  finned  tubes  by  solving  the  momentum  equation 
using  infinite  series  coefficients  matching  method  to  obtain  the  velocity 
profile  and  friction  factor.   They  reported  that  the  values  of  the  dimen- 
sionless  velocity  at  any  location  within  the  tube  was  found  to  depend 
mainly  on  the  number  of  fins  and  their  height,  and  to  a  lesser  degree  on 
the  fin  half-angle.   In  a  later  paper  [34],  they  analytically  solved  the 
energy  equation  of  fully  developed  laminar  flow  to  obtain  the  temperature 
profile  and  heat  transfer,  and  concluded  that  the  Nusselt  number  increases 
with  the  increase  of  the  number  of  fins  up  to  a  critical  fin  number  bevond 
which  a  reversal  of  trend  occurs. 

Two-Phase  Flow  Studies. 

A  limited  number  of  studies  on  augmentation  of  in-tube  condensation 
has  been  reported  in  the  literature.   Owing  to  the  complexity  of  the  conden- 
sation mechanism  in  internally  finned  tubes,  no  analytical  studies  have 
been  attempted.   The  correlations  developed  to  predict  the  heat  transfer 
and  pressure  drop  during  condensation  inside  internally  finned  tubes  are 
based  on  smooth  tube  correlations  and  are  semi-empirical  in  nature. 

Reisbig  [35]  studied  the  condensation  of  R-12  inside  horizontal  tubes 
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with  internal  longitudinal  fins.   He  reported  a  40%  increase  in  heat  trans- 
fer compared  to  smooth  tubes.   He  also  found  that  the  finned  tubes  per- 
formed best  in  two-phase  flow  and  that  these  tubes  were  not  advantageous 
at  extreme  qualities. 

Kroger  [36]  studied  the  laminar  condensation  heat  transfer  of  R-12 
inside  smooth  and  internally  finned  tubes.   His  experimental  study  included 
examination  of  the  effects  of  tube  inclination,  and  different  tube  diameters 
to  length  ratios.   He  noted  that  in  a  horizontal  position,  the  rate  of  heat 
transfer  for  a  large  diameter  tube  having  twelve  fins  was  almost  the  same 
as  that  for  a  smaller  tube  having  nine  fins.   The  heat  transfer  rate  at  an 
inclination  of  10  degrees  was  found  to  be  up  to  200%  greater  than  that  of 
an  equivalent  smooth  horizontal  tube. 

Vrable  et  al.  [37,38,39]  investigated  the  condensation  heat  transfer 
and  pressure  drop  during  condensation  of  R-12  inside  two  longitudinally 
finned  tubes.   They  reported  that  the  addition  of  extended  surfaces  on  the 
inner  surface  of  the  tubes  resulted  in  approximately  20%  increase  in  the 
average  tube  side  heat  transfer  coefficient.   They  also  proposed  a  heat 
transfer  correlation  based  on  Cavallini  and  Zecchin  [12,40]  correlation 
for  smooth  tubes. 

Royal  and  Bergles  [41-43]  conducted  an  experimental  investigation  on 
in-tube  condensation  of  steam  with  three  augmentation  techniques:   tube 
inclination,  twisted  tape  inserts,  and  internally  finned  tubes.   Their  con- 
clusions were:   (a)  the  effects  of  tube  inclination  on  in-tube  heat 
transfer  coefficients  were  small.   An  inclination  of  10  degrees  downwards 
in  the  direction  of  flow  increased  the  heat  transfer  coefficient  by  a  maxi- 
mum of  10%  over  the  horizontal  value,  on  a  nominal  area  basis.   (b)  twisted 
tape  inserts  enhanced  the  heat  transfer  coefficients  by  as  much  as  30% 
above  smooth  tubes,  on  a  nominal  area  basis.   (c)  internally  finned  tubes 
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enhanced  the  heat  transfer  coefficients  by  as  much  as  150%  above  smooth 
tube  values,  on  a  nominal  area  basis.   The  heat  transfer  and  pressure 
drop  correlations  developed  by  Royal  and  Bergles  for  internally  finned 
tubes  are  discussed  in  Chapter  V. 

Luu  and  Bergles  [5,44-46]  condensed  R-113  inside  rough  tubes,  tubes 
with  twisted  tape  inserts,  and  internally  finned  tubes.   They  reported 
that  twisted  tape  inserts  increased  the  heat  transfer  coefficients,  on  a 
nominal  area  basis,  by  approximately  30%  and  pressure  drop  by  250%,  over  by 
smooth  tube  results.   Internally  finned  tubes  increased  the  heat  transfer 
as  much  as  120%  over  that  of  smooth  tube  with  only  a  modest  pressure  drop 
increase.   They  concluded  that  the  improvement  in  heat  transfer  performance 
of  internally  finned  tubes  were  far  better  than  that  of  twisted  tape  inserts. 
They  developed  correlations  to  predict  the  heat  transfer  and  pressure  drop 
during  condensation  inside  internally  finned  tubes  and  tubes  with  twisted 
tape  inserts.   The  correlations  of  internally  finned  tubes  are  discussed  in 
Chapter  V. 

Said  and  Azer  [2,47]  investigated  the  effects  of  twisted  tape  inserts 
and  internally  finned  tubes  on  augmenting  the  condensation  heat  transfer  of 
R-113.   The  twisted-tapes  selected  were  of  different  pitch  to  diameter 
ratios.   The  four  internally  finned  tubes  tested  had  different  outside 
diameters  and  fin  geometries.   Their  results  showed  an  improvement  in 
heat  transfer  coefficients  by  internally  finned  tubes  as  high  as  51%  over 
smooth  tube  results,  on  a  nominal  area  basis.   The  heat  transfer  was  ac- 
companied by  a  modest  increase  in  the  pressure  drop.   The  twisted  tape 
inserts  enhanced  the  heat  transfer  coefficients  by  approximately  23%  over  the 
smooth  tube,  on  a  nominal  basis,  and  was  accompanied  by  a  significant 
increase  in  the  pressure  drop.   The  correlation  equations  proposed  for 
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heat  transfer  and  pressure  drop  during  condensation  by  Said  and  Azer  are 
discussed  in  Chapter  V. 
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Chapter  III 
EXPERIMENTAL  INVESTIGATION 

The  primary  objectives  of  this  experimental  investigation  were: 

1.  to  obtain  experimental  heat  transfer  and  pressure  drop  data  for  con- 
densation of  R-ll  inside  horizontal  tubes  with  and  without  internal 
fins  . 

2.  to  correlate  the  experimental  heat  transfer  and  pressure  drop  results 
with  the  existing  correlations  and  to  identify  among  these  correlations 
the  ones  that  correlate  best  with  the  data. 

3.1  TEST  FACILITY 

The  apparatus  used  in  the  study  was  originally  constructed  bv  Said  [2] 
for  studying  the  enhancement  of  in-tube  condensation  by  internally  finned 
tubes  and  tubes  with  twisted  tape  inserts.   He  used  R-113  as  the  condensing 
fluid.   A  minor  modification  was  made  on  the  test  facility  to  allow  for 
using  R-ll  as  the  condensing  fluid.   This  modification  required  the  instal- 
lation of  a  liquid  R-ll  cooler  between  the  liquid  receiver  and  the  circu- 
lating pump.   Figure  3.1  shows  a  photographic  view  of  the  test  facility. 

3.2  R-ll  FLOW  LOOP 

Figure  3.2  shows  a  schematic  diagram  of  R-ll  flow  loop.   It  includes 
the  following  components: 

1.  an  electrically  heated  boiler 

2.  a  superheater 

3.  a  liquid  circulating  gear  pump 

4.  a  liquid-vapor  separator 
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5.  an  after  condenser 

6.  a  liquid  receiver 

7.  liquid  flowmeters 

8.  a  liquid  drier 

In  addition  to  the  above  components,  the  test  facility  included  five 
test  condensers  which  were  mounted  horizontally  in  parallel  as  shown  in  Fig. 
3.2.   Each  test  condenser  consisted  of  four  subsections.   Each  subsection 
was  a  double  pipe  counter  flow  heat  exchanger  where  R-ll  condensed  inside 
the  inner  tube  and  the  coolant  flowed  in  the  annulus .   Among  the  five  inner 
tubes,  one  was  smooth,  the  second  had  straight  internal  fins,  the  third, 
fourth  and  fifth  had  spiral  internal  fins  with  different  helix  angles  and 
outside  diameters. 

3.2.1   Test  Condenser  Construction  and  Instrumentation 

The  geometric  parameters  of  the  smooth  tube  and  the  internally  finned 
tubes  tested  are  given  in  Table  3.1.  Figure  3.3  shows  a  photograph  of  the 
finned  tubes  tested. 

Each  test  condenser  was  constructed  out  of  one  single  inner  tube  where 
R-ll  condensed.   Each  -condenser  consisted  of  four  subsections .   These  sub- 
sections were  constructed  by  soldering  four  equal  lengths  of  copper  tubing 
around  the  outside  surface  of  the  test  condenser  using  proper  fittings  and 
connections.   Each  subsection  represented  a  double-pipe  counter-flow  heat 
exchanger.   The  width  of  the  annular  space  in  each  of  the  subsections  was 
2.54  cms  (1  in.)  in  all  the  test  condenser  tubes.   The  outer  tube  of  each 
of  the  four  subsections  was  provided  with  coolant  inlet  and  outlet  con- 
nections and  connections  to  measure  the  outside  surface  temperature  of  the 
inner  tube.   The  connection  between  two  subsections  of  test  condenser  was 
done  by  plastic  tubing  as  shown  in  Fig.  3.4.   An  observation  section  was 


TABLE  3-1.  Geometric  Parameters  of  the  Experimental  Tubes- 
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Tuoe  No. 


Type 


Smooch 


Straight  Spiral    Spiral    Spiral 
Finned    Finned    Finned    Finned 


Material 

No.  of 
fins,  n 

Outside  dia- 
meter, D 

o 

Inside  dia- 
meter, D. 

Equivalent 

diameter,  D 
e 

Hydraulic 

diameter,  D^ 
B 

Fin  height,  b 

Wall  thick- 
ness 

Fin  height/ 
inside  dia- 
meter 


IG 


16 


16 


1.3843    0.853     0.858     1.13 


16 


1.5875    1.5850    1.9012    2.2162    2.6568 


1.3843    1.4199    1.7145    2.0384    2.5375 


1.3843    1.36      1.666     1.987     2.457 


1.50 


0.1575  0.1S03  0.1981         0.2134 

0.1016  0.0826  0.0934  0.0889  0.0597 

0.1109  0.1052  0.0972  0.0841 


Actual   flow  1.5050  1.4527  2.1799  3.1009  4.7413 


Nominal    flow        1.5050 

area,    A.   ** 

fn 


Core  flow- 
area,  A.  ** 

re 

Actual  area, 
A  *** 


Nominal  area,   4.34 
A  *** 


Inter- fan 
spacing,  W 


Helix  angle, 


1.5835    2.3087    3.2634    5.0571 


0.9588    1.4397    2.1181    3.4987 


4.3489    6.68 


10.2      11.3 


13.1 


.4607    5.3863    6.4038    7.9718 


0.297  0.206  0.305  0.363 


9.71°  12.34°  19.36° 


TABLE  3-1.    Geometric  Parameters   of    the   Experimental   Tubes 
(continued) 


Tube  No. 


Pitch, 
cm/ 360° 


ST       22.2      20.3      15.2 


*  all  lengths  and  areas  are  in  cm  and  cm2  respectively,  as  appropriate. 

2 
**   area  in  cm 

2 
*';*area  in  cm  /cm 
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installed  at  the  exit  of  each  condenser,  made  of  high  pressure  clear 

glass,  to  observe  the  condensate  flow.   Each  test  condenser  had  a 
flexible  copper  vibration  eliminator  section  attached  to  it  to  accom- 
modate any  expansion  due  to  heating  during  operation. 

Each  test  condenser  was  instrumented  to  measure  the  inlet  and  outlet 
temperatures  of  the  cooling  water  in  each  subsection,  inlet  pressure  and 
temperature,  and  the  exit  temperature  of  R-ll.   All  temperatures  were 
measured  by  copper-cons tantan  B  and  S  2k   gage  thermocouples.   All  thermo- 
couples were  connected  to  a  data  acquisition  system  which  gave  a  printout 
of  all  temperature  readings.   All  thermocouples  were  calibrated  at  two 
points,  the  boiling  water  and  melting  ice  temperatures.   The  inlet  pressure 
was  measured  with  a  Heise  pressure  gage  type  H28832. 

There  were  five  pressure  taps  to  measure  the  pressure  drop  across 
each  subsection.   The  location  of  these  pressure  taps  are  shown  in  Fig. 
3.5.   All  the  pressure  taps  had  four  holes,  each  0.16  cms  (1/16  in.) 
diameter,  spaced  90°  apart  around  the  circumference  of  the  test  condenser. 
These  holes  were  covered  with  a  copper  sleeve  which  had  a  clearance  of 
about  0.48  cms  (3/16  in.).   The  sleeve  was  silver  brazed  to  the  condenser 
tube.   A  0.32  cms  (1/8  in.)  I.D.  copper  tube  was  attached  to  the  sleeve 
and  connected  to  the  pressure  measurement  device. 

The  pressure  drop  across  any  two  pressure  taps  was  measured  by  a 
properly  calibrated  Foxboro  Differential  Pressure  (D/P)  cell  type  13A. 
The  D/P  cell  pressure  drop  measurements  were  checked  against  the  measure- 
ments of  a  Pace  Wianco  pressure  transducer  (P/T)  model  KP15.   The  P/T  was 
connected  to  a  Pace  Wianco  digital  indicator,  model  CD25,  which  gave  the 
pressure  difference  across  any  two  taps  directly.   Figure  3.6  shows  a 
schematic  diagram  of  the  piping  circuit  for  pressure  measurements.   In 
general  the  D/P  cell  and  P/T  measurements  agreed.   Even  when  they  differed, 
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the  difference  was  small  and  then  the  average  of  two  readings  were 

taken. 

The  average  wall  temperature  of  each  test  condenser  tube  of  each  sub- 
section was  measured  by  four  thermocouples  attached  as  shown  in  Fig.  3.7. 
Two  thermocouples  were  attached  at  the  top,  the  other  two  at  the  bottom. 
An  axial  distance  of  0.61m  (24  in.)  was  maintained  between  the  two  thermo- 
couple stations  located  at  0.152m  (6  in.)  from  the  coolant's  inlet  and 
outlet.   The  construction  details  of  each  thermocouple  junction  is  shown 
in  Fig.  3.8.   Conical  grooves,  slightly  larger  than  the  bead  of  thermo- 
couple, were  drilled  and  the  beads  were  silver  brazed  in  these  grooves. 
Fast  drying  epoxy  (a  product  of  Armstrong  company  type  A-36/B-36)  was 
applied  to  the  silver  brazed  portion  of  thermocouple  and  a  short  length 
of  nylon  heat  shrink  tubing  was  slid  from  the  other  end  towards  the  junction 
until  it  touched  the  outside  surface  of  the  condenser  test  tube  and  covered 
the  junction.   More  epoxy  was  added  to  secure  the  nylon  tube  to  the  outside 
surface.   This  arrangement  ensured  the  elimination  of  contact  between  the 
thermocouple  junction  and  cooling  water.   The  outer  surface  of  the  condenser 
was  then  covered  with  fiber  glass  insulating  tape  near  the  thermocouple 
junction.   The  thermocouple  leads  were  then  tied  to  the  test  condenser 
tube  using  nylon  cable  tie.   Additional  epoxy  was  used  to  cover  the  thermo- 
couple leads  between  the  junction  and  nylon  tie.   The  fiber  glass  tape 
protected  the  thermocouple  leads  from  getting  heated  at  the  portion  where 
they  were  tied  to  condenser  tube  and  the  nylon  tie  prevented  the  detach- 
ment of  the  thermocouple  junction  by  accidental  pulling  of  the  thermocouple 
leads  from  the  other  end.   The  thermocouple  leads  were  then  led  through  the 
outside  tube  as  shown  in  Fig.  3.8.   The  space  between  the  leads  and  the 
copper  tube  attached  to  the  outer  tube  of  each  subsection  was  filled  with 
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Fig.  3.8  Thermocouple  Mounting  Detail. 
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epoxy  to  prevent  any  coolant  leaks.   The  thermocouple  junctions  were 

checked  for  good  thermal  contact  after  their  completion. 

The  flow  of  R-ll  leaving  the  circulating  pump  was  measured  by  two 
Fisher  and  Porter  variable  area  flowmeters.   The  flowmeters  were  mounted 
vertically  in  parallel  so  that  they  could  be  used  individually  or  simul- 
taneously to  measure  the  flow  rate.   One  had  a  range  of  0.0-0.35  g.p.m. 
and  the  other  0.0-0.55  g.p.m.   There  was  a  similar  flowmeter  connected 
between  the  vapor-liquid  separator  and  the  liquid  receiver  to  measure  the 
flow  rate  of  the  uncondensed  vapor  inside  the  test  condenser.   In  the 
present  study  complete  condensation  was  achieved  inside  the  test  condenser 
in  all  runs  and  hence  this  was  never  used.   Electric  heating  elements  were 
used  to  heat  R-ll  in  the  vapor  generator.   The  power  was  controlled  by  a 
powerstat  variable  transformer.   The  power  input  to  each  of  the  heating 
elements  was  measured  by  measuring  the  current  and  voltage  across  each 
heating  circuit.   Appendix  A  lists  additional  details  about  the  components 
of  the  refrigerant  flow  circuit. 

Individual  test  condensers  and  the  entire  refrigerant  circuit  was 
leak-proof  tested  under  pressurized  and  evacuated  conditions.   A  layer  of 
rubber  insulation,  0.064m  (2.5  in.)  thick,  prevented  heat  loss  from  the  test 
condenser  tube  to  the  ambient.   The  vapor  generator  was  also  insulated 
with  0.019m  (3/4  in.)  thick  layer  of  fiber  glass  insulation. 

3.3  Cooling  Water  Loop 

City  supply  water  was  used  as  a  coolant.   Figure  3.9  shows  a  schematic 
diagram  of  the  cooling  water  flow  loop.   The  city  water  was  supplied  to  a 
mixing  tank  from  the  water  main.   The  water  was  pumped  through  a  circu- 
lating pump  to  the  main  test  condenser,  after  condenser,  and  the  R-ll 
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liquid  cooler.   The  water  flow  rate  to  the  test  condenser  was  measured  by 

Brooks  Rotameter  (range  0.0  to  3.0  g.p.m.).   The  design  of  the  cooling 
water  loop  permitted  control  of  the  flow  through  test  condenser,  the  R-ll 
liquid  cooler  and  the  cooling  water  inlet  temperature.   The  R-ll  liquid 
cooler,  installed  between  the  liquid  receiver  and  the  refrigerant  circu- 
lating pump,  ensured  the  liquid  phase  of  R-ll  at  the  inlet  of  the  refrige- 
rant pump.   The  coolant's  inlet  temperature  was  achieved  by  mixing  appro- 
priate quantities  of  supply  water  from  the  mains  and  return  water  from  the 
test  condenser,  in  the  mixing  tank. 

3.4  Operation  of  Test  Facility  and  Data  Acquisition 

The  entire  system  was  evacuated  down  to  1000  microns  by  a  vacuum  pump 
and  then  charged  with  R-ll.   The  process  of  evacuation  and  charging  was 
repeated  two  to  three  times  to  rid  the  system  out  of  non-condensables. 

The  circulating  water  pump  was  put  into  operation  and  the  flow  rate 
was  adjusted  to  the  desired  level.   Before  applying  the  electric  power  to 
the  heating  elements  of  the  vapor  generator,  it  was  ensured  that  the 
heating  elements  were  well  immersed  in  the  liquid  by  checking  liquid  level 
through  the  sight  glass  attached  to  it.   If  the  liquid  level  was  low,  R-ll 
was  pumped  from  the  liquid  receiver  into  the  vapor  generator  until  the 
desired  level  was  achieved.   Electric  power  was  then  gradually  applied 
to  the  heating  elements  inside  the  vapor  generator.   When  the  liquid 
level  started  decreasing  in  the  vapor  generator,  the  refrigerant  circu- 
lating pump  was  put  into  operation.   The  balance  of  flow  into  and  out  of 
the  vapor  generator  was  achieved  by  adjusting  the  bypass  valve  after  the 
refrigerant  circulating  pump.   Once  the  flow  balance  was  achieved,  elec- 
tric power  was  then  applied  to  superheater.   The  desired  flow  rate  of  R-ll 
into  the  test  condenser,  its  inlet  pressure  and  superheat  for  any  experi- 
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mental  run,  was  achieved  by  adjusting  the  electric  power  to  the  vapor 

generator  and  the  superheater,  controlling  the  liquid  flow  rate  from 

the  refrigerant  pump,  and  by  controlling  cooling  water  flow  rate  and 

its  inlet  temperature.   For  each  experimental  run,  steady  state  condition 

was  established  when  the  readings  of  all  temperatures,  pressures,  flow 

rates,  and  liquid  level  in  the  vapor  generator  and  liquid  receiver  remained 

constant  for  a  period  of  least  two  hours. 

During  each  experimental  run,  the  following  measurements  were  taken. 

1.  Water  flow  rate  to  the  test  condenser  in  kg/s  (g.p.m.) 

2.  R-ll  inlet  gage  pressure  to  the  test  condenser  kPa  (psia) 

3.  The  inlet  and  outlet  temperatures  of  R-ll  of  the  test  condenser,  in 
deg  C  (deg  F) 

4.  The  temperature  of  liquid  R-ll  leaving  the  flowmeter,  in  deg  C 
(deg  F) 

5.  The  temperature  of  sixteen  thermocouples  attached  to  the  outer  wall 
of  the  test  condenser,  deg  C  (deg  F) 

6.  The  flow  rate  of  refrigerant  into  the  test  condenser  in  kg/s  (g.p.m.) 

7.  The  inlet  and  outlet  water  temperatures  of  each  subsection 

8.  The  pressure  drop  across  the  test  condenser  (for  tubes  1,  2  and  3) 
in  mm  Hg  (in.  of  Hg) 

The  pressure  drop  across  the  test  condenser  was  not  taken  for  tubes 
4  and  5,  as  these  tubes  had  larger  diameters  compared  to  the  other  tubes 
tested  and  consequently,  their  pressure  drops  were  too  low  to  be  measured 
accurately. 

The  flow  rates  of  R-ll,  measured  by  the  variable  area  flowmeters,  were 
corrected  according  to  the  manufacturer's  recommendation  to  account  for 
density  variations  of  the  fluid  tested,  since  these  flow  meters  were 
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originally  calibrated  to  handle  R-113. 

Heat  balance  error  of  ±5%  was  set  as  the  criterion  for  acceptability 
of  any  experimental  run  before  documentation.   The  heat  balance  was  calcu- 
lated as  the  ratio  of  the  difference  between  heat  gain  rate  of  the  cooling 
water  and  the  rate  of  heat  loss  of  the  condensing  fluid  to  the  rate  of 
heat  gain  of  the  cooling  water.  The  heat  gain  rate  of  the  cooling 
water  was  calculated  from  its  flow  rate  and  its  temperature  rise.   The 
rate  of  heat  loss  of  the  condensing  fluid  was  calculated  from  its  flow 
rate  and  its  enthalpy  change  between  its  inlet  and  outlet.-  In  all  the 
experimental  runs,  R-ll  was  slightly  superheated  at  the  inlet  of  the  test 
condenser  and  slightly  subcooled  at  the  exit  of  test  condenser. 

The  heat  balance  error  for  the  experimental  runs  rarely  exceeded  ±5%. 
However,  a  few  runs  with  heat  balance  error  slightly  higher  than  the  above 
limits  were  also  included  in  the  data  reported  in  this  study. 

The  ranges  of  experimental  parameters  covered  in  this  study  were  dic- 
tated by  the  various  limitations  of  the  existing  experimental  set  up.  The 
ranges  of  the  experimental  parameter  are  listed  in  Table  3.2. 


TABLE  3.2  Ranges  of  Experimental  Parameters  Covered 
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Refrigerant's  Mass  Flux 
(based  on  nominal  inside  area) 


Overall  Condensing  Heat 
Transfer  Coefficient 


Overall  Heat  Transfer 
Rate 

Inlet  Coolant  Temperature 


Test  Fluid  Inlet  Pressure 


Test  Fluid  Inlet  Superheat 


Outlet  Quality 


17.14(1.264x10*)  -  85.55(6. 308xl04) 
Kg   (lbm/hr-ft2) 
s-m2 

1044.2(183.9)  -  3456.4(608.7) 

W   (Btu/hr-ft2-°F) 
m2-°C 

1268(4331)  -  2924(9983) 
H  (Btu/hr) 

14.7(58.64)  -  25.0(75.56) 
°C  (°F) 

127.2(18.45)  -  172.7(25.04) 
kPa(psia) 

2.46(4.42)  -  9.86(17.74) 
°C  (°F) 
-.158  (Subcooled)  -  0.0 
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Chapter  IV 

EXPERIMENTAL  RESULTS 

4.1  Introduction 

During  the  present  investigation,  fifty-two  experimental  runs  were 
taken  for  the  smooth  tube  (tube  1),  fifty-two  runs  for  the  straight  finned 
tube  (tube  2),  fifty-nine,  thirty-three  and  thirty-six  runs  for  the  spirally 
finned  tubes,  tubes  3,  4  and  5,  respectively.   For  the  reasons  discussed 
earlier,  pressure  drop  data  were  taken  only  for  the  smooth  tube  (tube  1), 
the  straight  finned  tube  (tube  2)  and  the  spirally  finned  tube  (tube  3). 
Nineteen  pressure  drop  readings  were  taken  for  the  smooth  tube,  nineteen 
for  the  straight  finned  tube  and  eighteen  for  the  spirally  finned  tube. 
In  these  pressure  drop  measurements,  only  the  total  drop  across  the  entire 
test  condenser  tube  was  recorded. 

The  sectional  and  overall  average  heat  transfer  coefficients  were 
calculated  from  the  following  equation. 

fD 

c 
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"it  D.  Az  (T   -  T   )        D. 

1        s     wo       -   1  . 

2K  £n 


m  Cp    (T    -  T   .) 
wa  wa   wao    wai 


(4-1) 


In  the  above  equation,  Eq.  (4-1),  T  is  the  saturation  temperature 
corresponding  to  the  inlet  pressure  of  the  test  condenser  tube  since  the 
pressure  drop  in  all  experimental  runs  was  small  enough  to  have  no  signifi- 
cant effect  on  the  change  in  saturation  temperature. 

While  computing  the  sectional  heat  transfer  coefficients  for  section  1 
(refer  to  Fig.  3.4.),  for  each  test  condenser  using  Eq.  (4-1),  it  was  found 
that  the  resulting  heat  transfer  coefficient  was  very  high.   The  reason  for 
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such  results  was  due  to  the  fact  that  the  temperature  difference  (T  -T   ) 

s  wo 

was  very  small  to  yield  very  high  h.   It  was  then  reasoned  that  basing 

the  heat  transfer  coefficient  of  such  a  temperature  difference  was  not 

realistic  due  to  the  fact  that  R-ll  entered  this  section  superheated  and 

its  temperature  remained  above  the  saturation  temperature  T  over  a  major 

portion  of  this  section.   It  was  then  decided  to  base  h  for  this  section 

on  the  difference  [ (T  +T   )/2-T   ]  where  T    is  the  inlet  superheat 
s   sup     wo         sup 

temperature  of  R-ll. 

In  calculating  the  overall  heat  transfer  coefficient  of  the  entire 

test  condenser  from  Eq.  (4-1)  ,  iz  was  replaced  by  L,  the  combined  lengths 

of  four  sections.   T   was  taken  as  the  arithmetic  mean  of  the  sixteen 
wo 

thermocouples  attached  to  the  outside  wall  temperature  of  the  test  con- 
denser tube.   The  total  energy  transfer  was  taken  as  the  sum  of  energy 
gain  rate  of  the  cooling  water  for  the  entire  test  condenser. 

A  sample  of  the  calculation  procedure  is  given  in  Appendix  B. 
Appendix  C  gives  the  computer  program  for  the  data  reduction,  and  Appendix 
D  gives  the  reduced  data  of  all  the  experimental  runs.   The  estimation  of 
experimental  uncertainties  of  the  computed  overall  average  heat  transfer 
coefficient  is  given  in  Appendix  E.   The  uncertainty  was  estimated  to  be 
±16.15%. 

4.2  Heat  Transfer  Results 

From  the  reduced  data  of  the  experimental  runs,  the  following  general 
observations  can  be  made  on  the  heat  transfer  results. 

1.  The  average  sectional  heat  transfer  coefficients  of  each  subsection 
nearest  to  the  inlet  was  highest  among  all  subsections  and  decreased 
in  the  direction  of  flow  for  the  remaining  three  subsections. 

2.  Heat  gained  by  water  in  the  four  subsections  were  close  enough  in 
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majority  of  experimental  runs  which  indicated  that  the  quality  of 

condensing  vapor  decreased  linearly  in  the  direction  of  flow. 

3.   Also,  by  carrying  out  the  energy  balance  on  each  subsection,  starting 

from  the  one  nearest  to  the  inlet  and  then  proceeding  to  the  next 

subsection  in  the  direction  of  flow,  it  was  possible  to  estimate 

the  exit  quality  of  all  the  subsections. 

4.2.1  Smooth  Tube  and  Internally  Finned  Tubes 

The  experimental  runs  were  carried  out  at  three  nominal  inlet  pres- 
sures, 1.32  bar  (19.19  psia. ) ,  1.47  bar  (21.29  psia.),  and  1.67  bar  (24.17 

2 
psia.)(l  bar  =  10  kPa) .  The  mass  flux  based  on  the  inside  nominal 

area  ranged  from  41.33  (3.05xl04)  to  85.55  (6.31xl04)  kg/s.m2  (lbm/hr.ft2) 

for  tube  1,  54.00  (3.98xl04)  to  84.57  (6.24xl04)  kg/s.m2  (lbm/hr.ft2)  for 

tube  2,  37.20  (2.74xl04)  to  66.46  (4.90xl04)  kg/s.m2  (lbm/hr.ft2)  for  tube 

3,  26.35  (1.94xl04)  to  47.01  (3.47xl04)  kg/s.m2  (lbm/hr.ft2)  for  tube  4, 

and  17.14  (1.26xl04)  to  28.45  (2.10xl04)  kg/s.m2  (lbm/hr.ft2)  for  tube  5. 

Constant  inlet  pressures  of  the  condensing  fluid  could  not  be  main- 
tained precisely  because  of  the  practical  difficulties  in  trying  to  main- 
tain reasonable  inlet  superheat  and  exit  subcooling  for  the  desired  ranges 
of  mass  fluxes.  The  variation  of  the  inlet  pressures  of  all  the  tubes 
reported  for  different  mass  fluxes  of  the  condensing  fluid  were  within 
±0.028  bar  (+0.40  psia.)  from  the  nominal  pressures.   The  mass  flux  ranges 
for  different  tubes  were  dictated  by  the  practical  limitations  of  the  test 
facility. 

The  overall  average  heat  transfer  coefficients  at  three  nominal  inlet 
pressures  for  the  smooth  tube  are  plotted  versus  mass  flux  in  Figs.  4.1 
through  4.3.   Figures  4.5  through  4.7  are  the  plots  of  the  overall  heat 
transfer  coefficient  versus  the  mass  flux  at  different  nominal  inlet 
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4. 1  Experimental  overall  average  heat  transfer  coefficients  versus 
mass  flux,  Tube  1,  Pin=1.32  bar  (19.19  Psia) 
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Fig  4.2  Experimental  overall  average  heat  transfer  coefficients  versus 
mass  flux,  Tube  1,  Pin-1.47  bar  (21.29  Psia) 
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pressures  for  tube  2,  Figs.  4.9  through  4.11  are  similar  plots  for  tube  3, 

Figs.  4.13  through  4.15  for  tube  4,  and  Figs.  4.17  through  4.19  for  tube  5. 

Because  of  the  trend  of  the  data,  it  was  possible  to  correlate  the 

results  of  all  the  tubes  at  each  nominal  pressure  by  the  following  equation: 

h  =  C  Gn  (4-2) 

where  C  and  n  are  constants  obtained  by  the  least  square  regression  analysis. 
Values  of  C  and  n  are  given  in  Appendix  F,  Table  F.l,  for  all  the  tubes  at 
all  nominal  inlet  pressures.   Figure  4.4  shows  a  plot  of  the  regression 
equations  of  h  versus  G  at  all  inlet  pressures  for  tube  1.   Similarly,  the 
regression  equations  of  h  versus  G  at  all  inlet  pressures  are  shown  in  Fig. 
4.8  for  tube  2,  Fig.  4.12  for  tube  3,  Fig.  4.16  for  tube  4,  and  Fig.  4.20 
for  tube  5.   The  results  of  all  the  tubes  show  that,  for  the  same  mass  flux, 
h  decreased  with  the  increase  in  pressure.   This  is  due  to  the  fact  that, 
for  a  given  G,  the  density  decreases  as  the  pressure  decreases  and  as  a 
result  the  velocity  of  the  condensing  fluid  increases  resulting  in  the 
increase  of  the  h  value. 

4.3  Pressure  Drop  Results 

As  reported  earlier,  pressure  drop  data  were  taken  for  tubes  1,  2,  and 
3  only.   Nineteen  runs  of  pressure  drop  for  tube  1,  nineteen  for  tube  2, 
and  eighteen  for  tube  3  were  taken  at  two  nominal  inlet  pressures,  1.32 
bar  (19.19  psia.)  and  1.67  bar  (24.17  psia.). 

The  total  pressure  drop,  Ap  at  the  two  nominal  inlet  pressures  for 
the  smooth  tube  are  plotted  versus  the  mass  flux  in  Figs.  4.21  and  4.22. 
Figures  4.24  and  4.25  are  plots  of  the  total  pressure  drop  versus  mass 
flux  at  the  two  nominal  inlet  pressures  for  tube  2,  and  Figs.  4.27  and 
4.28  for  tube  3.   The  results  were  correlated  by  the  equation 
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Fig  4.4  Experimental  overall  average  heat   transfer  coefficients  versus 
mass  flux,      Tube   1,   at  all  pressures 
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Fig  4.5     Experimental  overall   average  heat   transfer  coefficients 
versus  mass   flux,    tube  2,    P.    =1.32  bar    (19.19  psia) 


41 


03 

c 

I 
_ 

X 

I 


\ 

3 


2.Sxl02L 


G       Kq/m2- 


50 


53 


3. BX IB4 


-J 1 1 1 1 '      '     ' 


sec 

70  80  30       100 


I 


I      '      I 


3.5 


- 


a 

"0 


£ 
\ 


|J= 


'    I    '    '    !    I    I   i 


2 


1.5x10 


bm/ft2-hr 


Fig  4.6  Experimental  overall  average  heat  transfer  coefficients 
versus  mass  flux,  Tube  2,  P.  =1.47  bar  (21.29  psia) 
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Fig  4.7  Experimental  overall  average  heat  Transfer  coefficients 
versus  mass  flux,  tube  2,  P.  =1.67  bar  (24.17  psia) 
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Fig   4.8  Experimental  overall  average  heat   transfer  coefficients 
versus  mass   flux,   Tube   2,    at   all  pressures 
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Fig  4.9  Experimental  overall  average  heat  transfer  coefficients 
versus  mass  flux,  Tube  3,  Pin-=1.32  bar  (19.19  Psia) 
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Fig  4.10  Experimental  overall  average  heat  transfer  coefficients 
versus  mass  flux,  Tube  3,  Pin- 1.47  bar  (21.29  Psia) 
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Fig  4.11  Experimental  overall  average  heat  ransfer  coefficients 
versus  mass  flux,  Tube  3,  Pin- 1.6 7  bar  (24.17  Psia) 
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Fig  4.12  Experimental  overall  average  heat  transfer  coefficients 
versus  mass  flux,  lube  3,  at  all  pressures 
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Fig  4.13  Experimental  overall  average  heat  transfer  coefficients 
versus  mass  flux,  Tube  4,  Pin-1.32  bar  (19.19  Psia) 
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g  4.14  Experimetnal  overall  average  heat  transfer  coefficients 
versus  mass  flux,  Tube  4,  Ptn-1.47  bar  (21.29  Psia) 
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Fig  4.15  Experimental   overall   average  heat   transfer  coefficients 
versus  mass  flux   ,    Tube  4,   Pin=1.67   bar    (24.17  Psia) 
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Fig  4.16  Experimental  overall  average  heat  transfer  coefficients 
versus  mass  flux,  Tube  4,  at  all  pressures 
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average  heat  transfer  coefficients 
versus  mass   flux,    Tube   5,    Pin-    1.32   bar    (19.19  Psia) 
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Fig  4.18  Experimental  overall  average  heat  transfer  coefficients 
versus  nass  flux,  Tube  5,  Pin=1.47  bar  (21.29  Psia) 
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Fig  4.19   Experimental  overall   average   heat   transfer  coefficients 
versus  mass   flux,    Tube  5,    Pin=1.67   bar   (24.17  Psia) 
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Fig.    4.20   Experimental   overall  average  heat   transfer  coefficients 
versus  mass   flux,    tube   5,    at   all  pressures 
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Lg  4.21  Experimental  total  pressure  drop  versus  mass  flux, 
Tube  1,  Pin-1.32  bar  (19.19  Psia) 
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Fig  4.22  Experimental  total  pressure  drop  versus  mass  flux, 
Tube  1,  Pin- 1.67  bar  (24.17  Psia) 
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Fig  4.23  Experimental  total  pressure  drop  versus  mass  flux, 
tube  1,  at  all  pressures 
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Fig  4.24  Experimental  total  pressure  drop  versus  mass  flux, 
Tube  2,  Pin=1.32bar  (19.19  Psia) 
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Fig  4.25  Experimental  total  pressure  drop  versus  mass  flux, 
Tube  2,  Pin=1.67  bar  (24.17  Psia) 
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Fig  4.26  Experimental  total  pressure  drop  versus  mass  flux,. 
Tube  2,  at  all  pressures 
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Fig  4.27  Experiment  total  pressure  drop  versus  mass  flux, 
Tube  3,  Pin=1.32  bar  (19.19  Psia) 
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Fig  4.28  Experimental  total  pressure  drop  versus  mass  flux, 
Tube  3,  Pln-1.67  bar  (24.17  Psia) 
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Fig  4.29  Experimental  total  pressure  drop  versus  mass  flux, 
Tube  3,  at  all  pressures 
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AP  =  C  G  (4-3) 


The  constants  C  and  n'  were  obtained  by  the  least  square  regression  tech- 
nique.  The  values  of  C  and  n'  for  the  tubes  1,  2  and  3  tested  at  all 
pressures  are  given  in  Appendix  F,  Table  F.2.   Figure  4.23  shows  the  least 
square  regression  lines  for  the  smooth  tube  at  two  different  pressures. 
For  the  smooth  tube,  it  was  observed  that  higher  inlet  pressures  resulted 
in  higher  pressure  drops.   The  regression  lines  at  the  two  different  inlet 
pressures  are  shown  in  Fig.  4.26  for  tube  2,  and  Fig.  4.29  for  tube  3.   For 
the  straight  finned  tube  (tube  2),  the  total  pressure  drop  was  high  at  the 
lower  inlet  pressure  and  higher  mass  flux  compared  to  that  at  higher  inlet 
pressure.   For  the  spiral  finned  tube  (tube  3)  higher  inlet  pressure 
yielded  marginally  lower  total  pressure  drop.   In  all  the  above  plots, 
results  show  that  Ap  increased  as  G  increased  for  the  two  nominal  inlet 
pressures. 

The  performance  of  all  the  tubes  tested  in  terms  of  their  heat  trans- 
fer and  the  associated  pressure  drop,  will  be  evaluated  in  Chapter  VI. 
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Chapter  V 
CORRELATION  OF  EXPERIMENTAL  DATA 

5.1  Introduction 

There  are  numerous  correlations  available  for  predicing  the  heat  trans- 
fer and  pressure  drop  during  in-tube  condensation.   Before  relying  on  these 
correlations  for  commercial  design  usage,  these  correlations  should  predict 
the  data  of  different  investigators  over  a  wide  range  of  flow  conditions 
for  different  fluids  within  reasonable  degree  of  accuracy.   Different  heat 
transfer  and  pressure  drop  correlations  were  tested  to  identify  the  cor- 
relations which  properly  predicted  the  data  of  the  present  study. 

5.2  Smooth  Tube  Correlations 

5.2.1  Heat  Transfer  • 

The  smooth  tube  experimental  data  were  compared  with  the  predictions 
of  the  correlations  of  Akers  et  al.  [6],  Eoyko  and  Kruzhilin  [7],  and  Shah 
[13].   A  summary  of  these  correlations  follows. 

The  overall  average  heat  transfer  coefficient  of  Akers  et  al.  [6]  is 
given  by 

ki  n    1/3 
h  =  C  57  Ree  PrJ  (3-D 

1 

where 

C  -  0.0265,  n  =  0.8  for  Re  >  5  x  10* 
e 

C  =  5.03,  n  =  1/3  for  Re  <  5  x  104 
e 

GD," 
Ree  =  IT"  (5-2) 


0.5 


+  U 


x) 


(5-3) 
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The  above  correlation  is  recommended  by  ASHRAE  Handbook  of  Fundamentals 
[48]. 

The  overall  average  heat  transfer  coefficients  of  Boyko  and  Kruzhilin 
[7]  is  given  by 

k. 


"£  _  0.8  _  0.43 
h  =  0.024  —  Re    Pr 


(p/Pm),-„  +  (P/Pm)„,„. 
mm        m  out 


where 


1  +  x 


(5-4) 


(5-5) 


The  local  heat  transfer  coefficient  of  Shah  [13]  is  given  by 

-0.38 


h  =0.023 


<■"   *f 


(1-x)   +3.Sx    (1-x)     P 


(5-6) 

All  the  properties  in  the  above  correlations  are  evaluated  at  the 
saturation  temperature.   Shah  also  suggested  that  an  approximate  expres- 
sion for  the  average  h  for  complete  condensation  can  be  obtained  by  sub- 
stituting average  x  in  Eq.  (5-6).   After  analyzing  the  experimental  results 
of  several  investigators,  Shah  suggested  that  an  approximate  expression  for 
the  average  h  for  complete  condensation  can  be  obtained  by  substituting 
x=0.5  in  Eq.  (5-6).   He  also  pointed  out.  that  h  obtained  from  this  expres- 
sion differed  by  no  more  than  5%  from  averaging  the  values  obtained  by  the 
predictions  of  Eq.  (5-6).   After  evaluation  of  Shah's  correlation,  Said 
selected  x=0.4  and  when  it  was  substituted  in  Eq.  (5-6)  the  following 
expression  resulted. 


h=0.023 


D. 


Re°-8  Pr°-4  [0.665  +  1.86  P^'38] 


(5-7) 


The  average  heat  transfer  coefficient  of  the  predictions  of  Eq .  (5-6) 
can  be  evaluated  by  integration  using  the  following  equation. 

/ln  p  (5-8) 


r    (x.   -  X   „)   X     h 

n     in    out    out 


Equation  (5-8)  is  based  on  energy  balance  with  the  assumption  that  the 
difference  between  the  condensing  fluid  and  the  inner  wall  temperature  is 
constant  for  entire  length  of  the  tube.   Said  [2],  in  his  investigation  of 
in-tube  condensation  with  R-113,  correlated  his  experimental  data  of  h  with 
the  predicted  values  from  Eq.  (5-7)  and  the  prediction  of  Eq.  (5-6)  averaged 
by  Eq.  (5-8).   Eq.  (5-7)  predicted  his  data  well.   Hence  in  the  present 
study,  it  was  decided  to  use  Eq.  (5-7)  for  correlating  the  data  of  the 
present  study. 

Figures  5.1  and  5.2  show  the  comparison  between  the  predictions  of 
Akers  et  al.  [6],  Eq.  (5-1),  and  the  overall  experimental  values  and  the 
experimental  overall  and  sectional  values,  respectively,  of  the  heat  trans- 
fer coefficients.   Figure  5.1  shows  that  the  predictions  of  the  overall 
heat  transfer  coefficients  were  within  ±30%  of  experimental  results  for 
83%  of  the  experimental  runs.   Also,  the  correlation  predicted  reasonably 
well  the  sectional  average  heat  transfer  coefficient  values  shown  in 
Fig.  5.2.   The  sectional  average  values  of  the  heat  transfer  coefficients 
of  section  1  were  not  included  in  Fig.  5.2  for  two  reasons.   First,  Eq. 
(5-1)  substantially  underpredicted  the  values  obtained  experimentally. 
Secondly,  the  temperature  difference  on  which  basis  the  overall  average 
and  sectional  averages  of  sections  2,  3,  and  A  were  calculated, were  dif- 
ferent than  those  of  section  1  as  outlined  earlier. 

Figure  5 . 3  shows  the  comparison  between  the  experimental  overall 
heat  transfer  coefficients  and  predictions  of  Boyko  and  Kruzhilin's 
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Fig  5.1  Comparison  of  experimental  overall  heat  transfer 

coefficient  with  predictions  of  the  correlation  of 
Akers  et  al .  [6J,    Eq.(S-l),  Tube  1 
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Fig  5.3  Comparison  of  experimental  overall  heat  transfer 
coefficients  with  predictions  of  the  correlation 
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correlation  [7],  Eq.  (5-4).   The  correlation  underpredicted  the  experi- 
mental overall  heat  transfer  values. 

The  comparison  between  the  experimental  overall  heat  transfer  results 
and  predictions  of  Shah's  correlation  [13],  Eq.  (5-7),  is  shown  in  Fig. 
5.4.   The  correlation  also  underpredicted  the  measured  values. 

The  results  of  the  above  comparisons  are  summarized  below: 

1.  The  correlation  of  Alters  et  al.  [6],  Eq.  (5-1),  correlated  83%  of 
experimental  overall  heat  transfer  values  within  ±30%.   It  also  cor- 
related the  sectional  heat  transfer  values  reasonably  well,  except 
section  1. 

2.  The  correlation  of  Boyko  and  Kruzhilin  [7],  and  Shah  [13]  predicted 
overall  heat  transfer  values  lower  than  measured  values. 

It  is  worth  noting  that  among  the  recent  investigations  of  in-tube  con- 
densation by  Royal  and  Bergles  [4,42],  Luu  and  Bergles  [5,45,46],  and  Said 
and  Azer  [2,47],  the  smooth  tube  data  obtained  by  Royal  and  Bergles,  and 
Said  and  Azer  correlated  best  with  predictions  of  Akers  et  al.  [6].   Luu 
and  Bergles  data  correlated  well  with  the  predictions  of  Boyko  and  Kruzhilin 
[7]. 

5.2.2  Pressure  Drop 

The  pressure  gradient  during  two-phase  flow  inside  a  horizontal  tube 
is  usually  divided  into  two  parts,  friction  and  momentum  terms. 

(5-9) 

The  main  factors  involved  in  the  frictional  pressure  drop  are  probably 

the  viscous  dissipation  in  the  liquid  film  and  the  pressure  drop  associated 

with  gas  flow  over  the  wavy  interface.   The  momentum  recovery  part  of  the 
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pressure  drop  is  due  to  the  deposition  of  faster  moving  vapor  on  the 
slower  moving  condensate  flow. 

The  literature  survey  revealed  that  two  basic  models  are  usually  used 
in  analyzing  the  two-phase  pressure  drop.   In  the  first  model,  the  two 
phases  are  assumed  homogeneous  with  the  same  gas  and  liquid  velocities. 
The  frictional  pressure  drop  in  this  model  is  developed  by  introducing 
modifiers  to  the  single  phase  friction  coefficient.   In  the  second  model, 
the  two  phases  are  considered  to  have  different  velocities  and  hence  are 
separated. 
1.    The  frictional  pressure  gradient  correlations 

Among  the  well  recognized  frictional  pressure  gradient  correlations 
are  those  of  Lockhart-Martinelli  [14]  and  Dukler  [16].   A  summary  of 
these  correlations  follows. 

Lockhart-Martinelli  [14]  defined  a  parameter  X  given  by 
0.50 


/ 


(5-10) 


where 


_  frictional  pressure  gradient  assuming  that  liquid  alone  is 
flowing  in  the  pipe 


_  frictional  pressure  gradient  assuming  that  gas  (or  vapor  in 
case  of  condensation)  alone  is  flowing  in  the  pipe 


In  two-phase  flow,  four  combinations  of  flows  are  possible.   (a)  Liquid 
viscous,  gas  viscous,  (b)  Liquid  viscous,  gas  turbulent,  (c)  Liquid  tur- 
bulent, gas  viscous,  and  (d)  Liquid  turbulent,  gas  turbulent.   Liquid  tur- 
bulent, gas  turbulent  flow  is  the  one  usually  encountered  during  in-tube 
condensation.   For  turbulent  flow  of  the  two  phases  during  condensation, 


the  parameter  X  can  be  replaced  by  X   and  written  as 
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(5-11) 


Lockhart  and  Martinelli  also  defined  the  parameter  41  as 

0.5 


/ 


For  condensation 


(5-12) 


2  f  (x  G) 

o 

P  D. 
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where 


f  =  0.045/ 


0.2 


(5-13) 


(5-14) 


f   is  the  friction  factor  for  adiabatic  two-phase  flow.   Lockhart  and 
o  r 

Martinelli  reported  the  relationship  between  $   and  X  graphically  by  ana- 
lyzing the  pressure  drop  data  of  simultaneous  adiabatic  flow  of  air  and 
various  liquids.   Soliman  et  al.  [10]  approximated  the  graphical  results 


of  $  versus  X   by 
v         tt  J 


[1  +  2.85  X 


0.523, 
'tt    i 


(5-15) 


The  frictional  pressure  gradient  in  the  Dukler  II  correlation  [16]  is 
given  by 
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where 


f  =  single  phase  friction  coefficient  evaluated  at  the  two  phase 
Reynolds  number 
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(5-17) 

(5-18) 
(5-19) 

(5-20) 
(5-21) 

(5-22) 


The  correlations  of  Lockhart-Martinelli  and  Dukler  were  originally 
developed  for  adiabatic  two-phase  flow.   For  the  condensation  process, 
Luu  and  Bergles  [5,45]  modified  the  friction  coefficients  in  Eq .  (5-13) 
and  Eq.  (5-16)  by  using  modifiers  suggested  by  Silver  and  Wallis  [49]. 
They  replaced  the  friction  coefficient  f  by  f   and  they  were  related  by 
(6/2fJ 


CO   0 


where 
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(5-23) 


(5-24) 


The  correlation  of  Lockhart-Martinelli  and  Dukler  II  are  based  on  the 

separated  two  phase  flow  model.   The  momentum  pressure  gradient  based  on 
the  same  model  as  outlined  by  Said  [2]  is  given  by 

,2        2 


2.   The  momentum  pressure  gradient  |-r*- 
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To  calculate  the  momentum  pressure  gradient,  one  needs  to  know  the  void 


dx 


fraction  i)J  and  the  quality  gradient  \^~\  ■      In  the  present  study,  the  void 
fraction  correlations  of  the  homogeneous  model  [18],  Lockhart-Martinelli 
[14],  Baroczy  [19],  Thorn  [20],  Zivi  [21],  and  Turner-Wallis  [22]  were 
selected.   Butterworth  [50]  suggested  a  generalized  expression  for  i|i 
given  by 
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(5-26) 


The  constants  A   q   r  ,  and  S  are  listed  in  Table  5.1  for  various  cor- 
relations. 

By  differentiating  Eq.  (5-26)  with  respect  to  x  and  substituting  in 
Eq.  (5-25)  the  following  expression  is  obtained. 
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Table  5.1  The  Constants  in  Eq.  (5-26)  for  the  Different  Void 
Fraction  Correlations 


Model 


Homogeneous  [18] 

1.0 

1.0 

1.0 

0 

Lockhart-Martinelli 

[14] 

0.28 

0.64 

0.36 

0.07 

Baroczy  [19] 

1.0 

0.74 

0.65 

0.13 

Thorn  [20] 

1.0 

1.0 

0.89 

0.18 

Zivi  [21] 

1.0 

1.0 

0.67 

0 

Turner-Wallis  [22] 

1.0 

0.72 

0.40 

0.08 
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To  predict  the  total  pressure  drop  during  condensation,  the  following 

procedure  was  adopted. 

a)  The  test  condenser  was  divided  into  a  number  of  small  increments  (80 
increments)  of  constant  lengths. 

b)  A  linear  change  in  quality  was  assumed  for  which  \~r~\  =  (-1/L)  where 
L  =  total  length  of  test  condenser. 

c)  Since  the  total  pressure  drop  involved  was  very  small,  the  properties 
of  liquid  phase  were  calculated  at  inlet  pressure  conditions. 

d)  The  frictional  pressure  gradient  was  calculated  separately  using 
Lockhart-Martinelli  and  Dukler  II  correlations,  and  for  momentum 
using  Eq.  (5-27)  for  various  void  fractions,  for  each  incremental 
leng  th . 

e)  The  pressure  drop  for  each  incremental  length  was  calculated  by  mul- 
tiplying the  pressure  gradients  of  friction  and  momentum  with  incre- 
mental length. 

f)  By  repeating  steps  (d)  and  (e)  for  all  incremental  lengths,  the  total 
pressure  drop  for  the  entire  test  condenser  tube  was  calculated.  The 
above  scheme  was  implemented  in  a  digital  computer  program. 

As  mentioned  in  Chapter  IV,  the  quality  at  exit  of  each  subsection 
could  be  estimated  by  energy  balance  of  the  condensing  fluid  and  cooling 
water.   This  information  can  be  used  to  find  MM  for  each  subsection  to 
calculate  the  momentum  pressure  gradient  from  Eq.  (5-27).   Said  [2]  used 
both  methods  of  quality  variation  to  calculate  the  momentum  and  hence,  the 
total  pressure  drop  in  the  test  condenser  tube.   He  concluded  that  the 
best  agreement  between  the  predicted  and  the  measured  pressure  drops  were  ob- 
tained by  assuming  the  quality  gradient  to  be  equal  to  (-1/L)  for  the  entire 
test  condenser  tube.   This  was  the  approach  followed  in  calculating  the 
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pressure  drop  m  the  present  study. 

In  the  present  study,  all  possible  combinations  of  the  frictional 
pressure  gradients  of  Lockhart-Martinelli  [14],  and  Dukler  II  [16]  with 
and  without  the  modification  of  Eqs.  (5-23)  and  (5-24),  and  the  momentum 
pressure  gradients  based  on  the  void  fraction  correlations  of  the  homo- 
geneous model  [18],  Lockhart-Martinelli  [14],  Baroczy  [19],  Thorn  [20], 
Zivi  [21],  and  Turner  and  Wallis  [22],  were  tried  to  predict  the  pressure 
drop  of  the  experimental  runs  of  tube  1.   The  best  agreement  between  the 
experimental  pressure  drop  and  the  predictions  were  obtained  without  apply- 
ing the  modified  friction  coefficient  given  by  Eq.  (5-23).   In  general,  the 
Lockhart-Martinelli  [14]  frictional  pressure  drop  correlation  predicted 
reasonably  well  the  experimental  pressure  drop  measurements  with  any  of  the 
six  void  fraction  correlations  used.   The  pressure  drop  predictions  using  the 
Dukler  II  [16]  frictional  pressure  drop  correlation  was  relatively  sensi- 
tive to  the  choice  of  the  void  fraction  correlations.   Only,  the  void 
fraction  correlation  of  Zivi  in  combination  with  the  Dukler  II  frictional 
correlation  predicted  well  the  experimental  overall  pressure  drop  results. 

Figure  5.5  shows  a  comparison  between  the  experimental  overall  pres- 
sure drop  and  the  predictions  of  the  combination  of  Lockhart-Martinelli 
correlation  for  frictional  pressure  drop  and  the  void  fraction  correlation 
of  Zivi.   The  results  show  that  the  predictions  were  within  ±40%  from  the 
experimental  overall  pressure  drop  for  64%  of  data  points.   Figure  5.6  shows 
a  comparison  between  the  experimental  and  predicted  overall  pressure  drops 
for  the  combination  of  Dukler  II/Zivi  correlations.   53%  of  data  points  of 
the  predicted  values  in  Fig.  5.6  were  within  ±40%  from  the  experimental 
overall  pressure  drop  measurements. 

It  is  of  interest  to  point  out  that  Royal  and  Bergles  [4,43]  reported 
that  the  best  predictor  of  their  pressure  drop  measurements  of  steam  was 
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Fig  5.5  Comparison  of  experimental  overall  pressure  drop  with 
predictions  of  the  correlations  of  Lockhart-Martinelli 
^14  7  for  friction  and  Zivi  [  21J   for  void  fraction,  tube  1 
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Fig  5.6  Comparison  of  experimental  overall  pressure  drop  with 

predictions  of  the  correlations  of  Dukler  II  f  16  J   for 
friction  and  Zivi  f 21 J  for  void  fraction,  tube  1 


the  combination  of  the  Dukler  II  [16]  and  the  Hughmark  [23]  void  fraction 
correlation.   Luu  and  Bergles  reported  that  Lockhart-Martinelli  [14]  and 
Dukler  II  [16]  correlations  for  frictional  pressure  drop  and  the  homo- 
geneous [18]  void  fraction  correlation  for  momentum  pressure  drop  predicted 
their  data  best.   They  also  reported  that  the  calculated  momentum  pressure 
drops  were  insensitive  to  the  void  fraction  values  predicted  by  different 
void  fraction  correlations.   Hence  they  used  the  homogeneous  void  fraction 
correlation  for  simplicity.   Said  and  Azer  [2,47]  tested  their  data  with 
various  frictional  and  momentum  pressure  drop  correlations  and  concluded 
that  the  Lockhart-Martinelli  correlation  [14]  for  frictional  pressure  drop 
and  the  homogeneous  [18]  void  fraction  for  momentum  pressure  drop  predicted 
their  smooth  tube  data  best. 

5.3  Internally  Finned  Tubes 

All  the  internally  finned  tubes  tested  gave  higher  heat  transfer  co- 
efficients compared  to  the  smooth  tube.   Also,  internally  finned  tubes 
yielded  higher  overall  pressure  drops  compared  to  the  smooth  tube  results. 
These  results  may  be  attributed  to  the  altered  flow  mechanism  in  in-tube 
condensation  due  to  the  presence  of  internal  fins  besides  the  increase  in 
heat  transfer  area.   The  performance  evaluation  of  the  internally  finned 
tubes,  taking  into  consideration  the  increased  heat  transfer  coefficients 
at  the  expense  of  the  overall  pressure  drop ,  will  be  made  in  the  next 
chapter. 

5.3.1  Heat  Transfer 
Finned  tube  correlations. 

The  development  of  internally  finned  tube  heat  transfer  correlations 
were  based  upon  applying  modifying  factors  to  the  smooth  tube  correlations 
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to  bring  about  the  best  agreement  between  the  predictions  and  the  experi- 
mental results.   Various  modifiers  based  on  certain  geometric  parameters 

i 
of  the  finned  tube  such  as  (W/D,  )  ,  (H/D,  ) ,  and  b  /WD,)  were  introduced 

h       h  1 

by  various  investigators. 

The  following  correlations  have  been  suggested  for  in-tube  conden- 
sation inside  internally  finned  tubes. 

The  correlation  of  Royal  and  Bergles  [42]  is  given  by 


h=0.0265  (k£/Dh) 
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where 
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(5-29) 


Royal  and  Bergles  used  this  correlation  to  correlate  their  heat  transfer 
data  of  steam  condensed  inside  horizontal  tubes  with  straight  and  spiral 
internal  fins.   This  correlation  was  based  on  applying  a  modifying  factor, 

9 

which  is  a  function  of  (b~/WE.),  to  the  smooth  tube  correlations  of  Akers 
et  al.  [6]. 

Luu  and  Bergles  [46]  proposed  the  following  correlation  which  best 
predicted  their  R-113  condensation  data  inside  internally  finned  tubes. 
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(5-31) 


This  correlation  was  based  on  modifying  the  smooth  tube  correlation  of 
Boyko  and  Kruzhilin  [7]  by  the  factor  (b  /WD,)     . 

Said  and  Azer  [47]  introduced  the  factors  F   F„,  and  F  as  modifiers 
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to  the  smooth  tube  correlations  of  Akers  et  al.  [6]  to  get  the  best  agree- 
ment between  the  experimental  and  predicted  values  of  the  condensation 
heat  transfer  coefficients.   Their  correlation  is  given  bv 


h=0.0265(k  /D  ) 


'G  D. 
e  1 


0.8 


Pr°-333  [1  +  0.93  F°-23  F°-58  F^17  Re°-054] 
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where 


G  =  G[(l  -  x)  +  x 


P*l0-5, 


(5-33) 


and 


Re^   =  G  D±(l  -  x)/w4 


(5-34) 


The  geometric  factors  F   F_,  and  F,  were  originally  introduced  by 
Carnavos  [28,29]  in  his  study  of  single  phase  heat  transfer  and  friction  co- 
efficient of  finned  tubes  to  get  the  best  agreement  between  the  experimental 
measurements  and  the  predictions.   These  modifying  factors  are  defined  as 
follows. 


Fn  =  A,  /A, 
1    fa  fc 


F,  =  A  /A 
2    n  a 


F,  =  sec  a 
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where 


ttD 


A   =  actual  free  flow  area,  — 7—   ,  cm 
ta  4 

<  2 

A   =  open  core  free  flow  area,  —7—  ,  cm 

A   =  nominal  heat  transfer  area  based  on  D.  as  if  fins  were 
n  9  ,  1 

not  present,  cm^/cm 


actual  heat  transfer  area,  cm*"/cm 
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The  equivalent  diameter  D  and  core  diameter  D  are  defined  bv 
e  c  J 

2  2 

ttD  ttDT 

-A  =  —T-  -  nbt/cosa  (5-38) 

~t  - !  (Di  - 2b)  (5"39) 

Hydraulic  diameter  is  defined  by 


D.  =  4A.  /A  (5-40) 

□     ra  a 


With   the  above   definitions,    it  can  be   shown   that 
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The  values  of  F  ,  F,,  and  F,  for  the  tubes  tested  in  the  present  study  are 
given  in  Table  5.2.   Also  the  various  geometric  parameters  of  the  finned 
tubes  are  shown  in  Fig .  5.7. 


Table  5.2  Computed  Values  of  F   F»,  and  F 


F 

Tub 

e  Tested 

Value 

2 

3 

4 

5 

Fl 

1.515 

1.514 

1.464 

1.355 

F2 

0.6678 

0.5281 

0.5667 

0.6085 

F3 

1.0 

1.0145 

1.0236 

1.0599 

So 


Fig.  5.7  Cross  Section  of  Internally  Finned  Tube. 
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Figures  5.8  through  5.16  show  comparisons  between  the  predictions  of 

Royal  and  Bergles  correlation  [42]  with  experimental  measurements. 
Figures  5.8,  5.10,  5.12  and  5.14  show  overall  average  heat  transfer  co- 
efficients.  Except  for  tube  3,  Fig.  5.10,  the  correlation  underpredicted 
h  for  tubes  2,  4,  and  5.   Similarly  Figs.  5.9,  5.11,  5.13  and  5.15  show 
sectional  averages  and  overall  heat  transfer  coefficients  for  tubes  2,  4 
and  5,  respectively.   Except  for  tube  3,  all  the  sectional  heat  transfer 
coefficients  were  underpredicted  by  the  correlation.   For  tube  3,  as  shown 
in  Fig.  5.11,  the  correlation  overpredicted  the  sectional  averages  of 
section  1  and  underpredicted  the  sectional  averages  of  section  4.   A  com- 
parison of  the  predictions  of  Royal  and  Bergles  with  the  experimental 
values  for  all  the  finned  tubes  is  shown  in  Fig.  5.16.   It  is  clear  from 
the  figure  that  the  correlation  underpredicted  h  for  all  finned  tubes 
except  tube  3. 

Figures  5.17  through  5.25  show  comparisons  between  the  experimental 
measurements  and  the  predictions  of  Luu  and  Bergles  correlation  [46]  for 
the  heat  transfer  coefficients.   Figures  5.17,  5.19,  5.21,  and  5.23  show 
comparisons  for  the  overall  heat  transfer  coefficients.   The  correlation 
predicted  h  within  ±30%  from  experimental  values  for  79%  of  the  data  points 
for  tube  2,  100%  for  tubes  3  and  4.   The  correlation  underpredicted  h 
for  tube  5.   Figures  5.18,  5.20,  5.22  and  5.24  show  the  comparisons  of 
the  experimental  sectional  averages  and  the  overall  heat  transfer  coef- 
ficients with  the  predictions.   The  correlation  had  a  tendency  to  over- 
predict  the  sectional  average  heat  transfer  coefficients  for  section  1, 
except  for  tube  5  where  they  were  underpredicted.   For  tubes  3,  4,  and  5, 
the  correlation  tended  to  underpredict  the  sectional  average  heat  transfer 

coefficient  for  section  4.   The  larcp  sfaft»r  <„  fk„ 
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Fig  5.8  Comparison  of  experimental  overall  heat  transfer 

coefficients  with  predictions  of  the  correlation  of 
Royal  and  Bergles  £42  J,  Eq.(5-28),  tube  2 
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Fig  5.9   Comparison  of  experimental  sectional  and  overall  heat 

transfer  coefficients  with  predictions  of  the  correlation 
of  Royal  and  Bergles  £  42  J,    Eq.(5-28),  tube  2 
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Fig  5.10  Comparison  of  experimental  overall  heat  transfer 

coefficients  with  predictions  of  the  correlation  of 
Royal  and  Bergles  [  42  J ,  Eq.(5-28),  tube  3 
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Fig  5.11   Comparison  of  experimental  sectional  and  overall  heat 

transfer  coefficients  with  predictions  of  the  correlation 
of  Royal  and  Bergles  [  42  J ',  Eq.(5-28),  tube  3 
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Fig  5.12   Comparison  of  experimental  overall  heat  transfer  coefficients 
with  predictions  of  the  correlation  of  Royal  and  Bergles 
[  <*2  J,    Eq.(5-28),  tube  4 
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5.13  Comparison  of  experimental  sectional  and  overall  heat 

transfer  coefficients  with  predictions  of  the  correlation 
of  Royal  and  Bergles  [  42  J ,    Eq.(5-28),  tube  4 
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Fig  5.14  Comparison  of  experimental  overall  heat  transfer  coefficients 
with  predictions  of  the  correction  of  Royal  and  Bergles 
/"42J,  Eq.(5-28),  tube  5 
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Fig  5.15   Comparison  of  experimental  sectional  and  overall  heat 

transfer  coefficients  with  predictions  of  the  correlation 
of  Royal  and  Bergles  [  42  J ,    Eq.(5-28),  tube  5 
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Fig  5.16   Comparison  of  experimental  overall  heat  transfer 

coefficients  with  predictions  of  the  correlation  of 
Royal  and  Bergles  [  42  J ,  Eq.(5-28),  all  finned  tubes 
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Fig   5.17      Comparison  of   experimental   overall  heat    transfer   coefficients 
with  predictions   of    the  correlation     of  LuU  and   Bergles 
CkbJ,    Eq.(5-30),    tube   2 
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Fig  5.18   Comparison  of  experimental  sectional  and  overall  heat  transfer 
coefficients  with  predictions  of  the  correlation  of  Luu  and 
Bergles  Z"46j,  Eq.(5-30),  tube  2 
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Fig  5.19  Comparison  of  experimental  overall  heat  transfer  coefficients 
with  predictions  of  the  correlation  of  Luu  and  Bergles 
f  46  J,  Eq.(5-30),  tube  3 


100 


I 


X 


B) 

O 

-a 


\ 

3 


a 


a      1 


exp 

IB 


W/m    -deg    c     0x10       ) 

IS  2C  25 


©  overall" 

*  1  section         aa 

D  2    SECTION    ±AiP* 

*  3    SECTION 
+  4    SECTION      *      4T^    C&f, 


30 


8 
X 


□J 

£ 
\ 

2 


18 


exp 


Btu/ft    -hr-deg    F     ( x 1 0 ~d ) 


Fig  5.20  Comparison  of  experimental  sectional  and  overall  heat 

transfer  coefficients  with  predictions  of  the  correlation 
of  Luu  and  Bergles  [ 46  J  ,    Eq.(5-30),  tube  3 
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Fig  5.21   Comparison  of  experimental  overall  heat  transfer  coefficients 
with  predictions  of  the  correlation  of  Luu  and  Bergles 
f  46.7,  Eq.(5-30),  tube  4 
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Fig  5.22   Comparison  of  experimental  sectional  and  overall  heat 

transfer  coefficients  with  predictions  of  the  correlation 
of  Luu  and  Bergles  [  46  J ,  Eq.(5-30),  tube  4 
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Fig  5.23  Comparison  of  experimental  overall  heat  transfer  coefficients 
with  predictions  of  the  correlation  of  Luu  and  Bergles 
T46J,  Eq.(5-30),  tube  5 
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Fig  5.24   Comparison  of  experimental  sectional  and  overall  heat 

transfer  coefficients  with  predictions  of  the  correlation 
of  Luu  and  Bergles  [  46  J ,    Eq.(5-30),  tube  5 
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Fig  5.25  Comparison  of  experimental  overall  heat  transfer  coefficients 
with  predictions  of  the  correlation  of  Luu  and  Bergles 
/"46J,  Eq.(5-30),  all  finned  tubes 
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coefficients  of  section  1  may  be  due  to  the  method  of  their  calculation 
adopted,  as  discussed  earlier.   In  general  Luu  and  Bergles  correlation 
predicted  h  values  fairly  well  for  all  the  tubes  as  shown  in  Fig.  5.25. 

Figures  5.26  through  5.34  show  comparisons  between  the  experimental 
heat  transfer  coefficients  with  the  predictions  of  the  correlation  of 
Said  and  Azer  [47].   The  correlation  underpredicted  the  heat  transfer  co- 
efficients for  all  the  tubes.   Figure  5.27,  5.29,  5.31  and  5.33  are  plots 
of  predicted  h  versus  experimental  measurements  for  tubes  2,  3,  4,  and  5, 
respectively.   Figure  5.35  includes  comparison  of  predicted  h  of  all  the 
finned  tubes  versus  the  experimental  values.   It  is  evident  from  Figs. 
5.26,  5.28,  5.30,  5.32,  and  5.34  that  the  correlation  underpredicted  the 
h  values.   Figures  5.27,  5.29,  5.31,  and  5.33  show  comparison  between  the 
experimental  and  the  predicted  sectional  and  overall  average  heat  transfer 
coefficients  for  all  tubes  tested.   It  is  interesting  to  note  that  although 
Said  and  Azer  [2,47]  developed  their  correlation  from  condensing  R-113 
inside  the  same  internally  finned  tubes  tested  in  the  present  study,  their 
correlation  underpredicted  the  measurements  for  R-ll.   It  should  be 
mentioned  that  the  range  of  flow  rates  with  which  Said  [2]  conducted  his 
experiments  with  R-113  was  higher  than  the  flow  rates  covered  in  the 
present  study. 

5.3.2  Pressure  Drop 
Finned  Tube  Correlations 

Royal  and  Bergles  [43]  had  a  limited  success  in  predicting  their  pres- 
sure drop  measurements  of  condensation  of  steam  inside  internally  finned 
tubes,  by  merely  replacing  the  inside  nominal  diameter  by  the  hydraulic 
diameter  of  the  finned  tube  in  the  Dukler  II  [16]  frictional  pressure 


hexp     W/m    -deS    C     (xl0~2  ) 


107 


I 


S 
x 


33 

D 
"0 
I 
L 

X 


\ 

3 

3 


"'""'""'""''■■■'nnh 


'■''•• 


hexp     Btu/ft    -hr-dcg    F     Cxl0 


-2 


Fig  5.26 


Comparison  of  experimental  overall  heat  transfer  coefficients 
with  predictions  of  the  correlation  of  Said  and  Azer  f47J, 
Eq.(5-32),  tube  2 
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Fig  5.27   Comparison  of  experimental  sectional  and  overall  heat 

transfer  coefficients  with  predcitions  of  the  correlation 
of  Said  and  Azer  [  47  J ,  Eq.(5-32),  tube  2 
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Fig  5.28  Comparison  of  experimental  overall  heat  transfer  coefficients 
with  predictions  of  the  correlation  of  Said  and  Azer  [klj , 
Eq.(5-32),  tube  3 
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Fig  5.29   Comparison  of  experimental  sectional  and  overall  heat  transfer 
coefficients  with  predictions  of  the  correlation  of  Said  and 
Azer  f  477,  Eq.(5-32),  tube  3 
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Fig  5.30  Comparison  of  experimental  overall  heat  transfer  coefficients 
with  predictions  of  the  correlation  of  Said  and  Azer  [  47  J, 


Eq.(5-32),  tube  4 
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Fig  5.31   Comparison  of  experimental  sectional  and  overall  heat  transfer 
coefficients  with  predictions  of  the  correlation  of  Said  and 
Azer  [ 47 J ,    Eq.(5-32),  tube  4 
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Fig  5.32   Comparison  of  experimental  overall  heat  transfer  coefficients 
with  predictions  of  the  correlation  of  Said  and  Azer  [  Wl  J, 
Eq.(5-32),  tube  5 
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Fig  5.33  Comparison  of  experimental  sectional  and  overall  heat  transfer 
coefficients  with  predictions  of  the  correlation  of  Said  and 
Azer  C  1*7  J,    Eq.(5-32),  tube  5 
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Fig  5.34   Comparison  of  experimental  sectional  and  overall  heat  transfer 
coefficients  with  predictions  of  the  correlation  of  Said  and 
Azer  [  47  J,    Eq.(5-32),  all  finned  tubes 
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drop  correlation  along  with  the  homogeneous  [18]  void  fraction  correlation. 

Luu  and  Bergles  [5,45,46]  applied  the  modification  of  Eq .  (5-23)  to  the 
friction  coefficient  in  the  Dukler  II  correlation  [16].  They  were  able 
to  correlate  reasonably  well  their  pressure  drop  data  of  R-113  with  the 
combination  of  the  modified  Dukler  correlation  and  the  homogeneous  model 
of  void  fraction.  They  used  also  the  hydraulic  diameter  in  the  Dukler 's 
correlation.  Said  and  Azer  [2,47]  introduced  the  following  modifier  to 
the  friction  coefficient  of  the  Dukler  correlation. 

3  37   0  521 
ffc  =  fo  (1+  F3    F4    >  (5'43) 

f,      replaces   f      in  Eq.    (5-16)    and   f      is   calculated  by  Eq .    (5-17).      F, 
tc  o  o  -n*'4 

is   a   factor  given  by: 


F.    -  A,    /A, 

4  fa      fn 


fDl    2 


D. 


[1  -    (4nbt)/(nD2   cosa)]  (5-44) 


A,  =  nominal  flow  area  based  on  D.  as  if  fins  were  not  present. 

fn     ,  l  r  ' 

*Di     2 

— ; —  ,  cm 
4 

F,  was  introduced  by  Carnavos  [28,29]  while  correlating  the  pressure 

drops  of  single  phase  flow  inside  internally  finned  tubes.   Table  5.3 

gives  the  value  of  F,  for  the  tubes  tested. 

Table  5.3  Computed  Values  of  F, 


F 

Tube   Tested 

Value 

2 

3                          4 

5 

F4 

0.9174 

0.9442               0.9502 

0.9376 
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Said  and  Azer  retained  D.  rather  than  replacing  it  by  D  in  the 


Dukler  II  [16]  correlation. 

After  careful  examination  of  these  existing  finned  tube  pressure 
drop  correlations,  it  was  decided  to  test  the  following  correlations  to 
predict  pressure  drop  in  the  present  study. 

A.  Lockhart-Martinelli  frictional  pressure  drop  [14]  correlation  with 
D  replacing  D.  along  with  a  combination  of  various  void  fraction  cor- 
relations tried  for  smooth  tube  predictions.   Lockhart-Martinelli  frictional 
pressure  drop  correlation  was  tried  with  and  without  the  friction  coef- 
ficient modification  of  Eq.  (5-23). 

B.  Dukler  II  [16]  frictional  pressure  drop  correlation  with  D  replacing 
D.  along  with  a  combination  of  various  void  fraction  correlations.  Dukler 
II  friction  pressure  drop  correlation  was  tried  with  and  without  the  fric- 
tion coefficient  modification  of  Eq.  (5-23). 

C.  Dukler  II  [16]  friction  pressure  drop  correlation  with  friction  coef- 
ficient f  calculated  from  Eq.  (5-43)  as  suggested  by  Said  and  Azer  [2,47] 
along  with  a  combination  of  various  void  fraction  correlations. 

The  results  of  these  correlations  are  summarized  in  the  following. 

1.  Lockhart-Martinelli  correlation  [14]  with  any  of  the  void  fraction 
correlations  predicted  the  measured  pressure  drop  fairly  well.   The 
introduction  of  friction  coefficient  modification  of  Eq.  (5-23)  in 
Lockhart-Martinelli  frictional  pressure  drop  correlation  overpredicted 
the  measured  pressure  drop. 

2.  Dukler  II  correlation  [16]  for  frictional  pressure  drop  predictions 
was  relatively  sensitive  to  the  void  fraction  correlation.   Only 
Zivi's  [21]  void  fraction  correlation  with  Dukler  II  correlation 
for  friction,  predicted  the  measured  pressure  drop  well.   The 


118 
modification  of  friction  coefficient  of  Eq.  (5-23)  yielded  better 
predictions  with  Dukler  II  correlation. 

3.  Dukler  II  [16]  fricLical  pressure  correlation  with  friction  coef- 
ficient f,   for  finned  tubes  as  suggested  by  Said  and  Azer  [2,47] 
in  combination  with  Zivi  void  fraction  correlation  predicted  the 
measured  pressure  drop  of  tube  3  within  ±40%  from  experimental 
values,  but  overpredicted  the  pressure  drop  for  tube  2. 

4.  As  in  smooth  tube,  momentum  pressure  drop  was  insensitive  to  the  void 
fraction  correlation  used. 

Figures  5.35,  5.36,  and  5.37  show  comparisons  between  the  predictions 
of  the  Lockhart-Martinelli/Zivi  void  and  experimental  measurements.   The 
results  show  that  58%  of  predicted  values  for  tube  2,  and  all  the  pre- 
dicted values  for  tube  3  are  within  ±40%  from  experimental  pressure  drop 
measurements.   Figure  5.37  is  a  plot  of  the  predicted  values  versus  experi- 
mental measurements  for  tubes  2  and  3.   Figures  5.38,  5.39,  and  5.40  show 
comparisons  between  modified  Dukler  II/Zivi  void  correlation  and  experimental 
pressure  drops.   The  results  show  that  79%  of  data  points  of  tube  2  and  94% 
of  data  points  of  tube  3  agree  within  ±40%  with  the  predicted  values.   It  is 
clear  from  Fig.  5.40  that  the  modified  Dukler  II/Zivi  void  correlation  pre- 
dicted the  finned  tube  condensation  pressure  drop  best.   Figures  5.41,  5.42, 
and  5.43  compares  predicted  pressure  drops  of  Dukler  II/Zivi  with  Said  and 
Azer  [2,47]  friction  coefficient  of  Eq.  (5-43),  and  experimental  measure- 
ments.  The  results  show  that  the  correlation  predicted  pressure  drops  of 
tube  3  well,  but  overpredicted  the  pressure  drop  of  tube  2.   Figure  5.44 
shows  comparison  between  the  predictions  and  experimental  measurements 
for  tubes  2  and  3. 
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Fig   5.35      Comparison  between  experimental   pressure  drop   and  predictions 
of   the  combination  of  Lockhart-Martinelli   [  14  J   and   Zivi 
[  21  J  void   fraction  correlations,    tube   2 
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Fig  5.36  Comparison  between  experimental  pressure  drop  and  predictions 
of  the  combination  of  Lockhart-Martinelli  [  14  J   and  Zivi 
C  21  J  void  fraction  correlations,  tube  3 
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Fig  5.37   Comparison  between  experimental  pressure  drop  and  predictions 
of  the  combination  of  Lockhart-Martinelli  f  14  J   and  Zivi 
f  21  J   void  fraction  correlations,  tube  2  and  3 
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Fig  5.38  Comparison  between  experimental  pressure  drop  and  predictions 
of  the  combination  of  Dukler  II  f 16  J   correlation,  modified 
by  Eq.(5-23)  as  suggested  by  Luu  and  Bergles,  and  Zivi  [  21  ] 
void  fraction,  tube  2 
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Fig  5.39   Comparison  between  experimental  pressure  drop  and  predictions 
of  the  combination  of  Dukler  II  £  16  J   correlation,  modified 
by  Eq.  (5-23), as  suggested  by  Luu  and  Bergles,  and  Zivi  [21] 
void  fraction,  tube  3 
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Fig  5.40  Comparison  between  experimental  pressure  drop  and  predictions 
of  the  combination  of  Dukler  II  [  16  7  correlation,  modified 
by  Eq.(5-23),  as  suggested  by  Luu  and  Bergles,  and  Zivi  [21] 
void  fraction,  tubes  2  and  3 


125 


A  P. 


a. 

< 


A  P 


cxp  PS' 


Fig  5.41   Comparison  between  experimental  pressure  drop  and  predictions 
of  the  combination  of  Dukler  II  f 16  J  correlation,  modified 
by  Eq.(5-43),  as  suggested  by  Said  and  Azer,  and  Zivl  [  21  J 
void  fraction,  tube  2 
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Fig  5.42  Comparison  between  experimental  pressure  drop  and  predictions 
of  the  combination  of  Dukler  II  [  16  J   correlation,  modified 
by  Eq.(5-43),  as  suggested  by  Said  and  Azer,  and  Zlvl  [  21 J 
void  fraction,  tube  3 
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Fig  5.43  Comparison  between  experimental  pressure  drop  and  predictions 
of  the  combination  of  Dukler  II  [   16  J   correlation,  modified 
by  Eq.(5-43),  as  suggested  by  Said  and  Azer,  and  Zivi  [  21J 
void  fraction,  tube  2  and  3 
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Chapter  VI 
PERFORMANCE  EVALUATION 

6,1  Introduction 

The  primary  objective  of  this  part  of  the  study  was  to  evaluate  the  per- 
formance of  the  internally  finned  tubes  in  agumenting  in-tube  condensation 
of  R-ll.   For  meaningful  performance  evaluation,  the  increase  in  heat  trans- 
fer coefficients  with  internally  finned  tubes  should  be  considered  along 
with  the  penalty  of  increased  pressure  drop.   In  order  to  compare  the 
increase  in  heat  transfer  coefficients  and  pressure  drops,  the  overall 
heat  transfer  coefficients  and  pressure  drops  versus  mass  flux  were  combined 
in  composite  plots  for  same  inlet  pressures,  as  shown  in  Figures  6.1,  through 
6.5,  for  all  the  tubes  tested.   Figures  6.1,  6.2  and  6.3  show  plots  of  the 
regression  equations  of  the  experimental  overall  average  heat  transfer  co- 
efficients h  versus  mass  flux  G  for  all  tubes  at  the  three  different  nominal 
pressures  tested.   The  best  performer  among  the  finned  tubes  was  tube  2 
followed  by  tube  3.   Due  to  the  limitation  of  the  test  facility,  it  was  not 
possible  to  cover  the  same  range  of  mass  flux  for  all  tubes.   As  the  inside 
diameter  of  the  test  condenser  increased,  the  range  of  mass  flux  decreased  as 
indicated  in  Figs.  6.1  through  6.5.   On  a  nominal  area  basis,  tubes  2  and  3 
enhanced  the  heat  transfer,  on  the  average,  by  55%  and  13%  respectively. 
Tube  5  performed  better  than  tube  4.   Since  the  range  of  mass  flux  for 
tubes  4  and  5  were  less  than  tube  1,  no  direct  comparison  could  be  made. 
Figures  6.4  and  6.5  show  plots  of  the  regression  lines  of  the  experimental 
pressure  drop  versus  the  mass  flux  for  tubes  1,  2,  and  3.   The  results  show 
that  tube  3  had  the  highest  pressure  drop  followed  by  tube  2.   So,  it  is 
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obvious  that  heat  transfer  was  enhanced  with  the  penalty  of  increasing 
the  pressure  drop. 

6.2  Performance  Indices 

There  have  been  numerous  attempts  to  develop  indices  to  evaluate  the 
performance  of  augmented  surfaces.   Bergles  et  al.  [51,52]  suggested  nine 
performance  indices  for  augmentation  of  heat  transfer  in  single  phase  flow. 
They  listed  the  objectives  of  augmentation  as  the  increase  in  heat  transfer, 
reduction  in  pumping  power,  and  reduction  in  size  of  equipment.   They  also 
listed  the  controlling  parameters  with  which  these  objectives  can  be 
achieved  such  as  basic  geometry,  flow  rate,  pressure  drops,  pumping  power, 
and  heat  duty.   In  developing  the  indices,  they  broadly  divided  the  per- 
formance indices  into  two  basic  groups.   In  the  first  group,  the  performance 
indices  were  directed  towards  enhancement  of  heat  transfer  for  existing  heat 
exchangers.   Therefore^ the  basic  geometry  was  fixed.   In  the  second  group, 
the  performance  indices  were  directed  towards  heat  transfer  augmentation 
in  the  design  of  new  heat  exchangers  with  the  objective  of  reducing  their 
size. 

Due  to  the  difficulty  of  imposing  some  of  the  constraints  on  the  per- 
formance indices,  Royal  [4]  used  two  of  these  indices  to  evaluate  the 
performance  of  internally  finned  tubes  and  tubes  with  twisted  tapes  during 
augmentation  of  condensation  heat  transfer  of  steam.   These  indices  were 
the  condenser  size  reduction  index  R,  and  the  pressure  drop  index  R   ,  and 
they  are  given  by 

R,  =  (A   /A  )  =  (h  /h   )  (6-1) 

n     aug  sm      sm  aug  *■   ' 

and 

R,   =  (AP    /AP   )  (6_7N 

Ap      aug   sm  *•"  *■' 
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In  the  above  indices  it  was  assumed  that  the  external  resistance  of 
condensing  surface  is  negligible  or  the  condensing  side  thermal  resistance 
controls  the  heat  transfer.   Also,  ft  and  R   were  evaluated  under  the  con- 
straints of  fixed  heat  duty,  nominal  diameter,  and  constant  temperature 
difference. 

Luu  [5]  used  the  same  indices  of  Eq.  (6-1)  and  Eq.  (6-2)  in  evaluating 
the  performance  of  internally  finned  tubes  and  tubes  with  twisted  tape 
inserts  in  augmenting  in-tube  condensation  of  R-113. 

In  the  present  study,  it  was  estimated  that  the  outside  heat  transfer 
coefficients  were  of  the  same  order  of  magnitude  of  the  inside  heat  transfer 
coefficient.   Hence  the  assumption  of  zero  external  resistance  is  unreal- 
istic and  the  use  of  performance  indices  of  Eq.  (6-1)  and  Eq.  (6-2)  with 
negligible  external  resistance  could  not  be  justified.   Said  [2]  arrived 
at  the  same  conclusion  while  evaluating  the  performance  of  the  internally 
finned  tubes  with  R-113  as  the  condensing  fluid. 

Azer  et  al.  [53]  and  Lin  [54]  used  the  ratio  of  the  pumping  power  to 
the  rate  of  heat  transfer  P/Q  as  an  evaluating  index  in  determining  the 
performance  of  static  in-line  mixers  while  augmenting  the  condensation 
heat  transfer.   Said  [2]  also  used  the  above  ratio  as  a  performance  index 
in  determining  the  performance  of  internally  finned  tubes  and  tubes  with 
twisted  tape  inserts.   The  index  (F/Q)  was  evaluated  under  the  constraints 
of  fixed  geometry,  same  mass  flow  rate  of  condensing  fluid,  and  same  tempe- 
rature and  mass  flow  rate  of  the  coolant.   The  pumping  power  was  obtained 
from  the  product  of  the  volume  flow  rate  of  the  liquid  at  the  circulating 
pump  and  the  pressure  drop  across  the  test  condenser.   The  rate  of  heat 
transfer  was  evaluated  from  the  temperature  rise  and  flow  rate  of  the 
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coolant.   The  same  index  was  used  in  the  present  study. 

Only  the  smooth  tube,  tube  1,  and  the  straight  finned  tube,  tube  2, 
satisfied  the  constraint  of  fixed  geometry,  which  required  that  the  length 
and  the  inside  diameter  of  the  tube,  augmented  or  unaugmented,  be  the  same. 
Therefore,  only  tubes  1  and  2  were  evaluated. 

The  constraints  of  same  temperature  and  flow  rate  were  replaced  by  the 

constraints  that  the  inlet  pressure  P.   and  the  flow  rate  of  condensing 

fluid  were  the  same.   The  constraints  of  the  constant  P.   and  mass  flow 

in 

rate  imply  a  constant  heat  load  when  condensation  is  complete  with  the  fol- 
lowing assumptions. 

1.  Pressure  change  during  condensation  has  negligible  effect  on  the 
latent  heat. 

2.  The  effect  of  inlet  superheat  and  exit  subcooling  on  total  energy 
transfer  is  negligible. 

After  careful  examination,  experimental  runs  of  tube  1  and  tube  2  with 
the  same  mass  flow  rate  and  the  same  inlet  pressure  were  selected  for  the 
evaluation  of  performance  of  tube  2.   Figures  6.6  and  6.7  show  plots  of 
1P/Q  versus  Re  for  tubes  1  and  2  at  two  nominal  inlet  pressures,  1.32  bar 
(19.19  psia)  and  1.67  bar  (24.17  psia),  respectively.   Re  was  based  on 
the  liquid  flowing  alone,  the  inside  diameter  D.,  and  the  dynamic  viscosity 
of  the  saturated  liquid.   The  lower  the  index,  the  lower  the  power  demand 
per  unit  heat  transfer.   It  is  clear  from  Figs.  6.6  and  6.7  that  power 
demand  for  tube  2  was  more  than  that  of  tube  1.   It  is  obvious  that  enhance- 
ment of  heat  transfer  was  accompanied  by  an  increase  in  pressure  drop. 

Azer  and  Shivakumar  [55]  used  another  performance  index  to  evaluate 
performance  of  internally  finned  tubes  during  saturated  boiling  heat 
transfer.   They  defined  the  performance  index  r  as  the  ratio  of  heat  trans- 
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Fig  6.6  Piping  power  per  unit  heat  transfer  rate  versus 
Reynolds  Number  at  P^-  1.32  bar  (19.19  psia) 
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Fig  6.7  "  Pumping  power  per  unit  heat  transfer  rate  versus 
Reynolds  dumber  at  P  -1. 67  bar  (24.17  psia) 


139 
fer  enhancement  to  the  ratio  of  pressure  drop  increase.   This  ratio  r 
was  expressed  by 

(  h  )    /(h) 

(Ap  )    /  (Ap  )  (6_3) 

aug        sra 

This  index  was  calculated  at  same  constraints  of  Eq.  (6-2).   The  higher  the 
ratio  r,  the  better  the  performance.   Figures  6.8  and  6.9  show  plots  of  r 
given  by  Eq.  (6-3)  versus  Reynolds  number  at  two  nominal  inlet  pressures, 
1.32  bar  (19.19  psia)  and  1.67  bar  (24.17  psia) ,  respectively.   The  plots 
indicate  that  the  values  of  r  are  less  than  1.  The  results  also  reinforce 
the  conclusions  drawn  from  Figs.  6.6,  and  6.7,  i.e.,  the  enhancement  in 
heat  transfer  was  accompanied  by  an  increase  in  the  pressure  drop.   However, 
the  trend  in  the  data  in  Figs.  6.8  and  6.9  indicate  that  r  increases  with 
the  increase  of  Re^  and  it  could  become  higher  than  1.   Because  of  the 
limitation  of  the  test  facility  this  observation  could  not  be  established 
experimentally. 

Besides  the  objectives  of  increasing  the  heat  transfer,  reducing  the 
pumping  power,  and  reducing  the  heat  transfer  area,  other  important  factors 
like  initial  cost,  maintenance  cost,  etc.,  are  to  be  considered  by  designer. 
The  choice  of  an  augmented  versus  a  smooth  tube  for  a  given  application 
must  be  based  on  the  compromise  between  the  improvements  in  heat  transfer 
and  the  penalty  of  increasing  the  pressure  drop. 
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Fig.  6.8  The  ratio  r  of  heat  transfer  enhancement  to  the  pressure  drop 
increase  versus  Reynolds  Number  at  P.  =1.32  bar  (19.19  Psia) 
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Chapter  VII 
SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 

Heat  transfer  and  pressure  drop  data  were  taken  during  condensation  of 
R-ll  inside  five  horizontal  tubes.   One  tube  was  smooth  and  four  were  in- 
ternally finned  with  different  geometric  parameters.   The  heat  transfer  and 
pressure  drop  data  were  compared  with  existing  smooth  tube  and  finned  tube 
correlations. 

The  results  are  summarized  in  the  following. 

1.  All  the  internally  finned  tubes  tested  yielded  higher  heat  transfer 
coefficients  than  the  smooth  tube  results.   Over  the  mass  flux  range 
tested,  an  enhancement  in  the  heat  transfer  coefficients  as  high  as 
55%  over  the  smooth  tube  was  obtained,  on  a  nominal  area  basis.   The 
heat  transfer  enhancement  was  accompanied  by  an  increase  in  the  pres- 
sure drop. 

2.  The  best  predictor  of  the  sectional  and  overall  average  heat  transfer 
coefficients  of  the  smooth  tube  was  the  correlations  of  Akers  et  al. 
[6],  Eq.  (5-1). 

3.  The  best  predictors  of  the  smooth  tube  and  internally  finned  tube 
pressure  drop  were  the  frictional  correlation  of  Lockhart-Martinelli 
[14],  and  Dukler  II  [16] ,  combined  with  the  void  fraction  correlation 
of  Zivl  [21]. 

4.  The  best  predictor  of  the  sectional  and  overall  average  heat  transfer 
coefficients  of  the  internally  finned  tubes  was  the  correlation  of 
Luu  and  Bergles  [5,45,46],  Eq .  (5-30). 


143 

5.   The  pumping  power  per  unit  heat  transfer  subject  to  the  constraints 
of  fixed  geometry,  same  inlet  pressure,  and  flow  rate,  was  used  to 
evaluate  the  performance  of  the  internally  finned  tubes.   The  ratio 
of  pumping  power  per  unit  heat  transfer  for  the  internally  finned  tube 
to  the  same  ratio  of  the  smooth  tube  was  greater  than  one.   This  indi- 
cated that  the  penalty  of  increased  pressure  drop  was  greater  than  the 
benefit  of  enhancement  of  the  heat  transfer  coefficients  for  the  in- 
ternally finned  tube.   A  similar  performance  index  r,  a  ratio  of  heat 
transfer  enhancement  to  the  pressure  drop  increase  yielded  values  less 
than  one  for  internally  finned  tube.  The  trend  of  the  variation 
of  r  with  the  mass  flow  rate  suggested  that  a  ratio  greater  than  1  is 
possible.   This  observation  could  not  be  verified  experimentally 
because  of  the  limitations  of  the  test  facility. 

Recommendation  for  Future  Studies 

1.  As  pointed  out  earlier,  before  commercial  usage  of  the  existing  heat 
transfer  and  pressure  drop  design  correlations,  additional  experimental 
data  for  condensation  inside  internally  finned  tubes  with  different 
geometric  parameters,  and  for  different  condensing  fluids  at  wide 

flow  ranges  are  needed.   What  is  currently  available  can  be  considered 
limited.   Also,  the  experimentation  of  future  studies  should  be  aimed 
at  optimizing  the  fin  geometry  to  obtain  best  heat  transfer  and  pres- 
sure drop  results . 

2.  More  reliable  pressure  drop  correlations  need  to  be  developed  for  con- 
densation inside  smooth  tubes  under  different  flow  patterns.   Once 
such  correlations  are  established,  better  pressure  drop  correlations 
for  internally  finned  tubes   could  also  be  developed. 
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3.   Quantification  of  the  performance  of  internally  finned  tubes  is  yet 

to  be  perfected.   Most  of  the  performance  indices  developed  for  single 
phase  flow  are  not  suited  in  two  phase  flow  as  the  constraints  with 
which  these  indices  are  evaluated  cannot  be  achieved  experimentally. 
More  reliable  and  attainable  performance  criteria  of  augmented  surfaces 
need  to  be  defined  for  two-phase  flow. 


145 
ACKNOWLEDGEMENTS 

I  wish  to  express  my  deepest  appreciation  and  gratitude  to 
Prof.  Nairn  Z,   Azer  for  his  guidance  and  valuable  advice  during  the 
course  of  this  study. 

Thanks  are  extended  to  Prof.  Leonard  E.  Fuller  and  Prof.  Robert  E.  Crank 
for  being  on  my  examination  committee. 

I  wish  to  thank  Dr.  Paul  L.  Miller,  Chairman,  Department  of 
Mechanical  Engineering,  for  providing  financial  support  during  the 
present  investigation . 

Finally,  I  wish  to  thank  my  family  for  their  patience,  understanding 
and  continuous  encouragement  during  my  graduate  study. 


146 


SELECTED  BIBLIOGRAPHY 

L.   Bergles,  A.  E. ,  "Research  Workshop  on  Augmentation  of  Convective  Heat 
Transfer",  Final  Report  HTL-8,  ISU-ERI-AMES-76026  (1975). 

2.  Said,  S.,  "Augmentation  of  Condensation  Heat  Transfer  of  R-113  by 
Internally  Finned  Tubes  and  Twisted  Tape  Inserts",  Ph.D.  dissertation, 
Dept.  of  Mechanical  Engineering,  Kansas  State  University  (1982). 

3.  Bergles,  A.  E. ,  Webb,  R.  L. ,  Junken,  G.  H. ,  and  Jensen,  M.  K. , 
"Bibliography  on  Augmentation  of  Convective  Heat  and  Mass  Transfer", 
Engineering  Research  Institute,  Iowa  State  University,  HTL-19-ISU-AMES-79206 
(1979). 

4.  Royal,  J.,  "Augmentation  of  Horizontal  In-Tube  Condensation  of  Steam", 
Ph.D.  dissertation,  Dept.  of  Mechanical  Engineering,  Iowa  State  University, 
Ames,  Iowa  (1980). 

5.  Luu,  M.,  "Augmentation  of  In-Tube  Condensation  of  R-113",  Ph.D. 
dissertation,  Dept.  of  Mechanical  Engineering,  Iowa  State  University,  Ames, 
Iowa  (1980). 

6.  Akers,  W.  W. ,  Deans,  H.  A.,  and  Crosser,  0.  K. ,  "Condensing  Heat 
Transfer  Within  Horizontal  Tubes",  Chemical  Engineering  Progress  Symposium 
Series  55  (1959):  171-176. 

7.  Boyko,  L.  D.  and  Kruzhilin,  G.  M. ,  "Heat  Transfer  and  Hydraulic 
Resistance  During  Condensation  of  Steam  in  a  Horizontal  Tube  and  in  a 
Bundle  of  Tubes",  International  Journal  of  Heat  and  Mass  Transfer  10 
(1967):  361-373. 

8.  Rosson,  H.  F.  and  Myers,  J.  A.,  "Point  Values  of  Condensing  Film 
Coefficients  Inside  a  Horizontal  Pipe",  Chemical  Engineering  Progess 
Symposium  Series  59  (1963):  190-199. 

9.  Travis,  D.  P.,  Rohsenow,  W.  M. ,  and  Baron,  A.  B. ,  "Forced-Convenction 
Condensation  Inside  Tubes:  A  Heat  Transfer  Equation  for  Condenser  Design", 
ASHRAE  Transactions  79,  Part  1  (1972):  157-165. 

10.  Soliman,  M. ,  Schuster,  J.  R. ,  and  Berenson,  P.  J.,  "A  General  Heat 
Transfer  Correlation  for  Annular  Flow  Condensation,"  Journal  of  Heat 
Transfer  90  (1968):  267-276. 

11.  Azer,  N.  Z.,  Abis,  L.  V.,  and  Swearingen,  T.  B.,  "Local  Heat  Transfer 
Coefficients  During  Forced  Convection  Condensation  Inside  Horizontal  Tubes," 
ASHRAE  Transactions  77,  Part  1,  (1971):  182-201. 

12.  Cavallini,  A.  and  Zecchin,  R. ,  "High  Velocity  Condensation  of  Organic 
Refrigerants  Inside  Tubes",  Reprint  from  the  Proceedings  of  the  13th 
International  Congress  of  Refrigeration,  Washington,  D.C.,  2,  Belgium: 


147 
International  Institute  of  Refrigeration,  (1971). 

13.  Shah,  M.  M. ,  "A  General  Correlation  for  Heat  Transfer  During  Film 
Condensation  Inside  Pipes",  International  Journal  of  Heat  and  Mass  Transfer 
22  (1979):  547-556. 

14.  Lockhart,-  R.  W.  and  Martinelli,  R.  C,  "Proposed  Correlation  of  Data 
for  Isothermal  Two-Phase,  Two-Component  Flow  in  Pipes",  Chemical  Engineering 
Progress  45,  (1)  (1979):  39-48. 

15.  Dukler,  A.  E. ,  Wicks,  III,  M. ,  and  Cleaveland,  R.  G. ,  "Frictional 
Pressure  Drop  in  Two-Phase  Flow:   A  Comparison  of  Existing  Correlations 
for  Pressure  Loss  and  Holdup",  AIChE  Journal  10,  January  (1964):  38-43. 

16.  Dukler,  A.  E. ,  Wicks  III,  M. ,  and  Cleaveland,  R.  G.,  "Frictional 
Pressure  Drop  in  Two-Phase  Flow:   An  Approach  Through  Similarity  Analysis", 
AIChE  Journal  10,  January  (1964):  44-51. 

17.  Hughmark,  G.  A.,  "Pressure  Drop  in  Horizontal  and  Vertical  C0- 
Current  Gas-Liquid  Flow",  I  &  EC  Fundamentals  2,  (4)  (1963):  315-321. 

18.  Collier,  J.,  Connective  Boiling  and  Condensation,  London:  McGraw 
Hill,  (1981). 

19.  Baroczy,  C.  J.,  "Correlation  of  Liquid  Fraction  in  Two-Phase  Flow 
with  Application  to  Liquid  Metals",  North  American  Avaiation  Report 
SR-8171,  El  Segundo,  CA. ,  (1963). 

20.  Thorn,  J.  R.  S. ,  "Predictions  of  Pressure  Drop  During  Forced  Circulation 
Boiling  Water",  International  Journal  of  Heat  and  Mass  Transfer  7  (1964): 
709-724. 

21.  Zivi,  S.  M. ,  "Estimation  of  Steady-State  Steam  Void  Fraction  by  Means 
of  Minimum  Entropy  Production",  Journal  of  Heat  Transfer  86  (1964): 
247-252. 

22.  Turner,  J.  M.  and  Wallis,  G.  B. ,  "The  Separate-Cylinder  Model  of  Two- 
Phase  Flow",  Thayer's  School  of  Engineering,  Report  No.  NYO-3114-6, 
Darthmouth  College,  (1965). 

23.  Hughmark,  G.  A.,  "Holdup  in  Gas-Liquid  Flow",  Chemical  Engineering 
Progress  58,  April  (1962),  62-65. 

24.  Watkinson,  A.  P.,  Miletti,  D.  L. ,  and  Kubanek,  G.  R. ,  "Turbulent  Heat 
Transfer  and  Pressure  Drop  of  Internally  Finned  Tubes",  AIChE  Symposium 
Ser.,  Vol.  69,  No.  131,  (1973):  94-103. 

25.  Watkinson,  A.  P.,  Miletti,  D.  L. ,  and  Kubanek,  G.  R. ,  "Heat  Transfer 
and  Pressure  Drop  of  Internally  Finned  Tubes  in  Turbulent  Air  Flow," 
ASHRAE  Transactions,  Vol.  81,  Part  1,  (1975),  330-349. 

26.  Watkinson,  A.  P.,  Miletti,  D.  L. ,  and  Kubanek,  G.  R. ,  "Heat  Transfer 
and  Pressure  Drop  of  Internally  Finned  Tubes  in  Laminar  Oil  Flow",  ASME 


148 

Paper  75-HT-41,  Presented  at  the  AIChE-ASME  Heat  Transfer  Conference, 
San  Francisco,  California,  August  (1975). 

27.  Carnavos,  T.  C. ,  "Cooling  Air  in  turbulent  Flow  with  Internally 
Finned  Tubes",  Presented  at  the  17th  Natl.  Heat  Transfer  Conf . ,  Salt  Lake 
City,  Utah,  AIChE  Paper  4,  Aug.  14-17  (1977),  also  Heat  Transfer  Engineering, 
Vol.  1,  No.  2,  Oct. -Dec.  (1979):  41-46. 

28.  Carnavos,  T.  C. ,  "Cooling  Air  in  Turbulent  Flow  Multi-Passage 
Internally  Finned  Tubes",  ASME  Paper  78-WA/HT-52,  Presented  at  the  Winter 
Annual  Meeting,  San  Francisco ,  CA. ,  Dec.  10-15  (1978). 

29.  Carnavos,  T.  C,  "Heat  Transfer  Performance  of  Internally  Finned  Tubes 
in  Turbulent  Flow",  Heat  Transfer  Engineering,  Vol.  1. ,  No.  4,  April-June 
(1980):  32-37. 

30.  Hu,  M.  H.  and  Chang,  Y.  P.,  "Optimization  of  Finned  Tubes  for  Heat 
Transfer  in  Laminar  Flow",  Transactions  ASME,  Journal  of  Heat  Transfer, 
Vol.  95,  Series  C,  (1973):  332-338. 

31.  Nandakumar,  K.  and  Masliyah,  J.  H. ,  "Fully  Developed  Viscous  Flow  in 
Internally  Finned  Tubes",  The  Chemical  Engg .  J..,  Vol.  10,  (1975):  113-120. 

32.  Masliyah,  J.  H.  and  Nandakumar,  K. ,  "Heat  Transfer  in  Internally 
Finned  Tubes",  Journal  of  Heat  Transfer,  Vol.  98,  (1976):  257-261. 

33.  Soliman,  H.  M.  and  Feingold,  A.,  "Analysis  of  Fully  Developed  Laminar 
Flow  in  Longitudinal  Internally  Finned  Tubes",  The  Chemical  Engineering 
Journal,  Vol.  14,  (1977):  119-128. 

34.  Soliman,  H.  M.  and  Feingold,  A.,  "Analysis  of  Heat  Transfer  in 
Internally  Finned  Tubes  Under  Laminar  Flow  Conditions",  Heat  Transfer 
(1978),  Vol.  2,  Hemisphere  Publishing  Corporation,  (1978):  571-576. 

35.  Reisbig,  R.  L.,  "Condensing  Heat  Transfer  Augmentation  Inside 
Splined  Tubes" ,  Paper  74-Ht-7 ,  Presented  at  the  AIAA/ASME  Therinophysics 
and  Heat  Transfer  Conference,  Boston,  July  (1974). 

36.  Kroger,  D.  G.,  "Laminar  Condensation  Heat  Transfer  Inside  Inclined 
Tubes",  Paper  No.  AIChE-29,  Presented  at  the  16th  National  Heat  Transfer 
Conference,  St.  Louis,  MO.,  Aug.  8-11,  (1976). 

37.  Vrable,  D.  L. ,  "Condensation  Heat  Transfer  Inside  Horizontal  Tubes 
With  Internal  Axial  Fins",  Ph.D.  dissertation,  University  of  Michigan 
(1974). 

38.  Vrable,  D.  L. ,  Yang,  Wen-Jei,  and  Clark,  John  A.,  "Condensation  of 
Refrigerant-12  Inside  Horizontal  Tubes  With  Internal  Axial  Fins",  In 
Heat  Transfer  (1974),  Fifth  International  Heat  Transfer  Conference  3, 
pp.  250-254.,  Tokyo:   Japan  Society  of  Mechanical  Engineers,  (1974). 

39.  Yang,  Wen-Jei,  "Equivalent  f  and  j  Factors  for  Condensation  Inside 
Horizontal  Smooth  and  Finned  Tubes",  Letters  in  Heat  and  Mass  Transfer 
1  (1974):  127-130. 


149 


40.  Cavallini,  Alberto,  and  Zecchin,  Reberto ,  "A  Dimenslonless  Correlation 
for  Heat  Transfer  in  Forced  Convection  Condensation",  pp.  309-313,  In  Heat 
Transfer  (1974),  Proceedings  of  the  Fifth  International  Heat  Transfer 
Conference  3.   Tokyo:  Japan  Society  of  Mechanical  Engineers,  (1974). 

41.  Royal,  J.  H.  and  Bergles,  A.  E. ,  "Experimental  Study  of  Horizontal 
In-Tube  Condensation",  Presented  at  the  ASHRAE  Dallas  Meeting,  Feb.  1-5, 
ASHRAE  Transactions  82,  Part  II  (1976):   919-931. 

42.  Royal,  J.  H.  and  Bergles,  A.  E.,  "Augmentation  of  Horizontal  In-Tube 
Condensation  by  Means  of  Twisted-Tape  Inserts  and  Internally  Finned  Tubes", 
Journal  of  Heat  Transfer  100,  (1)  (1978):  17-24. 

43.  Royal,  J.  H.  and  Bergles,  A.  E. ,  "Pressure  Drop  Performance  Evaluation 
of  Augmented  In-Tube  Condensation",  In  Heat  Transfer  (1978),  Sixth 
International  Heat  Transfer  Conference,  Vol.  3:  459-464.   Toronto: 
Hemisphere  Publishing  Corporation,  (1978). 

44.  Luu,  M.  and  Bergles,  A.  E. ,  "Experimental  Study  of  Augmentation  of 
In-Tube  Condensation  of  R-113",  ASHRAE  Transactions  Vol.  85,  Pt.  2,  (1979): 
132-145. 

45.  Luu,  M.  and  Bergles,  A.  E. ,  "Enhancement  of  Horizontal  In-Tube 
Condensation  of  Ref rigerant-113",  Paper  Presented  at  the  Session  on 
Condensation  Heat  Transfer,  National  Heat  Transfer  Conference,  ASME  (1979). 

46.  Luu,  M.  and  Bergles,  A.  E.,  "Enhancement  of  Horizontal  In-Tube 
Condensation  of  Ref rigerant-113" ,  ASHRAE  Transactions,  86,  Part  (1) 
(1980),  pp.  293-312. 

47.  Said,  S.  A.  and  Azer,  N.  Z. ,  "Heat  Transfer  and  Pressure  Drop  During 
Condensation  Inside  Horizontal  Finned  Tubes",  ASHRAE  Transactions  (1983) 
Vol.  89,  Pt.  (1). 

48.  ASHRAE  Handbook  of  Fundamentals,  New  York:  ASHRAE  (1981). 

49.  Silver,  R.  S.  and  Wallis ,  G.  B.,  "A  Simple  Theory  for  Longitudinal 
Pressure  Drop  in  the  Presence  of  Lateral  Condensation",  Proceedings  of 
the  Institute  of  Mechanical  Engineers  180,  (1),  (1965-66):  36-42. 

50.  Butterworth,  D. ,  "A  Comparison  of  Some  Void-Fraction  Relationships 
for  Co-Current  Gas-Liquid  Flow",  International  Journal  of  Multiphase 
Flow  1  (1975):  845-850. 

51.  Bergles,  A.  E. ,  Blumenkrantz ,  A.  R. ,  and  Taborek,  J.,  "Performance 
Criteria  for  Enhanced  Heat  Transfer  Surfaces",  Presented  at  13th  National 
Transfer  Conference,  AIChE-ASME,  Denver,  CO.,  Aug.  6-9,  (1972). 

52.  Bergles,  A.  E. ,  Blumenkrantz,  A.  R. ,  and  Taborek,  J.,  "Performance 
Evaluation  Criteria  for  Enhanced  Heat  Transfer  Surfaces",  In  Heat  Transfer 
(1974),  5th  International  Heat  Transfer  Conference,  Vol.  2,  pp.  234-238. 
Tokyo:   Japan  Society  of  Mechanical  Engineers,  (1974). 


150 

53.  Azer,  N.  Z.,  Fan,  L.  T. ,  and  Lin,  S.  T. ,  "Augmentation  of  Condensation 
Heat  Transfer  With  In-Line  Static  Mixers",  Proceedings  of  1976  HTFMI 
(1976):  512-526. 

54.  Lin,  S.  T. ,  "Augmentation  of  Two-Phase  Heat  Transfer  With  In-Line 
Static  Mixers",  Ph.D.  dissertation,  Dept.  of  Mechanical  Engineering, 
Kansas  State  University,  Manhattan,  KS.,  (1979). 

55.  Azer,  N.  Z.  and  Shivakumar,  V.,  "Enhancement  of  Saturated  Boiling  Heat 
Transfer  by  Internally  Finned  Tubes",  ASHRAE  Transactions  (1984),  V.  90, 
Pt.  1. 

56.  Kline,  S.  J.  and  McClintock,  F.  A.,  "Describing  the  Uncertainties  in 
Single-Sample  Experiments",  Mechanical  Engineering  75,  January  (1953):  3-8. 


151 
NOMENCLATURE 


Symbol 

2 

A        Flow  cross-sectional  area,  m 

A1       Constant  in  Eq. (5-26)  and  Table  (5.1) 

A        Actual  heat  transfer  area,  m  /m 

2 

A        Actual  free  flow  area,  m 

2 

A        Open  core  free  flow  area,  m 

Afn      Nominal  flow  area  based  on  inside  tube  diameter  as  if  fins 
were  not  present,  m 

A^       Nominal  heat  transfer  area  based  on  inside  tube  diameter  as 
if  fins  were  not  present,  m^ /m 

b        Fin  height,  m 

C        Constant  in  Eqs.(4-2)  and  (5-1) 

i 
C        Constant  in  Eq. (4-3) 

Cp  Specific  heat  of  liquid,  W.hr/Kg,°C 

Cpwa  Specific  heat  of  cooling  water,  W.hr/Kg,°C 

De  Diameter  of  internally  finned  tube  if  fins  melted  down,  m 

Dh  Hydraulic  diameter  (D  =  4Af  /A  )  ,  m 

D.  Inside  tube  diameter,  m 

D_  Outside  tube  diameter,  m 

Pressure  gradient,  bar/m 

Frictional  pressure  gradient,  bar/m 


®. 


jg-  Frictional  pressure  gradient  assuming  that  liquid  alone  is 

h  flowing  in  the  pipe,  Eq. (5-10) ,  bar/m 

ig£l  Momentum  pressure  gradient,  bar/m 

Igj  Frictional  pressure  gradient  assuming  that  gas  (or  .vapor)  alone 

Jv  is  flowing  in  the  pipe,  Eq.(5-10),  bar/m 

f  Friction  factor 
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S  Q  Friction  factor  in  the  presence  of  condensation,  Eq. (5-23) 

f jc  Friction  factor  developed  for  finned  tubes  by  Said  and  Azer, 
Eq.(5-43) 

f  Friction  factor  defined  in  Eqs.(5-16)  and  (5-17) 

F^  Parameter  defined  in  Eqs.(5-35)  and  (5-41) 

F  Parameter  defined  in  Eqs.(5-36)  and  (5-42) 

F  Parameter  defined  in  Eq.(5-37) 

F  Paramter  defined  in  Eq.(5-44) 

g  Acceleration  due  to  gravity,  g  =  9.81  m/sec 

G  Mass  flux,  Kg/hr.m2 

Ge  Adjusted  mass  flux,  Kg/hr.m2 

G{  Liquid  mass   flux    (G{    =   G(l-x)),    Kg/hr.m2 

Gv  Vapor  mass  flux  (Gv  =  Gx) ,  Kg/hr.m2 

h  Heat  transfer  coefficient,  W/m  °C 

h  Average  heat  transfer  coefficient,  W/m  °C 

h  Local  heat  transfer  coefficient,  W/m  °C 

H  Fin  pitch  (length  per  turn),  m/360°,  or  enthalpy,  Joule/Kg 

H  Enthalpy  of  saturated  liquid,  Joule/Kg 

H  Latent  heat  of  vaporization,  Joule/Kg 

k^  Thermal  conductivity  of  saturated  liquid,  W/m  °C 

K  Thermal  conductivity  of  condenser  tube,  W/m  C 

L  Length  of  condenser,  m 

™wa  Mass  flow  rate  of  the  cooling  water,  Kg/hr 

M;  Saturated  liquid  flow  rate  of  R-ll  (M,  =  M  (1-x) ) ,  Kg/hr 

M  Mass  flow  rate  of  R-ll,  Kg/hr 


Saturated  vapor  flow  rate  of  R-ll  (M  =  M  x)  ,  Kg/hr 


M 

V 

n        Constant  in  Eg. (4-2)  and  (5-1),  or  number  of  internal  fins 
n        Constant  in  Eq.(4-3) 
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P        System  pressure,  bar 

P        Critical  pressure,  bar 
cr 

AP       Pressure  drop,  bar 

P        Reduced  pressure  (P   =  P/P   ) 
re  re      cr' 

Pr  Prandtl  number  of  saturated  liquid 

I  Pumping  power,  W 

q  Constant  in  Eq. (5-26)  and  Table  (5.1) 

Q  Heat  transfer  rate  to  the  coolant,  W 

Qf  Heat  transfer  rate  from  the  condensing  fluid,  W 

r  Constant  in  Eq.(5-26)  and  Table  (5.1) 

R        Pressure  drop  index  for  the  performance  evaluation  Eq. (6-2) 
AP 

R,       Reduction  condenser  size  index  for  the  performance  evaluation, 
Eq.(6-1) 

Re       Reynolds  number 

Re       Reynolds  number  based  on  adjusted  mass  velocity  (Re  =  G  D./u„) 
e  q  e  i  <. 

Re?      Reynolds  number  based  on  liquid  or  liquid  phase  alone  flowing 
in  the  tube  (Re?  =  G(l-x)  D./u,  ) 

ReL      Reynolds  number  based  on  the  liquid  flowing  alone  (Re  =  GD./y.) 

S1       Constant  in  Eq. (5-26)  and  Table  (5.1) 

t        Fin  thickness,  m 

o 
T        Temperature ,   C 

—  o 

T        Average  saturation  temperature  of  R-ll,   C 

o 
T   .     Water  inlet  temperature,   C 
wai  r 

o 

T        Water  outlet  temperature,   C 
wao 

T        Average  wall  temperature,   C 
wo 

o 
AT       Temperature  difference,   C 

w  Uncertainity 

W  Channel  width  between  internal  fins 

x  Dryness  fraction  (ratio  of  vapor  mass  to  total  mass) 

x  Mean  dryness  fraction 
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X        Lockhart-Martinelli  parameter  defined  in  Eq. (5-11) 
Az      Length  of  the  subsection,  m 
Z        Axial  distance,  m 


Greek  Letters 

a        Spiral  fin  tube  helix  angle  (angle  between  fin  and  tube  axis) ,  degrees 

a(X)     Ratio  of  two-phase  friction  factor  to  single-phase  friction  factor 
at  two-phase  Reynolds  number,  Eq.(5-8) 

B       Ratio  of  two-phase  density  to  no-slip  density,  Eq.(5-19) 

5       Parameter  defined  in  Eq. (5-24) 

^       Ratio  of  liquid  volumetric  flow  to  the  total  volumetric  flow 
rate,  Eq. (5-22) 

U       Dynamic  viscosity,  Kg/hr.m 

U       Dynamic  viscosity  of  saturated  liquid,  Kg/hr.m 

u „_      Dynamic  viscosity  of  two-phase  homogeneous  mixture,  Eq. (5-21) 

Kg/hr.m 

u       Dynamic  viscosity  of  saturated  vapor,  Kg/hr.m 

2 
v       Kinematic  viscosity  of  saturated  liquid,  m  /hr 

0       Density,  Kg/m 

P       Desity  of  saturated  liquid,  Kg/m 

P       Density  of  two-phase  homogeneous  mixture,  defined  in  Eq. (5-5) , 
Kg/m3 

Pjjj.      Density  of  two-phase  homogeneous  mixture,  defined  in  Eq.  (5-20)  , 
Kg/m 

P        Density  of  saturated  vapor  phase  Kg/m 

*       Lockhart-Martinelli  parameter,  Eq. (5-12) 

ty  Void  fraction 


Subscripts 

1-4      Subsection  1-4 

aug      Augmented  surface 
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cal 

Calculated 

exp 

Experimental 

in 

Inlet 

NS 

No- Slip 

out 

outlet 

s 

Saturated 

SB 

Smooth  surface 

tt 

Turbulent-turbulent 

z 

Local,  or  axial 
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APPENDIX  A 

ADDITIONAL  INFORMATION  ON  THE  INSTRUMENTATION  AND  COMPONENTS 
USED  IN  THIS  STUDY 


1.    R-U   FLOW  CIRCUIT 


Components : 


1.  Refrigerant-11  Liquid  Circulating  Gear  Pump: 

Sherwood  Alear  Siegler  Company 

Bronze  Rocary  Gear  Pump 

Model:   S  and  V  series 

R.P.M. :   1725 

Pipe  Size:   1/4" 

Shaft  Diameter:   1/2" 

H.P.:   1/3 

Dripless  Mechanical  Shaft  Seal,  Self  Lubricated. 

2.  Refrigerant-11  Liquid  Circulating  Pump  Motor: 

Dayton  -  Electric  A.C.  Motor 
Model  No.:   5K991 
R.P.M.  :   1725 
H.P.:  1/2 
HZ:   60 

3.  Refrigerant-11  Filter: 

Sparlan  -  Catchall  Refrigerant  Filter 
Type:   C-304 

4.  Vapor  Generator  Heating  Elements  (Two): 

Chromalox  Immersion  Type  Heating  Elements 

Model:   AH2745 

Capacity:  240  V,  4.5  KW  each 

5.  Refrigerant-11  Liquid  Receiver: 

Midland-Ross  Refrigeranc  Type  Circular  Tank. 
Serial  No. :   2193 
Size:   3.5  gallons 

Working  Pressure:   Maximum  Allowable  Working  Pressure 
400  psi  at  650°F 
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6.  Refrigerant-11  Valves: 

Diaphram  Packless  Line  Valves 
Superior  Brand,   Solder  Co  Solder  Type 

A.  Model  No.:   214-4S  (1/4") 

B.  Model  No.:   216-10S  (5/8") 

7.  Refrigerant- 11  Vapor-Liquid  Separator: 

Penberty  Co.  V-L  Separator 
Serial  No.:   X- 503070003000 

8.  Refrigerant-11  After-Condenser: 

Dunham-Bush  Bundle  Type  Condenser 
Model  No.:   C1C-200-66-L 

9.  Refrigerant-11  Tube  Connectors: 

Standard  Copper  Tube 
Sweat  Fitting  Type 

10.  Vibration  Eliminators: 

Anacond  Vibration  Eliminators  Supplied  by  RECO 
Specification:  "Has  Fatigue-Resistance  Corrugated  in 

Bronze  Seamless  Tubing  Core  with  Bronze 
Braid  Covering.   Standard  Copper  Tube 
Fittings  are  Welded  on  Both  Ends." 

11.  Refrigerant-11  Vapor  Generator,  Locally  Constructed: 

Material: 

Steel  Tubing:   17.78  cm  O.D.  (16.51  cm  I.D.), 
length  78.74  cm 

Two  Circular  Steel  Plates  =   24.13  cm  Diameter, 

0.95  cm  thick. 

The  two  circular  steel  plates  were  welded  to  close  the 
two  ends  of  the  horizontal  steel  tubing  creating  a  tank.   Holes 
were  drilled  and  tapped  in  each  steel  plate  to  accommodate  the 
threaded  heater  elements  (one  element  in  each  plate,  and 
positioned  axially  opposite  to  each  other  toward  the  bottom  of 
of  the  tank) .   Holes  were  also  drilled  and  tapped  in  a  vertical 
plane  on  one  plate  Co  allow  the  connection  of  a  liquid  level 
indicator  for  the  purpose  of  knowing  if  heaters  were  covered 
completely  by  Che  liquid.   Two  holes  were  drilled  and  capped  on 
one  end  of  Che  cank  Co  allow  liquid  R-ll  Co  flow  inco  Che  cank 
chrough  the  middle  of  the  plate  and  vapor  to  exit  at  a  level 
above  the  center. 
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12.  Superheater,  Locally  Constructed: 

Material: 

Copper  Tubing:   1.59  cm  O.D.  (1.27  cm  I.D.), 

Length  101.6  cm 
Heating  Element:   Ribbon  type  chromel  of  0.204  SJ/ 

30.48  cm 
Teflon  Tape:   Saunder  type  S-17 
Epoxy:   Armstrong  A-68  and  B-68  Types 

The  teflon  tape  was  wrapped  around  the  copper  tube 
with  one  thickness.   The  heating  element  was  wound  uniformly 
around  the  tape  with  0.65  cm  distance  pitch.   The  epoxy  was 
applied  to  secure  the  heating  element. 

13.  Thermocouples: 

Copper-Cons tantan  thermocouples  of  type  B&S  24  gage. 


Instrumentation 


1.  Refrigerant- 11  Liquid  Level  Gauge: 

Brooks  Rotameter  View  Meter 

Type:   6-1355-VB 

Serial  No.:   6507-36340/4 

2.  Refrigerant- 11  Flow  Meters: 

A.  Fischer-Porter  Variable  Area  Type  Flow  Meter 

Range:   0^-0.35  GFM  Liquid 
Model:   10A3565S 
Serial  No.:  7207A4733A2 
Tube  No.:  FP-1/2-27-G- 10/55 

B.  Fischer-Porter  Variable  Area  Type  Flow  Meter 

Range:   0^0.5  GFM  Liquid 
Model:   10A3565S 
Serial  No.:  7207A4733A1 
Tube  No.:  FP-1/2-17-G- 10/55 

3.  Refrigerant- 11  Pressure  Gauge: 

Heise  Pressure  Gauge  of  Type  H28832 
Range:   0^200  psig 
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4.  Pressure  Transducer: 

Pace  Wiancko  Division  of  Whittaker  Corporation 
Model:   KP15  Pressure  Transducer 
Serial  NO.:   150330 

5.  Transducer  Indicator: 

Pace  Wiancko  Division  of  Whittaker  Corporacion 
Model:  CD25 
Serial  No.:   23449 

6.  Digital  Multimeter 

Model:   168  Autoranging  DMM 
Keithley  Instruments,  Inc. 

7.  Voltage  Regulator 

Superior  Electric  Co. 

Powerstat  Variable  Autotransformer 

Input:   240  V,  60  HZ 

Output:   0-280  V,  28A,  7.8  KW 

8.  A.C.  Ampere  Meter: 

Daystrom,  Incorporated  Weston  Instruments  Div. 
Weston  Instruments,  Inc. 
New  York,  New  Jersey 
Model:   433  No.  164330 

9.  A.C.  Volt  Meter: 

Daystrom,  Incorporated  Weston  Instruments  Div. 
Weston  Instruments,  Inc. 
New  York,  New  Jersey. 
Model:   433  No.  146652 

10.  Data  Acquisition  System: 

Esterline  Angus  an  Esterline  Company 

Model:   PD-2064 

Type:   Key  Programmable 

The  system  can  gather  analog  and  digital  data  from  up 
to  64  channels  under  the  control  of  tiny  microprocessor.   The 
system  outputs  the  measured  values  in  engineering  or  scientific 
units  through  various  output  devices.   The  solid-state  integrated 
circuit  microprocessor  is  combined  with  RAMs  (random  access  memory 
devices) ,  ROMs  (read-only  memory  devices) ,  and  PROMs  (program- 
mable ROMs)  to  provide  a  keyboard-programmable  system  that 
permits  the  instrument  to  scan,  measure,  collect,  identify, 
and  record  both  analog  and  digital  input  signals. 
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Accuracy: 

.With  Ambient  Temperature  ac  7  7°F  +  9°F 
+  0.01%  of  reading,  +  0.015%  full  scale,  +  1  count  on  4000  mV  range; 
+  0.01%  of  reading,  +0.03%  full  scale,  +  1  count  on  400  mV  range; 
+  0.01%  of  reading,  +  0.04%  full  scale,  +  1  count  on  40  mV  range. 

Over  Full  Operation  Ambient  Temperature  Range  of  32°F  to  122°F 
+0.5  uV  per  °C,  +  0.01%  of  reading,  +  .04%  full  scale, 
+  1  counc  on  all  ranges . 

11.  Manoraecer: 

Meriam  Instrument  Co. 

Type:   H 

Model:   30EC10 

Serial  No. :   B23131 

Range:   40" 

Manometer  with  mercury  as  the  indicating  fluid. 

12.  Differential  Pressure  Cell: 

Foxboro  Type  13A 

Range:   -0.034*^0.184  bar  (-.5^2.67  psia) 

D/P  cell  connections,  locally  constructed:   Copper 

tubing  of  0.64  cm  O.D.  was  used  to  connect  the 

pressure  tap  and  the  D/P  cell. 
Calibration:   The  differential  pressure  cell  was  calibrated 

according  to  the  manufacturer's  recommended  calibration 

procedure  before  it  was  connecced  to  the  test  condenser. 

A  linear  least  squares  regression  correlation  was  used 

to  obtain  the  calibration  curve  given  by, 

Y  =  0.00021074  x  -  0.031749  (A-l) 
where : 

Y  =  pressure  drop  in  bars 

X  =  D/P  cell  output  in  mm  of  mercury. 

A  total  of  about  100  points  were  used  to  calibrate 
the  pressure  cell. 
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13.   Vaccum  Pump 

Matheson  Scientific 

Division  of  Will  Ross,  Inc. 

Serial  No. :   1173 

Power:   115  V,  60  HZ 

Connections:   3  conductor  power  cord  with  2-prong  adapter. 

Inlet  and  outlet  connector  to  3/8"  I.D.  hose. 
Function:   Portable  A.C.  powered  source  of  vaccum  (to  686mm/ 
27"  Hg)  or  pressure  (to  1,7  kg/cm  ,  25  psig) . 


R-ll 
[48] 


All  theraophysical  and  transport  properties  of 
were  obtained  from  ASHRAE  Handbook,  of  Fundamentals 


2.   HATER  FLOW  CIRCUIT 


1.   Cooling  Water  Flow  Meters: 

Brooks  Rotameter 

Type:   1110-09H3A1B 

Serial  No.:   7201-74650/1 

Tube  No.:   R-9M-25-1  BR-3/4-14G10 

Range:   0^3  GPM 


2.   Cooling  Water  Pump: 


A.O.  Smith  Co.  Pump 
Model  No.:   C48L2DA11A4 
Serial  No. :   J69 
H.P.:   1 
R.P.M. :   3450 
HZ:   60 
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SAMPLE  OF  DATA  REDUCTION  AND  CALCULATION  PROCEDURE  OF 
HEAT  TRANSFER  COEFFICIENTS 

The  raw  data  of  all  experimental  runs  were  fed  into  a  digital 
computer  prorgram  which  was  written  to  reduce  the  data  into  useful 
forms.  The  computer  program  listing  is  given  in  Appendix  C.  The  reduced 
data  are  given  in  Appendix  D.  The  experimental  runs  in  the  reduced  data  are 
coded  in  three  digit  numbers.  The  first  digit  represents  the  tube  number 
and  the  remaining  two  digits  represent  the  run  number.  A  sample  of  the 
calculation  procedure  for  run  number  146  is  given  below.  Since  all  the 
experimental  measurements  like  temperatures,  pressures,  flow  rates  etc. 
were  in  British  units,  the  calculations  for  data  reduction  were  done  in 
British  units  and  quantities  of  interest  were  converted  into  SI  system 
of  units.  During  the  calculation  procedure  reference  is  made  to  the 
thermocouple  stations  and  pressure  tap  locations  in  Fig  B.l 


Recorded  Experimental  Data  for  Run  number  146 

Reference  (ambient)  temperature  =  74°F 

Atmospheric  pressure  =  14.10  psia 

Tube  I.D.,  D±  =  0.545  in. 

Tube  O.D.,  Dq  =  0.625  in. 

Tube  length  of  each  subsection,  flz  =  27.0  in. 

Number  of  subsections  in  the  tube  =  4 

Tube  thermal  conductivity,  K  =  220 

Test  fluid  flow  rate,  GPMF  =  0.105 

Water  flow  rate  through  the  test  condenser,  GPMW  -  1.58 


Btu 


hr  ft2  °F 


gal 


gal 
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•j 

3 


- 

a 


Test  section  inlet  pressure,  P. 

In 

Total  pressure  drop,  AP 


24.10  psia 
0.066  psi 
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Coolant  temperatures 


Coolant 
Coolant 
Coolant 
Coolant 
Coolant 
Coolant 
Coolant 
Coolant 


inlet  temperature  at  subsection  4,  T. 


25 


outlet  temperature  at  subsection  4,  T 


20 


inlet  temperature  at  subsection  3,  T 


19 


outlet  temperature  at  subsection  3,  T 


14 


inlet  temperature  at  subsection  2,  T 


13 


outlet  temperature  at  subsection  2,  T 

8 

inlet  temperature  at  subsection  1,  T 
outlet  temperature  at  subsection  1,  T 


66.56  °F 
68.00  °F 
67.28  °F 
68.90  °F 
69.62  °F 
72.50  °F 
73.14  °F 
73.94  °F 


R-ll  Temperature 

Test  fluid  inlet  temperature,  T 

Test  fluid  outlet  temperature,  1 

26 

Test  fluid  temperature  at  flowmeter,  T, 


27 


110.02  °F 
96.98  °F 
86.72  °F 


Wall  Temperatures  in  °F 
Subsection  1       Subsection  2 


T3  =  96.08 

T,  =  102.74 
4 

T5  =  94.46 

T,  =  96.98 

o 


31.14 


TJ0  =  90.32 
T  =  75.02 
TJ2  =  82.76 


Subsection  3 


T   =  74.84 


T.,  =  90.68 
Id 


T^  =  73.22 


Tlg  =  84.56 


Subsection  4 

T       : 
21 

=   73 

.58 

T 
22 

=   83 

.66 

T 
23 

"   69 

.80 

T 
24 

■   84, 

.92 

Calculation  Procedure 

Sectional  heat  transfer  area,  A  : 


•V  D.  A  z 
l 


=  Tr(0.545)(27) 
144 


=  0.321  ft 
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Cross  sectional  flow  area,  A,  : 

fn 

2 
A.     =      TrDl     =        (0.545)2      =    1.62002   x   10"3   ft2 
4  (4) (144) 

Coolant  flow  rate  into  the  test  condenser,  m   : 

wa 

.   ,  1-58  gallons   6Q  minute  .j£_         Ibm 

wa        minute        hour  gallon  3 

=  789.52  lbm/hr 


Test  fluid  flow  rate,  M  : 
T 

Since  the  flow  meter  was  calibrated  for  R-113,  a  correction  factor 

of  1.036  was  used  for  R-ll  flow  rate. 

L   -  (0.105  x  1.036)  3allons  x  60  ^nute   Q      _|£ 
T  minute         hour  gallon 

x  91.33  ^ 
ft 

=  79.70  lbm/hr 

Test  fluid  mass  flux,  G  : 

G  =  M/A   =  ii^Z° =  0.49197  x  105  lbm/hr. ft2 

1.62002x10 

Heat  transfer  rate  at  each  subsection,  Q. : 


Q,  =  m   Cp    (T   -  T   . ) 
i    wa  rwa   wao   wai 

Q1  =  (789. 52)(1)(73. 94-72. 14)  =  1421.14  Btu/hr 

Q2  =  (789. 52)(1)(72. 50-69. 62)  =  2273.82  Btu/hr 

Q3  =  (789. 52)(1)(68. 90-67. 28)  =  1279.02  Btu/hr 

Q4  =  (789. 52)(1)(68. 00-66. 56)  =  1136.91  Btu/hr 


Total  heat  transfer  rate,  Q  : 

Q  "  Qj  +  Q2  +  Q  +  Q   =  6110.88  Btu/hr 
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Test  fluid  heat  transfer,  Q  : 

Qf  -  Mt[hi(Pi;Ti)  -  Hf5] 

=   79.70    [105.36   -   28.07] 

=  6160.18  Btu/hr 

where 

H^  is  the  specific  enthalpy  of  the  test  fluid  at  station  1,  Btu/lbm 

H   is  the  specific  enthalpy  of  the  saturated  R-ll  at  T-.,    Btu/lbm 
tD  26 

Energy  balance  error  in  percentage,  Error  %  : 

Error  Z   =  100  x  Q  ~  Qf 
Q 

6110.18-6160.18 


=  100  x 


6110.18 
=  -0.81% 

Quality  at  each  station,  x,  : 

l 

The  quality  was  calculated  from  the  energy  balance  for  each  substation 
in  which  the  energy  gained  by  coolant  and  the  energy  lost  by  condensing 
fluid  were  assumed  to  be  equal. 

Qj  =  MT(H1-H2) 

Q2  ■  MT(H2-H3) 

Q3  -  MT(H3-HA) 

Q4  =  MT(VV 

Qj  ~  Q.  are  known  and  proceeding  from  inlet  to  exit  the  above  equations 
were  solved  for  enthalpies  H  at  each  station.  The  results  were 
H  =  87.53  Btu/lbm  H,  =  42.95  Btu/lbm 

H3  =  59.00  Btu/lbm  H5  =  28.69  Btu/lbm 
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Hl  "  "fl  _  105.36-29.03  =  ,  Q 


1    H,  ,  75.02 

fg2 

H3  "  Hf3 
x  =  -~ -  =  59.00-29.03  =  0.339 

fg3         75.02 

_  H4  "  Hf4  _  42.95-29.03  _niafi 
fg4 

fg5 

where 

Hf  and  H   are  the  saturated  liquid  enthalpy  and  the  latent  heat, 
respectively,  corresponding  to  the  saturation  pressure  at  each  station. 

In  the  present  study,  the  magnitude  of  the  pressure  drops  was  low 
and  hence  the  effect  of  pressure  drop  on  saturated  liquid  enthalpy  and 
latent  heat  was  neglected.  It  was  assumed  that 

Hfl=  Hf2=  Hf3=  Hf4=  Hf5 


and 


H   =  H,  ,=  H.  ,-  H.  .-  H.  . 

tgl   fg2   fg3   fg4   fg5 


Heat  Transfer  Coefficients 


Average  outside  wall  temperature,  T   : 

wo 


_       T  +  T  +  T.  +  T. 

Subsection  1   T   ,  -  — — -, =  97.57  °F 

wol  4 


Subsection  2   T  ,  =  T9  +  T10  +  Tll  +  T12  =  82.31  °F 

Wo2     ; 

4 
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Subsection   3        T 


T        +  T        +  T        +  T 
15        x16        ^7         18 


wo3 


80.83     F 


Subsection  4        T 


21         22        x23         24 


wo  4 


=   77.99     F 


Sectional   average  heat   transfer  coefficient,    h. 

-1 

TlB.iz    (T  ,-T      .)        DJ  D 

h.   = 


TTD.Az(T.-T      .)        D.  D 

- si       woi  1      ,    ,     o, 

Q~ 2K     la<V? 


(B-l) 


T   and  T   .  are  the  average  saturation  and  average  outside  wall 
si      woi  & 

temperatures  for  the  corresponding  subsection  respectively. 

As  discussed  earliear  in  Chapter  IV,  the  sectional  heat  transfer 

coefficient  of  section  1  was  calculated  by  replacing  T  .  in  Eq.(B-l) 

si 

by  (Tj  +  T2)/2 
where 

T  is  the  inlet  superheat  temperature  of  R-ll 
Since  the  effect  of  pressure  drop  on  saturation  temperature  was 
neglected 


T.=T,=T„=T.=T.=  101.53  F 
si    si    s2    s3    s4 

For  this  run  T  =  110.02  °F 


hl   = 


IT  (0.545)  (27)  [(101. 53+110. 02)/2   -   97.57]  0.545  0.625, 

(144) (1421.14)  "      (24) (220)         (0.545) 


-,  -1 


=  543.41  Btu/hr.  ft2.°F 


h2  = 


TT(0.545)(27)(101.53-82.31)  _    0.545      0.625," 
(144)(2273.82)  (24)  (220)  ln(0.545:) 


370.49  Btu/hr. ft2. °F 
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TT(0.545)(27)(101.53-80.83)    _        0.545        u  ,0.625 


(1A4) (1279.02) 


(24) (220)        '0.545 


192.97  Btu/hr.ft2.°F 


TT  (0.545)  (27)  (101 .53-77  .99)    _        0.545 


(144) (1136.91) 


(24) (220) 


.     ,0.625,1 
ln(0^5)J 


150.78  Btu/hr.ft2.°F 


Overall  average  heat  transfer  coefficient  (tube),  h  : 

-1 

TTD.  L(T  -  T   )    D,    D 
.       l    s   wo     i  ,  ,  o. 


T   is  the  average  saturation  temperature  of  the  test  condenser 
corresponding  to  the  average  pressure. 

T   is  the  average  wall  temperature  of  the  test  condenser 
wo 

=  (T   ,+  T  ,+  T  ,+  T   . )/4.0 
wol   wo 2   wo 3   wo4 

L        is    the   length  of   four  subsections   =4      z. 
Therfore, 

-1 
TT(0.545)(108)(101.53-84.68)    _        0.545  ,0-625, 

(144)  (6110. 18)  (24)  (220)    Xrn0.545; 


=  283.46  Btu/hr.ft2.°F 
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APPENDIX  C 

COMPUTER  PROGRAM  USED  IN  DATA  REDUCTION 

C     *********»****»******»*********»»*»**»«.»»*»»».*,»*»»»»„»»**,„ 

c 

c               *  DEFINITION    OF    SYMBOLS  > 

C 

C               AFLOW  -TEST    CONDENSES    CEOSS    SECTIONAL    ABEA,FT»*2 

C               AISF  -SECTIONAL    HEAT   THANSFES    AREA,FT**2 

C              CPW  -COOLANT    SPECIFIC    HEAT, BTU/LBH    DBG    F 

C               DPT  -TOTAL    PBESSUEE    DEOP,    PSI 

C               DPTS  -TOTAL    PRESSUSE    DEOP,    BAR 

C  DTD1-DTDU-SECTIONAL    DRIVING    TEMPERATURE    DIFFERENCE , DEG    F 

C               DTDA  -TEST   CONDENSES    DRIVING    TEMPERATURE    DIFFERENCE,     DEG    F 

C               DTDAS  -TEST    CONDENSES    DBIVING    TEMPERATURE    DIFFFEFNCE,    DEG    C 

C              ERROR  -PERCENTAGE    DEVIATION    BETWEEN    THE    HEAT    GAINED    BY    COOLAI"* 
AND    THE    HFAT    LOST    BY    CONDENSING    FLUID 

GPMF  -CONDENSING    FLUID    FLOW    PATE   THROUGH    THE    CONDENSES,     GPH 

GPMW  -COOLANT    FLOW    SATE    THEOUGH    THE    CONDENSES,    GPM 

r               ™E!S»-  ,-C0OLA,'T    F10H    MT2    THBOUGH    THE    CONDENSES    II!    SI    UNITS, KG/HS 

C  HTC1-HTCU-SECTIONAL    AVESAGE    HEAT    TBANSFEE    COEFFICIENTS,  BTU/HR    FT«2    DEG    F 

C               HTCA  -OVEBALL    AVERAGE    HEAT    TRANSFER    COEFFICIENT,     B^U/HB    FT**2    DEG    F 

C               HTCAS  -OVEBALL    AVERAGE    HEAT   TRANSFER    COEFFICIENTS,    W/M**2    DEG    F 

C  HF  -ENTHALPY    OF    SATURATED    LIQUID    AS    A    FUNCTION    OF    AVERAGE    SATUSATED 

<~  TEHPESATUBE,     BTU/LBM 

C  H3  -ENTHALPY    OF    SATURATED    VAPOUR    AS    A    FUNCTION    OF    AVERAGE    SAT 

C  UEATED    TEMPEEATUBE     ,BTU/LBM 

C  BIN  -ENTHALPY    OF    INLET    TEST    FLUID,    BTU/LBH 

C  BO  1-HOI  -EXIT    SECTIONAL    ENTHALPY, BTU/LBM 

C  !C  -THEEHAL    CONDUCTIVITY   OF    THE    SATUSATED    LIQUID,  BTU/HB    FT 

C  no  -INSIDE    DIABETES    OF    THE    CONDENSES    TUBE,    FT 

IOD  -OnTSIDE    DIAMETER    OF    THE    CONDENSES    TUBE,    FT 

C  L  -SECTIONAL    TEST    CONDENSES    LENGTH,    FT 

C  LDEN  -DENSITY    OF    SATUSATED   LIQUID    AS    A    FUNCTION    OF    SATURATED 

C  TEMPEEATUBE,     LBM/F?*»3 

C  LDENFM  -DENSITY    OF    LIQUID    E-ll     AT    FLOW    METEE    TEMPEEATUBE,    LBM/*T*»3 

C  N  -NUMBER    OF    POINTS    WITH    DIFFFENT    BASS    VELOCITIES 

C  PI  -A    CONSTANT 

C  PA  -AVESAGE    PEESSnEE    INSIDE    THE   CONDENSES,    PSIA 

C  Pj.N  -INLET    PBESSUEE    OF    THE    CONDENSING    FLUID,    PSIA 

C  PINS  -INLET    PBESSUEE   OF   THE   CONDENSING    FLUID    IN    SI    UNITS,    BAB 

C  QF  -TOTAL    HEAT    TRANSFEBED    FEOH    THE    TEST    CONDENSES,     BTU/HR 

C  QFS  -TOTAL    HEAT    TSANSFEBED    FBOB    THE    TEST    FLUID,     W 

C  QW  -TOTAL    HEAT    TBANSFEE    RATE    FROM    THE    TEST    CONDENSES,     BTU/HS 

C  QWS  -TOTAL    HEAT    TBANSFEE    RATE    FEOM    THE    TEST    CONDENSER       W 

C  QW  1-QUH  -SECTIONAL    HEAT   TBANSFEE    RATE,    BTU/HS 

C  T1-T27  -TEMPEEATUBE    AT    DIFFERENT    LOCATIONS,     DEG    F 

C  T1S  -INLET   TEMPERATURE    OF    THE    CONDENSING    FLU^D,     DEG    C 

C  TG  -TEST    FLUID    TOTAL    H ASS    VELOCITY,     LBM/HB    FT**2 

C  TGS  -TEST   FLUID    TOTAL    BASS    VELOCITY,     KG/S    MT*»2 

C  THF  -TEST    FLUID    BASS    FLOW    SATE,    LBB/HB 

C  TKW  -COOLANT    MASS    FLOW    EATE,     LBM/HE 

C  TSA  -SATURATION    TEBPEBATUBE    IN    THE    TEST    CONDENSES,     DEG    F 

C  TSAS  -SATUEATION    TEBPEBATUBE    IN    THE    TEST    CONDENSER,     DEG    C 

C  TfcAI-TWAH-AVESAGE    SECTIONAL    WALL    TEMPEEATUBE,     DEG    F 

C  TCA  -OVEBALL    AVESAGE    WALL    TSMPEBATUEE,    DEG    J 

C  TWAS  -OVEBALL    AVERAGE    WALL    TEBPEBATUBE,    DEG    C 

C  TVAP  -FUNCTION    USED    TO    CALCULATE    SATUBATION    TENPERATUR" 

C  HI  -INLET   TEST    FLUID    BASS    QUALITY 

C  Xl-JU  -EXIT    SECTIONAL    TEST   FLUID    MASS    QUALITY 
REAL   IID, IOD, LDEN, LDENFM 
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2  DIMENSION    PIS  (65)     DPT<65)  ,PINS(65)  ,DPTS  (65)  ,  GPHW  (65)  ,GPMWS  (65)  ,QW1 
'  (65)  ,QK2(65)  ,QW3  (65)  ,gn  (65)  ,QW(65)  ,QWS  (65)  .1  1  (65)  ,T    S  (t 

3  DIMENSION    DTD  1  (65)',  DTD2  (65)  ,  DTD3  (65)  ,DTD4  (65)  ,TTDA  (651  ,  DTDAS  (651 

1  GPMF  (65)  ,  EBHOR(65)  ,  HTC1  (65)  ,  HTC2  (65)  ,HTC3  (65)  ,KTC4    65    ! 

\  ""(t5'"HTCAS<fc5),XIl!(t5),Xl(6S),I2(65),X3(t5),X»    65 

4  2r^U?PS  (tS'  'GT'fc5>  -GTS(65)  ,TMW(65)  .POUT  (65)  ,  TCO  (65)  , 

«  DIMENSION    PIN2(5,65),PINS2(5,65)  ,DPT2(5,65)  ,DPTS2(5    65) 

1  GPMW2(5.65).GPMWS2(5,6  5)  ,0. HI  2  (5,  65)  ,QW22  (5,6  5)  , 

2  2W32<5,65),QW42(5,65|  ,0.802(5.65)  ,QWS2  (5, 65|  ,T012  (5 ,65)  , 

3  T1S2(5,65)  ,TSA02(5,65) 

5  DIMENSION   TSAS2<5,65),TWA12(5,65),TWA22(5,65),TWA32(5,65). 

1  T»A42(5,6S),TWA02(5,65)  ,TWAS2(5,65)  ,DTD12(S,65 \  . 

2  DTD22(5.65),DTD32(5,65)  ,  DTD42  (5  ,  65)  ,  DTDJ2  (5,6  5)  , 

3  DTDAS2(5,65)  ,GPMF2  (5,65)  ,  ERROR  2  (  5,  65)  ,HTC12(5,65)  , 
»  HTC22(5,65|  ,HTC32(5.65)  ,HTCU2(5,65)  ,  HTCA2  (  5,  65)  . 

5  HTCAS2(5,65) 

6  DIMENSION    XIN2(S,65),X12(5,65),X22(5,65|,X32(5,65)  ,X42(5,65), 

1  QF2(5,65),QFS2(S,65),GT2(5,65)  ,GTS2  (5  ,65)  ,  TMW2  (5,  65)  , 

2  POOT2(5,65),TC02(5,65).TCOS2(5,65)  ,NCOUNT(10) 

C   ***»»*«    JOEHAT    STATEMENTS 
C 
1  10    FORMAT(4F10.5,3I10) 

8  11    FORMAT  (4F10.  5) 

9  12    FORMAT(  10F8.  2) 
10               112    FORMAT  (1H1) 


11  21    FORMAT (/////,T13,  'TABLE     (D-  1)  :-',/,'-•,  1  IX  ,80  ('-' )     /) 

12  22    F0RMAT(T13, 'RUN', T18, 'INLET    PBESS. ', T3 1 ,' INLET    TEMP,  '.TM 

.,  1'COND.FLUID    MASS    TEL. ' ,166, • COOLANT    TIN ' ,T80, ' COOLABt'rATE'     /) 

13  23    FDRHAT(T14,<#',T18,'PSIA         BAR' ,T3  1 , 'DEG    F    DEG    C  ,T44,  '  L3M/HR-FT2 

1       KG/S-M2       DEG    F      DEGC         GPM  KG/HR  •  ,/,  121 ,  80  (•-'),  /) 

14  25,FORMAT(T13,I3,T17,F5.2,T24,FS.3.T30,?6.1,T37,F5.1.'*44,E11    5    "57 

1F6.2,T66,FS.2,T73,F5.2,T80,F5.3,T87;F6.2)  ' 

15  26    FORMAT (/////, T 10,  'TABLE     (D-2)  :-',/,' -'ax, 97  (■-' )     /l 

16  27     PDHMAr(T10, 'RUN', T27, 'QUALITIES', T53, 'SAT.     TEMP. ',  T66.  •  AVE  AGE    SSL 

1L    TEMPEPATTJHES    IN    DEG    F       DEG    CM  AV-RAGE    HAL 

17  50    FORMAT (T56, 'TUBE') 

18  28    FOFKAT(T11,'#',T17,'IIN',T2S,'I',T31,'II',T38,'III'    T45    '17'    T52 

1'DEG    F',T59,'DEG    C  ,  T6  9,  •  I '  ,  T75  ,  'II '  ,T82  ,  •  III'  ,T89  ,  '  IV  ^95.  'TDBE- 
2,1102, 'TUBE', /, 9X, 97('-<)  ,/) 

19  30    K1RMAT(T10,I3,T15,F6.3,T22,F6.3,T29,F6.3,T36,F6.3,T43,F6.3,T51 

lF6.2,T59,F5.2,T66,F6.2,T73,F6.2,T80,F6.2,T87,F6.2,T94,Fe.2,T10  2, 

20  31    FORMAT  (/////,  T13,  'TABLE     (D-3)  :-',/,'-',  1  IX  ,9  3  ('-' )     /) 

21  32    F0RMAT(T13,'H0N       HEAT    TRANSFER    COEFFS.     IN    STU/HR-SQFT«F '    T58 

1'W/SQM*C         DRI7ING   TEMP.     DIFFERENCES    DEG    F    DEG    C  ,/) 

22  33    F0SMAT(T14,'*<,T21,'I  II  In  17  TUBE' ,T59,  '  TUBE' 
1,T68,'I               II              in            17            TUBE         TUB  E'  ,/,  I  2X  ,Q  3  [ '  -  ' )     /) 

FORMAT(T13,I3,T17,F7.2,T25,F7.2,T33,F7.2,T41,F7.2,T49,F7   2 ''58 
1F6.  1,T66,F5.2,T73,FS.2,T80,F5.2,T87,F5.2,T94,F'5.2,T10  l,P5    21*       ' 

24  36    FORMAT  (/////, T33,  'TABLE     (D-4)  :-',/,'-',  3  U, 66  ('-•)  ,/) 

25  37    FORMAT(T33, 'RON', T49, 'INLET    PRESSURE' , T77, "TOTAL    PRESSURE    DROP', 

26  38    FORMAT  (T34,'*', 147,  'PSIA-.T62.' BAR', T78,'PSIA ',193,' PASCAL',/, 

132X,  66  ('-•)  ,/) 

27  40    F0RMAT(T33,I3,T46,F5.2,T60,F6.3,T74,F8.3,T9I,F8.21 

28  41     FORMAT (/////, T13,'TABEL     (D-5)  :-',/,  12X,  39  ('-•)  ,/) 


23  35 


29 


31 
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"2l^?"'Ea(IoRJ,T''T30',CCOLA'''r  BTa/HE,<T62',TnBE'<"8,.TEST  '"»«  , 

30  u3  raaaw(Ti«,''i',«i,'i  n  m  iv,ts5,-btu/fr 

I    W,T76.-BT0/HE        I  )!',/,12Z,89(«-.)  " 

45  FOBMAT(T13,I3,T18.F7.2,T27,F7.2,T36,F7.2,T45,P7.2.TS4  F8  2  T65 
1F7.2,T75,F8.2,T86.F7.2,T96,F6.2)  '-»,«•,«.  2,165, 

C    *"*****"*******»*»*****»*»**»*«.«»*»»«*.,,.»,»,»,»„»„»»„„„ 

C     CALCULATION  OP  SATURATED  LIQUID  ENTHALPY  IN  BTU/L3M  AS  A  FUNCTION  OP 
C      SATUBATION  TEHP.  DBG  F  (CURVE  FIT  FBOM  ASHP.AE  TABI  ES  )   ?n"CTI°''  °F 

32         HFr  (X)=8.  0732+0.  203944*X+1.495E-05*X*:r+8.29B-0a*X*X*X 

r      ™1«J:^™^2Fm!"t'S*!!D  7AP°°a  ««■»"«  I"  BTU/IBM  AS  A  FUNCTION 


C      OF  SATURATION  TEHP.  DEG  F  (CUBVE  FIT  PROM  ASHBAE  TABLES 


C 


•MMMMttWMMMMntMMMMtMMIMMMtMtnMMI 


********* 


33  HGF(X)  =91.942+0.  1 21*X+9E-06 *X*X- 1 . 92E-07*X*X*X 

C  CALCULATION    OF    SATURATED    LIQUID    DENSITY    LBM/CU    FT    AS    A    FUNCTION    OF 

C  SATURATION    TEMP    DEG    F     (CUBVE    FIT    FROM    ASHBAE    TABLES    )        ™N<-li0B    °F 

34  WEN  (X)=96.  734-0.  02066*X-7.068E-04*X*X+2.  61 1E-06*X*X*X 

C  CALCULATION    OF    SATUBATED    TEMP.    IN    DEG    F    CORRESPONDING    JO 

SATURATION    PB.    IN    PSIA     (CURVE    FIT    FROM    ASHRAE   TABLES    | 

35  TVAP(X)  =5.  4895+6. 3303*X-0.  12726  l*X*X+  I.  243  22F-0  3*X*X*X 

C 

c 

C   ******    5E1D    INPUT    DATA 
C 

36  ISTABT=  1 

37  800    CONTINUE 

38  READ  (5,  10]     PI,IID.IOD,CPW,N,L,K 
C  »RITE(6,10)     PI,IID,IOD,CPW,N,L.K 

39  NCOUNT(ISTART1  =N 
10                         DO     1000    1=1, N 

"'  READ(S,  II)     PIN(I)  ,DPT(I1,GPMH(I)  ,UPMF(I) 

C  WRITE(G,11)     PIH(I),DPT(I)  ,GPMW(I)  ,GPMF(I) 

42  SEAD(5,12)    71  (I),  12,13,14,15,16,17,18,19,110,111,112,113,  "14    T1S 

1T16,T17,T18,T19,T20,T21,T22,T23,T24,T25,T26,T27,HIN 

c 
c 

C   ******    CALCULATION    OF    THE    AREAS 

C 

13  AFLOH=PI*IID*IID/576 

"4  AISF=PI*IID*L/144 

C 

C    ******    CALCULATION    OF    HEAT    GAINED    BY    WATER 

C 

45  TC0(I)=T25 

**  THW(I)=GPSH  (I)  *60*0.  13368*62.3 

*7  Q»1  (I)  =  TSW(I)  *CPW*(T2-T7) 

48  QW2  (I)=TMSJ(I)*CPH*  (T8-T13) 

*9  QH3  (I)  =TMS(T)  *CPW*(T14-T19) 

50  QH4  (I)  =TMW(I)  *CP«*(T20-T2S) 

51  QW(I)»QB1  (I)  +O.W2  (D+QIJ3  (I)  +QW4(I) 


WRITE  <6, 12}T1  (I)  ,  12,13,  T4,  T5,  T6,  T7  ,T8,  T9,T  10  ,T1  1  ,T12  ,  T13  .11  4    T15 
1T16,T17,T18,T19,T20,T21,T22,T23,T24,T25,T26,T27,HIN 
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C    ******    CALCULATION    OP    SATURATION    TEtiPESATUSE 

52  ?OOT(I)=PI!T(I)-DPT(I) 

53  TSA(I)=TVAP(PIS(I)) 
C 

C  ******  CALCULATION  OF  HEAT  TBANSFER  COEFFICIENTS 

54  B=Iin/2H./K*ALOG(IOD/IID) 

55  TSUI  (I)  =  (T3+T4+T5*T6)/4. 

56  TSA2(I)  =(T9+TI0+T1HT12)/4. 

57  T"A3(I)=(T15+T16+T17-fT18)/U. 
5a  TWA4(I}=(T21+T22+T23+T24)/4. 

«  TBA(I)  =  (TWA1  (I)+TWA2(I)+TWA3  (H +TWA4  (I)  ) /n. 

60  DTD1  (I)  =  (TSA<I)+T1(I))*0.5-TWA1(I) 

61  DTD2(I)=TSA(I)-TSA2(I) 

62  DTD3  (I)  =TSA  (I)  -TU  A3  (I) 
13  DTU4  (I)  =TSA  (I) -TWAU  (I) 

64  DTDA(I)  =TSA  (I) -TWA  (I) 

65  C1=PI*IID*L*DTD1(I)/(QW1  (I)*l«4.) 

66  C2=PI*IID*L*DTD2(I1  /(QW2(I)  •  HK».  I 

67  C3  =  PI*IID*L*DTD3(I)/(QW3(I)  *144.| 
f  a  C1=PI*IID*I,*DTD<1<I)/(QW4(I)  *1U4.) 

6  9  CA=4.*PI*IID*L*DTDA(I)/(0.S(H  *144.) 

70  HTC1  (I)=1.0/  (C1-B) 

71  HTC2(I1  =  I.0/(C2-B) 

72  HTC3  (I)=1.0/(C3-B) 

73  HTC4(I1  =I.0/(C4-B) 
7"  HTCA(I)=1.0/(CA-B) 

C 

C    ******    CALCULATION    OF    QDAIITIES 

C 

75  HF=SFF(TSA(I)  ) 

76  HG=HGF(TSA(I)  ) 

77  n»(I)=(HIH-HF)/(HG-HF5 

78  LDENFM=LDEN  (T27) 
T!tF=fcO.0*O.  1336a*GPKF(I)*LJ3BNF11 

8  0  GT  (I)=THF/AFLOW 

8  1  HO  I=HIN-QW1  (I)  /TWF 

82  X1  (I)=(H01-HF)/(HG-HF) 

83  H02=HO  1-QW2(I) /T.1F 

8"  X2  (I)=  (H02-HF)/(HG-HF) 

8  5  H0  3=KO2-QW3(I)/TNF 

86  X3  (I)=(H03-HF)/(HG-HF) 

87  H04=H03-QW4(I) /TKF 

88  M  (I)=(H04-HF)/(HG-HF) 
C 

C    ******    OVEEALL    HEAT    BALANCE 
C 

89  HFO=HFF  (T26) 

90  QF  ( I)  =TWF*  ( HIN-HFO) 

91  EEEOS(I)=  IOO.0*!QW(I)-QF(I))/QW(I) 

C    ******    C0NVE2SI0N    TO    SI    UNITS 
C 

92  PINS(I)=PIN<I)/U1.5 

93  TlS(I)  =  (TI(I)-32.0)/l.a 

94  TCOS  (I)  =  (TCO(I)-32.0)/1.8 

95  GTS  (I)  =GT  (I)  /737.339 

96  GPawS(I)=GPMf!  (I)  *226.662 

97  TSAS  (I)  =(TSA(I)-32.  01/1.8 


177 


?,®  TWAS(I)=(TWA(I)-32.0)/1.8 

99  HTCAS(I)  =  HTCA(I)*5.6783 

100  DTDAS(I)=DTDA  (IJ/1.8 

!01  DPTS(I)  =DPT(I)  *t.894757E+03 

]°l  QWS(I)=QW(I)*0.2928751 

J03  QPS(I)=Q?(i)  *0.292875l 

104  PIH2(ISTABT,I)=PIN(I) 

105  PINS2(ISTART.I)=PIHS(I) 

106  TO  12  (ISTABT,!)  =T1  (I) 

107  T1S2  (ISTABT,!)  =T1S  (I) 

108  GT2(ISTABT,I)=GT(I>- 

109  GTS2  (ISTART,I)=GTS(I) 
1,0  TCn2(ISTART,I)  =TCO(I) 

TCOS2  (ISTABT,!)  =TCOS  (I) 

•2  GPHW2(ISTAPT,I)=GPMS(I) 

113  GPHWS2  (ISTABT, I)  =GPMRS  (I) 

■J*  TSA02  (ISTABT,!)  =TSA  (I) 

TSAS2(ISTART,I)  =TSAS  (I) 

11*  TWA12(ISTABT,I)=TWA1  (I) 

""  TWA22(ISTABT,I)  =TWA2  (I) 

118  TWA32  (ISTABT, I)=TSA3  (I) 

1,9  TWA42(ISTABT,I)=TVA4(I) 

120  TUA02  (ISTABT,!)  =TWA  (I) 

,2'  TWAS2(ISTAET,I)  =TWAS  (I) 

122  HTC12(ISTAET,I)=HTC1  (I) 

>23  HTC22(ISTAHT,I)=HTC2(I) 

1211  HTC32(ISTABT,I)=HTC3(I) 

125  HTC42  (ISTABT, I)=HTC4  (I) 

\lt  HTCA2(ISTABT,I)=HTCA(I) 

127  HTCAS2(ISTAHT,I)=HTCAS(I) 

DTD12  (ISTABT, I)=DTD1  (I) 

129  DTD22(ISTART,I)=DTD2(I) 

130  DTD32  (ISTABT,!)  =DTD3  (I) 

131  DTD42(ISTABT,I)=DTI>4  (I) 

132  DTDA2  (ISTABT,!)  =  DTDA  (I) 

133  DTDAS2(ISTART,I)  =DTDAS  (I) 

134  DPT2  (ISTABT,!)  =DPT  (I) 

135  DPTS2(ISTART,I)  =DPTS(I) 

136  ETK2  (ISTABT, I)=XI!I  (I) 

137  -      X12  (ISTABT,  I)  =X  1(1) 

X22  <ISTART,I)=X2(I) 

139  X32(1STABT,I)=X3(I) 

1"0  X42(ISTABT,I)=X4(I) 

I1"  Q"  12(ISTABT,I)  =QW1  (I) 

11,2  QW22  (ISTABT, I)=QH2(I) 

143  QH32  (ISTABT,!)  =QW3  (I) 

144  Q»t2  (ISTABT, I)=QWH(I) 

145  Q»0  2(ISTABT,I)=Q»  (I) 

146  Q»S2  (ISTABT,I)=QWS(I) 

147  QF2  (ISTABT,!)  =QP(I) 

148  Q?S2  (ISTABT,I)=QPS(I) 

149  EBBOR2(ISTABT,I)  =ERB08  (I) 

150  1000    COHTIHUE 

'SI  ISTART=ISTABT+ 1 

152  I?  (ISTABT-5)     800,800,900 

153  900    CONTIHDE 

154  IPRIHT=0 
'55  «BITE{6,112) 

156  WRITE(6,21) 

157  3BITE(t,22) 
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158  WRITE(fc,23) 

159  DO    241     ISTART=1,5 
'6  0  II=ISTART* 100 

161  N=NCOU11T(ISTAST) 

162  DO    2tt    1=1, N 

163  III=II+I 

16a  PIN(I)=PIN2  (ISTABT,!) 

165  PINS  (I)  =PINS2  (ISTART, I) 

"*  T1(D=T0  12(ISTART,I) 

167  T1S  (I)  =T1S2  (ISTABT,!) 

,ea  GT(I)=GT2(ISTART,I) 

,69  GTS(I)=GTS2  (ISTAP.T,I) 

,7°  TC(MI)=TC02(ISTAaT,I) 

1?1  TCOS  (I)=TCOS2(ISTAHT,I) 

172  GPHB(I)  =GPt1W2  (ISTART, I) 

173  GPHWS(I)=GPMWS2  (ISTABT, I) 
I7*  IPRINT=IPRINT+1 

175  XF(IPRI!!T.E0..2t)     GO    TO     1001 

176  GO    TO    1002 

177  100  1    CONTINUE 

178  WRITE(6,112) 

179  WRITE(6,21) 

180  WHIT'?  (6, 22) 
l8'                        WRITE(6,23) 

182  IPRINT=0 

183  1002    CONTINUE 

185  241    CONTINUE 

186  IPBINT=0 

187  WRITS(6,1I2) 

188  WRITE  (6  ,26) 
'89  WRITE(6,27) 
190  WRITE  (6  ,50) 
19  1  WRITE  (6,  28) 

192  DO    291     ISTART=1,5 

'?3  H  =  ISTART»  100 

I9*  N=NCOONT(ISTAiiT) 

195  DO    29    1=1, N 

196  HI=H+I 

1  97  H  a  (I)  =  XI N2  (ISTART,  I) 

198  X1  <I)=X12  (ISTART, I) 

199  X2(I)=X22(ISTART,I) 

200  X3  <D  =  X32  (ISTART, I) 
20  1  X»  (I)  =  XQ2  (ISTART,!) 

202  TSA(I)=TSA02  (ISTART, I) 

203  TSAS(I)  =TSAS2(ISTART,I) 

204  T»»1(I)=TW»12(ISTABT,I) 
llS  TWA2(I)  =TWA  22  (ISTART,!) 
2°6  TWA3  (I)=TWA32(ISTABT,I) 

207  T*ia(I)=TWA42(ISTART,I) 

208  TWA  (I)  =TWA02  (ISTART,!) 

209  TWAS(I)  =TWAS2  (ISTART, I) 
21°  IPRINT=IPRINT+1 

2"  IT(IPRINT.EQ.26)     GO    TO     1011 

212  GO    TO    1012 

213  10  I  I    CONTINUE 
21"  WRITE(6,112) 
2'5  WRITE(6,2e) 
21  6  WRITE  (6, 27) 
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217  TOITE(6,50) 

218  WRITE(6,28) 

219  IPRIHT=0 

220  10  12    CONTIKOE 

221  "mlVm  'J2i,^'2? &■ '"  m  'n  (I>  'X3<1'  '*»«>  •*»  (I)  -"AS  (I,  , 

222  jg^SJI^W^^'^Wm.MMm.WAtD.TWlSd) 

223  IPRINT=0 

224  WRITE(6,  1  12) 

225  WRITE(6,31) 

226  WRITE(6,32) 

227  TOITE(6,33) 

228  do    341    ISTAPT=1,5 

229  II  =  ISTART*100 

230  H'RCOONT(ISTART) 

231  30    34    1=1,8 

232  III=II-H 

233  HTC1(I)=HTC12(ISTART,I) 
23ti  HTC2(I)  =HTC22(ISTART,I) 

235  HTC3  (T)  =HTC32  (ISTART,!) 

236  HTC4(I)  =HTC42(ISTART,I) 

237  HTCA(I)  =HTCA2(ISTART,I) 

238  HTCAS(I)  =  HTCAS2(ISTAaT,I) 

239  DTD  1(1)  =DTD  12(ISTART,I) 

21.0  DTD2(I)  =DTD22(ISTART,I) 

21.1  DTD3(I)=DTD32  (1ST  AST.  I) 
212  DTD4(I)  =DTD42(ISTART,I) 

21.3  DTDA  (I)  =DTDA2  (ISTAST.I) 
2"'|  DTDAS(I)=DTDAS2(ISTART,I) 
2«5  IPRIHT=TPRINT+1 

21.4  IT(IPRI)1T.EQ.26)     GO    TO     1013 
217  GO    TO    1014 

248  10  13    CONTIHDE 

219  WRITE  (6, 112) 

250  WRITE(6,3I) 

251  WRITE(6,32) 

252  WRITE(6, 33) 

253  IPRIHT=0 

254  10  14    CONTINUE 

255  3"    ™"E'P  '35)     "I.HTC1  (I)  ,HTC2  (I)  ,HTC3  (I)  ,  HTC4  (I)  ,HTCA  (I)  ,HTCAS(II 

256  34,    SmSl  m'DTI)3  (IJ  'DTDa(I)  'DTDA(I'  'DTD^(I)'  m  'HTCiS<I)  ' 

257  IPRI1IT=0 

258  SRITE(6,112) 

259  WRITE  (6  ,36) 
26  0  SPITE  (6, 37) 

261  TOITE(6,38) 

262  DO    391    ISTART=1,5 

263  II=ISTAST*100 

26  4  N=NCOUNT(ISTART) 

26  5  DO    39    1=1, H 

266  111=11*1 

2t7  ?IN(I)=PIN2(ISTAST,I) 

268  PINS  (I)  =PINS2(ISTART,I) 

269  DPT(I)=DPT2  (ISTART.I) 

270  DPTS(I)  =DPTS2(ISTART,I) 

111  I?  (DPT  (I).  GT.  0.0)     IPHIHT=IPRINT+1 

272  I?(IPPI1'T.EQ.26)     GO    TO     1015 

273  GO    TO    1016 

274  10  15    CONTINUE 
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275  WP,ITE(6,112) 

276  WRITE  (6  ,36) 

277  WRITE  (6,  37) 

278  WRITE  (6  ,38) 

279  IPHIHT=0 

280  1016    CONTINUE 

2882  3"  Sktihoe'""'0-0'  mTE(fc'40'  vxi*nm.n»m.mm.wam 

283  IPRINT=0 

281  WPITE(6,112) 

285  WRITE  (6  ,4  1) 

286  WRITE(6,42) 

287  WRITE(6,43) 

288  DO     441    ISTABT=I,5 

289  II  =  ISTART*100 

290  N=NCOOTT(ISTART) 

291  DO    44   1  =  1, N 

292  111=11  +  1 

293  Q"'(D=QW12(ISTAET,I) 
25"  QW2(I)=Q.W22(ISTABT,I) 

295  Q"  (I)=QW32  (ISTABT.I) 

296  QW«(ri=QW42(ISTAST,I) 

297  O.H(I)=QW02(ISTABT,I) 

298  Q"S(D=Q¥S2(ISTART,I) 

299  QF  (I)=QF2  (ISTAST.I) 

300  QFS(I)  =QFS2(ISTABT,I) 

301  EBBOB(I)=EHBOB2  (ISTAST.I) 

302  IPRIIIT=IPRIHT+ I 

303  IF  (IPHINT.EQ.26)     GO    TO    1017 
301  so    TO     1018 

305  1017    CONTINUE 

-306  WBITE(6,  1  12) 

3<>7  WRITE  (6,41) 

308  WPITE(6,42) 

3"9  WRITE(6,43) 

3  10  IPRINT=0 

311  1018    CONTINUE 

iSSfffi^,m8;8.SS&alB1'Bfla,'fl,"w'0,,m'w«'Bf« 

314  44    CONTINUE 

3  15  441    C3NTIH0E 

316  WBITE(6,112) 

317  ST0p 

313  END 
SENTRY 
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ooincx)^inuDcj^(N^cN=rcor~u3fv)cri^u3^>cofNc^<N^r 
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confnh^^oicNnoocooj^i^crtMNcccNa^tapncrim 
*-^is-(^mo4in  — Mr-  ^ruii-cNm  —  ^-cncn- r^vu^-©,«oiA 

n*-ocDi?a,oomku^o^inin(ri3,ocoa,'^inyDn(N(^ 
Ln^r*co^i7i'-fN!N3ijiu^fva3oomujr'(k*cricTi-i-fri 
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COOd,*Of*H'-<7'OMhfnintfh\o^^>-rnoi?-m- 

«-  —  <-  —  *-*-*-—  t-  —  t-  —  *-—CN-  «-  —  *-  —  i-  —  CN  — 

o^pn^^iriinoD^oja3CNO(^r^^c,n^cT''CocNOcoLO^rco 
©p-mLnmLficoLOir.  su  ^  ^  co  rn  no  —  v^omcoor*o*~'u".  cn 

tf^r^OC7^OCNaDC0C^©^\yt^Or^^v£C7*C0—  CTiOt-**JCOC"> 
<-CNr-3-fninojCTi»£l^c^iri3'C0cDCNr-   CN  (Ji   K   ^  O   r  <N*"^T 

p^r^^*ijr^oaD^^C7icrr^^oLn^r^.^Ln<r*  —  cncolt)cnct>co 
in=3-^-r-r*-  —  crcoyDrsiin  —  oaoco^'^-or^cocN  —  ^r-cNro 

U3^or*U3-  c3-CNOCDmO^CO^CN--rOCT\ira'COOC7>r*~OiT 

r-«i>^--  n-£iT-oa30comr,*^inuimnvocornha'>',oc^c 
m  —  *—  trivo>urocDrvcD3-rfc^rino(N(rio^a'>'^r'fricritNr'to 

CN^a^LOCO^r^-   COtNr-O^ySrOcNr^t—   0~ir-CDCN^^OfNO 
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moo^r-coLDomcNr-LnLnom^'cN  —  r-r-^rt^^otN^LTj 
mfNu?cnococD(r>in  —  r-a-cn  —  a-o^or^o  —  *-roCTiomr- 

CNtNcTirocor-co^co  —  m  —  Lftsga3SOt0O*Qorh***4?<7v0tcofN 

c->  cn  o  —  crit*-a)—  oiKD^nOfN-  iflaicnotNvo^-  *-  a- 
nnrna-^^cocDJT^voincNCDaitNfMCNr-  i--  =r  cr*  vo  hj  c<  u"> 

<-  —  <-  —  *-  —  <-  —   r-  —  *-  —  t-  —   i-CN<N  —  <--  t-—   <-  —  ,-fN 
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APPENDIX  E 

UNCERTAINITY  ANALYSIS  IN  EXPERIMENTAL  MEASUREMENTS 
OF  OVERALL  HEAT  TRANSFER  COEFFICIENTS 

Due  Co  the  limits  in  the  accuracy  with  which  the  measurements  were 
taken,  the  final  result  of  intrest  is  associated  with  certain  uncertainity . 
The  attempt  was  made  to  determine  the  uncertainity  in  the  experimental 
measurements  of  the  condensation  heat  transfer  coefficients.  The  method 
outlined  by  Kline  and  McClintock  f  56  J   was  used  in  estimating  the 
uncertainity  associated  with  overall  heat  transfer  coefficient  h  for  the 
smooth  tube,  run  144. 

The  overall  average  heat  transfer  coefficient  is  calculated  by 

h=  9 (E-l) 

A  ( T  -  T  . ) 
n  s   wi 

The  experimental  uncertainity  for  the  overall  average  heat  transfer 

coefficient  is  given  by 


,Bh   .2  ,  .ah      ,2    .Bh     .2 

L  n    n        s   s 

1/2 


T 

'W.  WT.) 

Wl           Wl 

Eq 

(E-l) 

*h 

1 

*"        A 
n 

( T   -   T    . ) 

S           Wl 

3h 

c-A 

n 

Q 

A2(T  -  f   .) 
nx   s        wi 

Vh_  . 

Q 

(E-2) 


3t      A  (T  -  T  .)2 
s      n  s   wi 


(E-3) 


(E-4) 


(E-5) 
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oh m  0 

oT  .   A  (T  -  T  .)2  (E_6) 

»i    n  s   wi 

For  this  run 

A  -1TD.L  =   (%|^-)(9.0)  =  1.284  ft2 
n     l       Iz 

T  =  103.16  °F 

s 

From  the  one  dimensional  heat  conduction  equation,  the  inside  wall 
temperature  can  be  estimated  as 

\±  -  \o +  rii ln  «5?  (E-7) 

i 

Also,  neglecting  the  temperature  drop  in  the  tube  wall,  it  can  be 
approximated  as 

T  .  =  T   =  87.08  °F 
wi    wo 

For  this  run,  the  heat  gained  by  the  cooling  water 

Q  =  5491.16  Btu/hr 

Substituting  in  Eqs.(E-3)  to  (E-6)  we  get 

oh       1  1 


5Q        1.284(103.16-87.08)        20.67 


(E-8) 


g_. i542Llfi =_207.12  (e-9) 

n        (1.284)    (103.16-87.08) 

3h -5491.16 

— — j  =  -16.52  (c-10) 

oT  (1.284)(103. 16-87. 08)" 

a£_  . 5491^ _  16_52  (E_U) 

3T  ,         (1.284)(103. 16-87. 08)" 
wl 

The   estimation  of    the  uncertainity  of   energy   transfer,   w    : 
The   energy   transfer   is   given  by 

Q  =  m       Cp        (T       -  T      .)  (E-12) 

wa       wa       wao        wai 


The  uncertainity  in  Q  is  given  by 


h 


m    m 
.*-  wa   wa 
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I  3Cp    Cp 
^    wa    waj 


22. 


ot        t 

^  wao   wao 


*L 


1  2 


-,1/2 


oT   .  "T 
L.  wai   waiJ  _ 


+  heat  exchange  with  the 
environment 
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(E-13) 


For  this  particular  run 


m   =  554.66  lbm/hr 
wa 

Cp    =1.0  Btu/lbm.°F 
wa 

T    =  76.64  °F 
wao 

T   .  =  67.28  °F 

wai 

Therefore 


f*—  =  Cp   (T   -T   .)  =  1.0(76.64-67.28) 
dm      wa  wao   wai 
wa 


=9.36  Btu/lbn 


|5 =  ii   (T   -  T  J  =  554.66(76.64-67.28) 

dCp      wa  wao   wai 
wa 

=  5191.63  lbm.°F/hr 


|2 -  i  Cp   =  (554.66)(1.0) 

dT       wa  wa 
wao 

=  554.66  Btu/hr.°F 


-2 _  A   Cp   --(554.66X1.0) 

dT        wa   wa 
wai 

=  -554.66  Btu/hr.°F 

The  uncertainity  associated  with  the  water  flow  was  assumed  to  be 
within  +  1%  of  the  flow  rate.  Therefore, 


w.   =  0.01  x  554.66  =  5.54  lbm/hr 
m 
wa 
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and  the  uncertainity  of  the  specific  heat  w„    was  assumed  to  be 

Cp 

0  wa 

w„    =  0.004  Btu/lbm.  F 
Cp 
wa 

The  uncertainities  in  T   ,  and  T   .  were  determined  from  the 
wao       wai 

uncertainities  of  thermocouples  and  data  acquisition  system  reading 
as  follows  : 

Uncertainity  due  to  thermocouple  wire  inaccuracies  +  0.75  F 
Uncertainity  due  to  data  acquisition  system +.  0.01  %  of  the  reading 


Therefore, 


1/2 


wT    =  [(0.75)   +  (76.64  x  0.0001)  1    =  0.75  °F 


and 


[(0.75)2  +  (67.28  x  0.0001)2 


,1/2 


»  0.75  F 


For  this  particular  run,  heat  balance  error  of  -0.3%  obtained  and 
it  was  considered  to  be  the  estimate  of  the  heat  transferred  to  the 
ambient.  Substituting  the  values  in  Eq.(E-13)  yields 

w  =  [(9.36  x  5.54)2  +  (5191.63  x  0.004)2  +  (554.66  x  0.75)2 
4  1/2 

(-554.66  x  0.75)2]    +  (0.003  x  5491.16) 

=  607.43  Btu/hr 


The  estimation  of  the  uncertainity  of  heat  transfer  area,  w 


A 


A     =1TD.L 


(E-14) 


Thus,    the  uncertainity 

f3A              I2        f3A         " 
n  n 

357   V     +   3i7wl 


1/2 


2  2 

(tfLw      )      +    (irD.w   ) 


nl/2 


(E-15) 


For  this  particular  run, 


D  =  0.0454  ft 
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L  =  4Az  ■  9.0  ft 
The  uncertainities  in  measuring  the  inside  diameter  and  the  length  of 
the  tube  were  estimated  by 

w   =  0.001  ft 

Di 

»L  -  0.01  ft 
Therefore, 

r  2  2-1 1/2 

WA  =  L(TrX  9'°  X  °-001)   +  (Trx  0-0*54  x  0.01)  J 
=  0.028  ft2 
The  estimation  of  the  uncertainity  for  the  average  saturation 

temperature,  w_  : 

T 
s 

The  average  saturation  temperature  is  a  function  of  the  saturation 

pressure  P 
s 

VVV  (E-16) 

The  uncertainity  is  estimated  by 

pT 

Ts   L  s     sj 

The  uncertainity,  w   ,  results  from  the  following  uncertainities 
s 

Uncertainity  due  to  the  pressure  gage  reading:  +0.2  psia 
Uncertainity  due  to  the  pressure  gage  resolution:  +0.2  psia 
Uncertainity  associated  with  the  pressure  fluctuations  in  the 
system:  7.  0.25  psia 

Therefore, 

1/2 
wp  =  [(0.2)2  +  (0.2)2  +(0.25)2]    =  0.377  psia 
s 

From  the  thermophysical  property  tables  given  in  ASHRAE  Handbook  of 
Fundamentals,  and  at  the  pressure  range  of  interest  to  this  experimental 
run,  it  was  estimated  that 
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At 

d?     '    AP        2-27      F/psla 
s 

Therfore, 

w_     =["(2.27   x  0.377)]=   0.856 

T 
s 

The  estimation  of  the  uncertainity  of  the  average  inside  wall 

temperature,  w_  : 

T  . 
wi 


As  shown  earliear  in  this  Appendix 

T  .  =  I 
wi    wo 

The  factors  which  contribute  to  the  uncertainity  of  the  average 

outside  wall  temperature,  were 

Uncertainity  due  to  thermocouple  wire,  inaccuracies  +0.75  F 
Uncertainity  associated  with  the  attatchment  method  Z  1.5  f 
Uncertainity  due  to  the  data  acquisition  system    +  0.01%  of  the 

reading 

Therefore, 

1/2 
w    S  w    =  [(0.75)2  +  (1.5)2  +  (87.08  x  0.0001)2"] 


T  .    T 
wi     wo 


=  1.67  F 


The  values  of  w  ,  w   ,  w_  ,  and  w_   are  substituted  in  Eq.(E-2) 

4    n   T        T  . 
s        wi 


to  yield 


w_  =  R2O67  X  607-43)2  +  (-207.12  x  0.028)2  +  (-16.52  x  0.856)2 
h  1/2 

+  (16.52  x  1.67)2"] 

=  +  43.10  Btu/hr.ft2.°F 

Thus  for  this  run, 


h  =  266.87  +  43.10     Btu/f t2 .hr .°F 
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Therefore,  the  uncertainity  for  the  experimental  condensation  overall 

average  heat  transfer  coefficient  of  this  run  is  about  +  16.15%.  It  can 
also  be  assumed  that  this  uncertainity  holds  true  for  all  other  runs. 

The  main  reason  for  such  a  high  value  of  uncertainity  can  be  attributed 

to  the  fact  that  the  driving  temperature  difference  (T  -  T  ,)  was  low 

s    wi 

for  this  run. 
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APPENDIX  F 


Table  F.l:  Tabulation  of  Stastistical  Information  for  the  Curve 

Fits  to  the  Experimental  Overall  Heat  Transfer  Coefficient 

8 

h  =  C  G 


3  0.372 

4  0.320 

5  0.519 


1.47                             1  0.736 

2  0.675 

3  0.614 

4  1.015 

5  0.666 


1.67  1  0.531 

2  0.7U 


3  0.330 

4  1.031 

5  0.360 


Nominal  Tube  Correlation 

Inlet  No.  Coefficient 

Pressure  n  C 

(bar) 


1.32  1  0.654 

2  0.507 


0.281 

0.934 

2.197 

0.908 

6.465 

0.664 

6.233 

& 

-02 

0.924 

1.681 

0.870 

0.104 

0.972 

0.309 

0.976 

0.420 

0.906 

7.638 

E- 

-03 

0.929 

0.361 

0.796 

0.905 

0.938 

0.198 

0.961 

3.420 

0.673 

6.172 

E- 

■03 

0.935 

4.964 

E-02 

0.942 
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Table  F.2:  Tabulation  of  Scastical  Information  for  the  Curve 
Fits  to  the  Experimental  Overall  Pressure  Drop 


C  Ga 


Nominal 

Tube 

Correlation 

Inlet 

So. 

Coefficient 

Pressure 

n 

C 

(bar) 

1.32 

1 

1.200 

1.617 

E-ll 

0.759 

2 

1.052 

1.199 

E-06 

0.985 

3 

1.253 

2.440 

E-07 

0.973 

1.67 

1 

1.357 

2.668 

E-08 

0.809 

2 

1.700 

1.184 

E-09 

0.986 

3 

1.242      ■ 

2.671 

E-07 

0.973 
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AUGMENTATION  OF  CONDENSATION  HEAT  TRANSFER  OF  R-U 
BY  INTERNALLY  FINNED  TUBES 

ABSTRACT 

The  continuous  demand  for  more  efficient  heat  exchangers  has 
stimulated  the  interest  of  investigators  to  explore  the  use  of  different 
techniques  to  augment  the  heat  transfer.  Before  incorporating  any  of  the 
augmentation  techniques  into  a  full  scale  industrial  use,  a  considerable 
amount  of  testing  is  necessary.  The  ultimate  objective  is  the  development 
of  reliable  design  correlations  for  predicting  the  heat  transfer  and 
pressure  drop  in  heat  exchangers  with  heat  transfer  augmented  surfaces. 
The  present  study  was  concerned  with  augmentation  of  in-tube  condensation 
of  R-ll  by  internally  finned  tubes. 

Heat  transfer  and  pressure  drop  data  were  taken  during  condensation 
of  R-ll  inside  five  horizontal  tubes.  One  tube  was  smooth,  one  had 
straight  internal  fins  and  the  remaining  three  had  spiral  internal  fins 
with  different  outside  diameters  and  fin  geometries.  The  smooth  tube 
was  used  as  a  basis  for  comparison  with  the  results  of  the  finned  tubes. 
The  experimental  results  of  all  tubes  tested  were  compared  with  predictions 
of  numerous  existing  correlations  in  order  to  identify  the  correlation 
which  best  predicted  the  results  of  the  present  study. 

The  results  are  summarized  below. 

1.  On  a  nominal  area  basis,  an  enhancement  in  heat  transfer  as  high  as 
55%  over  the  smooth  tube  results  was  obtained  by  internally  finned  tubes, 
over  the  mass  flux  range  tested.  The  heat  transfer  was  accompanied  by  a 
modest  increase  in  the  pressure  drop. 

2.  The  best  predictor  of  the  smooth  tube  sectional  and  overall  heat 
transfer  coefficients  was  the  correlation  of  Akers  et  al . [6J .   For  internally 


finned  tubes,  the  correlation  of  Luu  and  Bergles  /"45,467  predicted  the 
sectional  and  overall  heat  transfer  coefficient  best. 

3.  The  best  predictor  of  the  pressure  drop  of  the  smooth  tube  was  the 
combination  of  the  frictional  drop  correlation  of  Lockhart-Martinelli  [\kj 
and  the  void  fraction  correlation  of  Zivi  [2\] .   Good  correlation  was  also 
obtained  by  replacing  the  correlation  of  Lockhart-Martinelli  by  the 
correlation  of  Dukler  [li] ,    for  the  frictional  pressure  drop. 

4.  The  same  combinations  of  frictional  pressure  gradient  and  void  fraction 
correlations  of  the  smooth  tube  predicted  the  finned  tubes  pressure  drop 

by  replacing  the  nominal  inside  diameter  by  the  hydraulic  diameter. 

5.  The  ratio  of  the  pumping  power  per  unit  heat  transfer,  subject  to  the 
constraints  of  fixed  geometry,  same  inlet  pressure  and  flow  rate  was  used  to 
evaluate  the  performance  of  one  of  the  finned  tubes  in  relation  to  the 
smooth  tube.  The  results  showed  that  over  the  mass  flow  rate  tested,  the 
smooth  tube  required  less  pumping  power  per  unit  heat  transfer. 

6.  The  ratio  of  heat  transfer  enhancement  of  the  finned  tube  to  smooth 
tube  to  a  similar  ratio  for  the  pressure  drop  was  also  used  to  evaluate 
the  performance  of  the  same  finned  tube.  The  ratio  was  also  subject  to  the 
same  constraints  as  the  ratio  of  the  pumping  power  per  unit  heat  transfer. 
This  ratio  was  lower  than  1  but  increased  with  the  increase  in  mass  flow 
rate.  The  result  suggested  that  this  ratio  could  exceed  1  at  higher  mass  flow 
rates.  However,  such  observation  could  not  be  confirmed  because  of  the 
limitation  of  the  test  facility.  A  ratio  greater  than  one  indicates  that 

the  increase  in  heat  transfer  of  the  internally  finned  tube  is  accompanied 
only  by  a  modest  increase  in  pressure  drop.  Hence,  in  the  flow  range  where 
this  ratio  could  be  greater  than  one,  internally  finned  tubes  could  be 
advantageous  compared  to  smooth  tubes. 


